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We present a theory of a degenerate atomic Fermi gas, interacting through a narrow Feshbach 
resonance, whose position and therefore strength can be tuned experimentally, as demonstrated 
recently in ultracold trapped atomic gases. The distinguishing feature of the theory is that its 
accuracy is controlled by a dimensionless parameter proportional to the ratio of the width of the 
resonance to Fermi energy. The theory is therefore quantitatively accurate for a narrow Feshbach 
resonance. In the case of a narrow s-wave resonance, our analysis leads to a quantitative description 
of the crossover between a weakly-paired BCS superconductor of overlapping Cooper pairs and a 
strongly-paired molecular Bose-Einstein condensate of diatomic molecules. In the case of pairing 
via a p-wave resonance, that we show is always narrow for a sufBciently low density, we predict a 
detuning-temperature phase diagram, that in the course of a BCS-BEC crossover can exhibit a host 
of thermodynamically-distinct phases separated by quantum and classical phase transitions. For an 
intermediate strength of the dipolar anisotropy, the system exhibits a. -\- ipy paired superfluidity 
that undergoes a topological phase transition between a weakly-coupled gapless ground state at large 
positive detuning and a strongly-paired fully-gapped molecular superfluid for a negative detuning. 
In two dimensions the former state is characterized by a Pfafhan ground state exhibiting topological 
order and non-Abelian vortex excitations familiar from fractional quantum Hall systems. 

PACS numbers: 03.75.Ss, 67.90.+Z, 74.20.Rp 
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I. INTRODUCTION 

A. Weakly- and strongly-paired fermionic 
superfluids 

Paired superfluidity in Fermi systems is a rich sub- 
ject with a long history dating back to the discovery 
of superconductivity (charged superfluidity) in mercury 
by Kamerlingh Onnes in 1911. Despite considerable 
progress on the phenomenological level and many exper- 
imental realizations in other metals that followed, a de- 
tailed microscopic explanation of superconductivity had 
to await seminal breakthrough by Bardeen, Cooper and 
Schrieffer (BCS) (for the history of the subject see, for 
example, Ref. [1] and references therein). They discov- 
ered that in a degenerate, flnite density system, an ar- 
bitrarily weak fermion attraction destabilizes the Fermi 
sea (primarily in a narrow shell around the Fermi en- 
ergy) to a coherent state of strongly ovcrlaping "Cooper 
pairs" composed of weakly correlated time-reversed pairs 
of fermions. 

In contrast, superfluidity in systems (e.g., liquid ^He), 
where constituent fermions (neutrons, protons, electrons) 
are strongly bound into a nearly point-like bosonic atom, 
was readily qualitatively identifled with the strongly in- 
teracting liquid limit of the Bose-Einstein condensation 



of composite bosonic *He atoms (for a review, see for 
example Ref. [2]). 

While such weakly- and strongly-paired fermionic s- 

wave superfluids were well understood by early 1960's, 
the relation between them and a quantitative treatment 
of the latter remained unclear until Eagles's [3] and later 
Lcggctt's [4], and Nozicrcs and Schmitt-Rink's [5] sem- 
inal works. Working with the mean-field BCS model, 
that is quantitatively valid only for a weak attraction 
and high density (a superconducting gap much smaller 
than the Fermi energy), they boldly applied the model 
outside its quantitative range of validity [6] to fermions 
with an arbitrarily strong attraction. Effectively treat- 
ing the BCS state as a variational ground state, such 
approach connected in a concrete mean-field model the 
two types of s-wave paired superfluids, explicitly demon- 
strating that they are two extreme regimes of the same 
phenomenon, connected by a smooth (analytic) crossover 
as the strength of attractive interaction is varied from 
weak to strong. This lack of qualitative distinction be- 
tween a "metaUic" (BCS) and "molecular" (BEC) s-wave 
superfluids, both of which are characterized by a com- 
plex scalar (bosonic) order parameter 5*, was also antic- 
ipated much earlier based on symmetry grounds by the 
Ginzburg-Landau theory [1]. 

Nevertheless, the two types of superfluids regimes 
exhibit drastically (quantitatively [7]) distinct phe- 
nomenologies [5, 8]. While in a weakly-paired BCS su- 
perconductor the transition temperature Tc nearly co- 
incides with the Cooper-pair binding (dissociation) en- 
ergy, that is exponentially small in the pairing poten- 
tial, in the strongly-paired BEC superfluid Tc is deter- 
mined by the density, set by the Fermi temperature, and 
is nearly independent of the attractive interaction be- 
tween fermions. In such strongly coupled systems the 
binding energy, setting the temperature scale T* above 
which the composite boson dissociates into its constituent 
fermions (e.g., of order eV in ^He) can therefore be or- 
ders of magnitude larger than the actually condensation 
temperature Tc <C . This large separation between Tc 
and T« is reminiscent of the phenomenology observed in 
the high-temperature superconductors (with the range 
Tc < T < T^, referred to as the "pseudo-gap" regime), 
rekindling interest in the BCS-BEC crossover in the mid- 
90's [8] and more recently [9]. 

With a discovery of novel superconducting materials 
(e.g., high-Tc's, heavy fermion compounds), and super- 
fluids (^He), that are believed to exhibit finite angular 
momentum pairing, the nature of strongly- and weakly- 
paired superfluids has received even more attention. It 
was soon appreciated [10-12] that, in contrast to the 
s-wave case, strongly and weakly paired states at a fi- 
nite angular momentum are qualitatively distinct. This 
is most vividly illustrated in three dimensions, where 
for weak attraction a two-particle bound state is absent, 
the pairing is stabilized by a Fermi surface and there- 
fore necessarily exhibits nodes and gapless excitations 
in the finite angular momentum paired state. In con- 



trast, for strong attraction a two-particle bound state 
appears, thereby exhibiting a fully-gapped superfluid- 
ity with concomitant drastically distinct low tempera- 
ture thermodynamics. Other, more subtle topological 
distinctions, akin to quantum Hall states, between the 
two types of paired grounds states also exist and have 
been investigated [10, 11]. Consequently, these qualita- 
tive distinctions require a genuine quantum phase transi- 
tion (rather than an analytic crossover, as in the case of 
s-wave superfluid) to separate the weakly and strongly- 
paired states. This transition should be accessible if the 
pairing strength were experimentally tunable. 

B. Paired superfluidity via a Feshbach resonance 

The interest in paired superfluidity was recently re- 
vived by the experimental success in producing degen- 
erate (temperature well below Fermi energy) trapped 
atomic Fermi gases of^Liand^^K [13 16]. A remarkable 
new experimental ingredient is that the atomic two-body 
interactions in these systems can be tuned by an external 
magnetic field to be dominated by the so-called Feshbach 
resonant (FR) [17, 18] scattering through an intermedi- 
ate molecular (virtual or real bound) state. 

As depicted in Fig. 1, such tunable Feshbach reso- 
nance [19] arises in a system where the interaction poten- 
tial between two atoms depends on their total electron 
spin state, admitting a bound state in one spin chan- 
nel (usually referred to as the "closed channel" , typically 
an approximate electron spin-singlet). The interaction 
in the second [20] "open" channel (usually electron spin 
triplet, that is too shallow to admit a bound state) is then 
dominated by a scattering through this closed channel 
resonance [21]. Since the two channels, coupled by the 
hyperfine interaction, generically have distinct magnetic 
moments, their relative energies (the position of the Fes- 
hbach resonance) and therefore the open-channel atomic 
interaction can be tuned via an external magnetic field 
through the Zeeman splitting, as depicted in Fig. 1. 

In the dilute, two-body limit the low-energy s-wave 
Feshbach resonant scattering is characterized by an ,s- 
wave scattering length, that, as illustrated in Fig. 2, is 
observed to behave according to [22, 23] 

a(ff)=a6,(l--^^), (1.1) 

diverging as the magnetic field passes through a (system- 
dependent) field Hq, corresponding to a tuning of the 
resonance through zero energy. (Analogously, a p-w&ve 
resonance is characterized by a scattering volume v{H), 
as discussed in detail in Sec. II B). In above, the ex- 
perimentally measurable parameters atg and Hyj are, re- 
spectively, the background (far off-resonance) scattering 
length and the so-called (somewhat colloquially; see [24]) 
magnetic "resonance width" . 

An ,s-wavc Feshbach resonance is also characterized by 
an additional length scale, the so-called effective range 
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FIG. 1: Interactions between tlie atoms generically depends 
on their mutual spin state. This figure depicts two potentials 
corresponding to two spin states of the pairs of atoms. One of 
them (usually referred to as an "open channel" ) is too weak to 
support a bound state, while the other (a "closed channel") 
supports a bound state, but is energetically unfavorable at 
large distances. The closed channel potential can bo moved 
vertically with respect to the open channel potential via the 
Zeeman effect, by changing an external magnetic field. 

ro, and a corresponding energy scale 

ro = — , (1.2) 

mrQ 

that only weakly depend on H. This important scale 
measures the intrinsic energy width of the two-body reso- 
nance and is related to the measured magnetic-field width 

H^^, via 

Fo « AmfilalgHl/h\ (1.3) 

with the Bohr magneton. Fq sets an energy crossover 
scale between two regimes of (low- and intermediate- 
energy) behavior of two-atom s-wave scattering ampli- 
tude. 

A key observation is that, independent of the nature of 
the complicated atomic interaction leading to a Feshbach 
resonance, its resonant physics outside of the short mi- 
croscopic (molecular size of the closed-channel) scale can 
be correctly captured by a pseudo-potential with an iden- 
tical low-energy two-body scattering amplitude, that, for 
example, can be modeled by a far simpler potential ex- 
hibiting a minimum separated by a large barrier, as il- 
lustrated in Fig. 5. The large barrier suppresses the de- 
cay rate of the molecular quasi-bound state inside the 
well, guaranteeing its long lifetime even when its energy 
is tuned above the bottom of the continuum of states. 

Although such potential scattering. Fig. 5 is microscop- 
ically quite distinct from the Feshbach resonance. Fig. 1, 
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FIG. 2: A schematic of a typical, experimentally observed 
behavior of an s-wavc scattering length a{H) as a function of 
magnetic field H in & vicinity of a Feshbach resonance. 



this distinction only appears at high energies. As we 
will see, the low energy physics of a shallow resonance 
is controlled by a nearly universal scattering amplitude, 
that depends only weakly on the microscopic origin of 
the resonance. Loosely speaking, for a large barrier of 
a potential scattering depicted on Fig. 5 one can asso- 
ciate (quasi-) bound state inside the well with the closed 
molecular channel, the outside scattering states with the 
open channel, and the barrier height with the hyperfine 
interactions-driven hybridization of the open and closed 
channels of the Feshbach resonant system. The appro- 
priate theoretical model was first suggested in Ref. [19], 
and in turn exhibits two-body physics identical to that 
of the famous Fano- Anderson model [25] of a single level 
inside a continuum (see Appendix B). 

A proximity to a Feshbach resonance allows a high 
tunability (possible even in "real" time) of attractive 
atomic interactions in these Feshbach-resonant systems, 
through a resonant control of the s-wave scattering length 
a{H), Eq. (1.1) via a magnetic field. As we will discuss 
in Sec. VII, a p-wave Feshbach resonance similarly per- 
mits studies of p-wave interacting systems with the in- 
teraction tunable via a resonant behavior of the scatter- 
ing volume v{H). This thus enables studies of paired 
superfluids across the full crossover between the BCS 
regime of weakly-paired, strongly overlapping Cooper 
pairs, and the EEC regime of tightly bound (closed- 
channel), weakly interacting diatomic molecules. More 
broadly, it allows access to interacting atomic many- 
body systems in previously unavailable highly coherent 
and even noncquilibrium regimes [26-28], unimaginable 
in more traditional solid state systems. 



C. Narrow vs wide resonances and model's validity 

An atomic gas at a finite density n (of interest to us 
here) provides an additional length, n~^/^ ~ and 
corresponding energy, ep = h'^kp/'^fm scales. For the s- 
wave resonance, these scales, when combined with the 
length vq or the resonance width Fo, respectively, allow 
us to define an s-wave dimensionless parameter (with nu- 
merical factor chosen for later convenience) 

..^^JH^IJl. (1.4) 

that measures the width of the resonance or equivalcntly 
the strength of the Feshbach resonance coupling (hy- 
bridization of an atom-pair with a molecule) relative to 
Fermi energy. For a p-wave (and higher angular momen- 
tum) resonance a similar dimensionless parameter can be 
defined (see below) . The key resonance- width parameter 
7 [24] naturally allows a distinction between two types 
of finite density Feshbach-resonant behaviors, a narrow 
(7 <C 1) and broad (7 ^ 1). Physically, these are dis- 
tinguished by how the width Fq compares with a typical 
atomic kinetic energy ep. Equivalcntly, they are con- 
trasted by whether upon growth near the resonance, the 
scattering length a{H) first reaches the effective range 
|ro| (broad resonance) or the atom spacing i (narrow res- 
onance). 

Systems exhibiting a narrow resonant pairing are ex- 
tremely attractive from the theoretical point of view. As 
was first emphasized in Ref. 28 and detailed in this pa- 
per, such systems can b(^ accurately modeled by a sim- 
ple two-channel Hamiltonian characterized by the small 
dimensionless parameter 7, that remains small (corre- 
sponding to long-lived molecules) throughout the BCS- 
BEC crossover. Hence, while nontrivial and strongly 
interacting, narrow Feshbach resonant systems allow a 
quantitative analytical description, detailed below, that 
can be made arbitrarily accurate (exact in the zero res- 
onance width limit), with corrections controlled by pow- 
ers of the small dimensionless parameter 7, computable 
through a systematic perturbation theory in 7. The 
ability to treat narrowly resonant systems pcrturbatively 
physically stems from the fact that such an interaction, 
although arbitrarily strong at a particular energy, is con- 
fined only to a narrow energy window around a resonant 
energy. 

As we will show in this paper [28], such narrow res- 
onant systems exhibit a following simple picture of a 
pairing superfiuid across the BCS-BEC crossover, illus- 
trated in Fig. 3. For a Feshbach resonance tuned to a 
positive (detuning) energy the closed-channel state is a 
resonance [29] , that generically leads to a negative scat- 
tering length and an effective attraction between two 
atoms in the open-channel. For detuning larger than 
twice the Fermi energy, most of the atoms are in the 
open-channel, forming a weakly BCS-paired Fermi sea, 
with exponentially small molecular density, induced by a 
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weak Feshbach resonant (2-)atom-niolecule coupling (hy- 
bridization). The BCS-BEC crossover initiates as the 
detuning is lowered below 2ep, where a finite density 
of atoms binds into Bose-condensed (at T — {)) closed- 
channel quasi-molecules, stabilized by the Pauli principle. 
The formed molecular (closed-channel) superfluid coex- 
ists with the strongly-coupled BCS superfluid of (open- 
channel) Cooper pairs, that, while symmetry-identical 
and hybridized with it by the Feshbach resonant coupling 
is physically distinct from it. This is made particularly 
vivid in highly anisotropic, one dimensional traps, where 
the two distinct (molecular and Cooper-pair) superflu- 
ids can actually decouple due to quantum fluctuations 
suppressing the Feshbach coupling at low energies [30]. 
The crossover to BEC superfluid terminates around zero 
detuning, where conversion of open-channel atoms (form- 
ing Cooper pairs) into closed-channel molecules is nearly 
complete. In the asymptotic regime of a negative de- 
tuning a true bound state appears in the closed-channel, 
leading to a positive scattering length and a two-body 
repulsion in the open-channel. In between, as the po- 
sition of the Feshbach resonance is tuned through zero 
energy, the system is taken through (what would at 
zero density be) a strong unitary scattering limit, cor- 
responding to a divergent scattering length, that is nev- 
ertheless quantitatively accessible in a narrow resonance 
limit, where 7<( ~ l/ikpro) plays the role of a small 
parameter. This contrasts strongly with systems inter- 
acting via a featureless attractive (e.g., short-range two- 
body) potential, where due to a lack of a large poten- 
tial barrier (see Fig. 12) no well-defined (long-lived) res- 
onant state exists at positive energy and a parameter 
7s (proportional to the inverse of effective range tq) is 
effectively infinite. For such broad-resonance systems, 
a gas parameter v}/^ \a[H)\ is the only dimensionless 
parameter. Although for a dilute gas {n^/^abg <C 1) a 
controlled, perturbative analysis (in a gas parameter) 
of such systems is possible away from the resonance, 
where v}/^ \a{H)\ ^ 1, a description of the gas (no 
matter how dilute), sufficiently close to the resonance, 
such that ■n}/'^\a{H)\ > 1 is quantitatively intractable in 
broad-resonance systems [31]. This important distinc- 
tion between the narrow and broad Feshbach resonances 
and corresponding perturbatively-(in)accessible regions 



in the kp 



plane around a Feshbach resonance are 



illustrated in Fig. 4. 

Nevertheless, because of their deceiving simplicity and 
experimental motivation (most current experimental sys- 
tems are broad), these broad- resonance systems (exhibit- 
ing no long-lived positive energy resonance [29]) were a 
focus of the aforementioned earlier studies [3-5, 8] that 
provided a valuable qualitative elucidation of the BCS- 
BEC crossover into the strongly-paired BEC supcrflu- 
ids. However, (recent refinements, employing enlight- 
ening but uncontrolled approximations notwithstand- 
ing [9, 19, 32, 33]) these embellished mean- field descrip- 
tions are quantitatively untrustworthy outside of the BCS 
regime, where weak interaction (relative to the Fermi en- 
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FIG. 3: An illustration of the BCS-BEC crossover in the limit 
of a narrow Feshbach resonance width, 7^ ^ 1. The evolution 
with detuning loo is illustrated, with (a) the BCS regime of 
uoo > 2eF, where particles are predominantly open-channel 
atoms forming a Cooper-paired Fermi sea, (b) the crossover 
regime of < wo < 2eF, where a fraction of atoms between ujo 
and ep have converted into a BEC of bosonic (closed-channel) 
molecules, with the rest forming a Cooper-paired Fermi sea 
at a chemical potential ^, and (c) the BEC regime of uoo < 0, 
where (to order 7s ^ 1) only Bose-condensed molecules are 
present. 



ergy) provides a small parameter justifying a mean-field 
treatment - and outside of the BEC regime where, al- 
though mean-field techniques break down, a treatment 
perturbative in rt^/^|a| <C 1 is still possible [6]. The in- 
ability to quantitatively treat the crossover regime for 
generic (non-resonant) interactions is not an uncommon 
situation in physics, where quantitative analysis of the 
intermediate coupling regime requires an exact or numer- 
ical solution [34]. By integrating out the virtual molec- 
ular state, systems interacting through a broad (large 7) 
resonance can be reduced to a nonresonant two-body in- 
teraction of effectively infinite 7, and are therefore, not 
surprisingly, also do not allow a quantitatively accurate 
perturbative analysis outside of the BCS weak-coupling 
regime [31]. 

The study of a fermionic gas interacting via a broad 
resonance reveals the following results. If a < (the in- 
teractions are attractive but too weak to support a bound 



state) and n 



1/3 I 



<C 1, such a superfluid is the stan- 



dard BCS superconductor described accurately by the 
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FIG. 4: An illustration of pcrturbativcly accessible and inac- 
cessible (grey) regions in the inverse particle spacing vs in- 
verse scattering length, plane around a Feshbach 
resonace, where a diverges. Note that outside the grey region, 
even for a broad Feshbach resonance there is a small param- 
eter that is either the gas parameter or Feshbach resonance 
coupling, or both, and hence the system can be analyzed per- 
turbative. 



mean-field BCS theory. If a > (the interactions are 

attractive and strong enough to support a bound state) 
and n^/'^a <C 1, the fcrmions pair up to form molecular 
bosons which then Bosc condense. The resulting molecu- 
lar Bose condensate can be studied using n^/'^a as a small 
parameter. In particular, in a very interesting regime 
where a ^ |ro| (even though a <C nT^/^) the scattering 
length of the bosons becomes approximately « 0.6a 
[35] , and the Bose condensate behaves as a weakly inter- 
acting Bosc gas with that scattering length [36] , as shown 
in Ref. [6]. Finally, when |a|n-^/'^ 1, the mean-field the- 
ory breaks down, the superfluid is said to be in the BCS- 
BEC crossover regime, and its properties so far could 
for the most part be only studied numerically, although 
with some encouraging recent analytical progress in this 
direction [31]. Much effort is especially concentrated on 
understanding the jajn""^/^ — > oo unitary regime [34] (so 
called because the fermion scattering proceeds in the uni- 
tary limit and the behavior of the superfluid becomes 
universal, independent of anything but its density). 

In this paper we concentrate solely on resonantly- 
paired superfluids with narrow resonances, amenable to 
an accurate treatment by mean-field theory regardless 
of the scattering length a. The identification of a small 
parameter [28, 31], allowing a quantitative treatment of 
the BCS-BEC crossover in resonantly-paired superfluids 
in itself constitutes a considerable theoretical progress. 
In practice most s-wave Feshbach resonances studied up 




FIG. 5; A potential with a low-energy bound state whose en- 
ergy is shown by a dashed line. If the potential is rnodifiod to 
make it more shallow, the bound state disappears altogether, 
replaced by a virtual bound state. If the potential is made 
even more shallow, a resonance - a state with positive energy 
and finite lifetime - appears. 



to now correspond to 7s — 10, which is consistent with 
the general consensus in the literature that they are wide. 
Yet one notable exception is the very narrow resonance 
discussed in Ref. [14] where we estimate 7^ ~ 0.1; for a 
more detailed discussion of this, see Sec. IX. 

Even more importantly is the observation that the per- 
turbative parameter 7 is density- (Fermi energy, (p) de- 
pendent, scaling as 7s ^ l/y^Fp, 7p ^ ^Jep for an s- 
wave and p-wave Feshbach resonant pairing respectively. 
Hence, even resonances that are classified as broad for 
currently achievable densities can in principle be made 
narrow by working at higher atomic densities. 



D. Finite angular momentum resonant pairing: 
p-wave superfluidity 

We also study a p-wave paired superfluidity driven by a 
p-wave Feshbach resonance, where the molecular (closed- 
channel) level is in the angular momentum £ = 1 state. 
While in degenerate atomic gases a p-wave superfluidity 
has not yet been experimentally demonstrated, the exis- 
tence of a p-wave Feshbach resonance at a two-body level 
has been studied in exquisite experiments in ^"^K and ^Li 
[37, 38] . Recently, these have duly attracted considerable 
theoretical attention [39-43]. 

One might worry that at low energies, because of the 
centrifugal barrier, s-wave scattering will always domi- 
nate over a flnite angular momentum pairing. However, 
this is easily avoided by working with a single fermion 
species, in which the Pauli exclusion principle prevents 
identical fermionic atoms from scattering via an s-wave 
channel, with a p-wave scattering therefore dominat- 
ing [44] . Being the lowest angular momentum channel in 
a single-species fermionic gas not forbidden by the Pauli 
exclusion principle, a p-wave interaction is furthermore 
special in that, at low energies it strongly dominates over 
the higher (than 1=1) angular momentum channels. 

There is a large number of features special to p-wave 
resonant superfluids that make them extremely inter- 
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csting, far more so than their s-wave cousins. Firstly, 
as we will show in Sec. II B and VB2, p-wave (and 
higher angular-momentum) resonances are naturally nar- 
row, since at finite density a dimensionless measure of 
their width scales as 7p = 71 ~ t^"^ (7^ ~ in the 

t angular momentum channel), that in contrast to the 
s-wave case can be made arbitrarily narrow by simply 
working at sufficiently low densities (small tp). Conse- 
quently, a narrow p-wave Feshbach-resonant superfluid, 
that can be described arbitrarily accurately [45] at suf- 
ficiently low densities for any value of detuning, is, in 
principle, experimentally realizable. 

Secondly, superfluids paired at a finite angular- 
momentum are characterized by richer order parameters 
(as exemplified by a p-wavc paired '^He, heavy-fermion 
compounds, and d-wave high-Tc superconductors) corre- 
sponding to different projections of a finite angular mo- 
mcntiim and distinct symmetries, and therefore admit 
sharp quantum (and classical) phase transitions between 
qualitatively distinct ^-wave paired superfluid gromid 
states. In fact, as we will show, even purely topologi- 
cal (non-symmetry changing) quantum phase transitions 
at a critical value of detuning arc possible [11, 41, 42, 46]. 
This contrasts qualitatively with a smooth (analytic) 
BCS-BEC crossover (barring an "accidental" first-order 
transition), guaranteed by the aforementioned absence 
of a qualitative difference between BCS and BEC paired 
superfluidity. 

Thirdly, some of the p-wave (and higher angular mo- 
mcntiim) paired states are isomorphic to the highly non- 
trivial fractional quantum Hall effect ground states (e.g., 
the Pfaffian Moore-Read state) that have been demon- 
strated to display a topological order and excitations 
(vortices) that exhibit non-Abelian statistics [11]. Since 
these features are necessary ingredients for topologi- 
cal quantum computing [47], a resonant p-wave paired 
atomic superfluid is an exciting new candidate [42] for 
this approach to fault-tolerant quantum computation. 

Finally, a strong connection to unconventional finite 
angular momentum superconductors in solid-state con- 
text, most notably the high-temperature superconduc- 
tors provides an additional motivation for our studies. 



E. Outline 

This paper, while quite didactic, presents considerable 
elaboration and details on our results reported in two re- 
cent Letters [28, 42]. The rest of it is organized as follows. 
We conclude this Introduction section with a summary 
of our main experimentally relevant results. In Sec. II we 
present general, model-independent features of a low and 
intermediate energy s-wave and p-wave scattering, with 
and without low energy resonances present. In Sec. Ill 
we discuss general features of the microscopic models of 
scattering, tying various forms of scattering amplitudes 
discussed in Sec. II to concrete scattering potentials. We 
introduce one- and two-channel models of s-wave and p- 



wave Feshbach resonances [19, 32] in Sec. IV and Sec. V, 
compute exactly the corresponding two-body scattering 
amplitudes measured in experiments, and use them to fix 
the parameters of the two corresponding model Hamil- 
tonians. These models then by construction reproduce 
exactly the experimentally-measured two-body physics. 
In the Sec. VI we use the resulting s-wave Hamilto- 
nian to study the T = narrow resonance BCS-BEC 
crossover in an s-wavc resonantly-paired superfluid, and 
compute as a function of detuning the molecular con- 
densate fraction, the atomic (single-particle) spectrum, 
the Oth-so\.md velocity, and the condensate depletion. 
In the Sec. VI B 4 contained within the Sec. VI, we ex- 
tend these results to a finite temperature. In Sec. VII 
we use the p-wave two-channel model Hamiltonian to 
analytically determine the p-wave paired ground state, 
the spectrum and other properties of the corresponding 
atomic gas interacting through an idealized isotropic p- 
wave resonance. We extend this analysis to a physically 
realistic anisotropic p-wave resonance, split into a dou- 
blet by dipolar interactions. We demonstrate that such 
a system undergoes quantum phase transitions between 
different types of p-wave superfluids, details of which de- 
pend on the magnitude of the FR dipolar splitting. We 
work out the ground-state energy and the resulting phase 
diagram as a function of detuning and dipolar splitting. 
In Sec. VIII we discuss the topological phases and phase 
transitions occurring in the p-wave condensate and re- 
view recent suggestions to use them as a tool to observe 
non-Abelian statistics of the quasiparticles and build a 
decoherence-free quantum computer. In Sec. IX we dis- 
cuss the connection between experimentally measured 
resonance width Hyj and a dimenionless parameter 7^ 
and compute the value of 75 for a couple of prominent 
experimentally realized Feshbach resonances. Finally, we 
conclude in Sec. X with a summary of our results. 

Our primarily interest is in a many-body physics of 
degenerate atomic gases, rather than in (a possibly in- 
teresting) phenomena associated with the trap. Conse- 
quently, throughout the manuscript we will focus on a 
homogeneous system, confined to a "box" , rather than an 
inhomogeneous (e.g., harmonic) trapping potential com- 
mon to realistic atomic physics experience. An extension 
of our analysis to a trap are highly desirable for a more 
direct, quantitative comparison with experiments, but is 
left for a future research. 

We recognize that this paper covers quite a lot of ma- 
terial. We spend considerable amount of time studying 
various models, not all of which are subsequently used 
to understand the actual behavior of resonantly paired 
superfluids. This analysis is important, as it allows us to 
choose and justify the correct model to properly describe 
resonantly interacting Fermi gas under the conditions of 
interest to us. Yet, these extended models development 
and the scattering theory analysis can be safely omitted 
at a flrst reading, with the main outcome of the anal- 
ysis being that the "pure" two-channel model (without 
any additional contact interactions) is sufficient for our 
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purposes. Thus, we would like to suggest that for ba- 
sic understanding of the s-wave BCS-BEC crossover one 
should read Sections II A, VA, VIA, and VI B 1. 



F. Summary of results 

Our results naturally fall into two classes of the s-wave 
and p-wave Feshbach resonant pairing for two and one 
species of fermionic atoms, respectively. For the first 
case of an s-wave resonance many results (see [9] and 
rcforenc;c;s therein) have appeared in the literature, par- 
ticularly while this lengthy manuscript was under prepa- 
ration. However, as described in the Introduction, most 
of these have relied on a mean-field approximation that is 
not justified by any small parameter and is therefore not 
quantitatively trustworthy in the strong-coupling regime 
outside of the weakly-coupled BCS regime. One of our 
conceptual contribution is the demonstration that the 
two-channel model of a narrow resonance is character- 
ized by a small dimensionless parameter 7, that controls 
the validity of a convergent expansion about an exactly- 
solvable mean- field 7 = limit. For a small 7. the por- 
turbative expansion in 7 gives results that are quantita- 
tively trustworthy throughout the BCS-BEC crossover. 
For s-wave and p-wave resonances these key dimension- 
less parameters are respectively given by: 



Is 



2 2 



1 



ni/3 (37r8)i/3' 

2 2 1/3 1^^^ 



(1.5) 
(1.6) 



where n is the atomic density, Qs and Qp are the closed- 
open channels coupling in s-wave and p-wave resonances, 
controlling the width of the resonance and m an atom's 
mass. The numerical factors in Eqs. (1.5), (1.6) are cho- 
sen for purely for later convenience. 

The many-body study of the corresponding finite den- 
sity systems is expressible in terms of physical parame- 
ters that are experimentally determined by the two-body 
scattering measurements. Hence to define the model we 
work out the exact two-body scattering amplitude for the 
s-wave [19, 48] and p-wave two-channel models, demon- 
strating that they correctly capture the low-energy reso- 
nant phenomenology of the corresponding Feshbach res- 
onances. We find that the scattering amplitude in the 
s-wave case is 



fs{k) 



1 



—a~^ + ^rok^ — ik 



(1.7) 



where Uq « (if — i?o) is the magnetic field- 

controlled detuning (in energy units), E = k^/m, and Fq, 
introduced in Eq. (1.3), is the width of the resonance, a 
and ro, which can be expressed in terms of Fq and lvq, 



represent standard notations in the scattering theory [49] 
and are the scattering length and the effective range. We 
note that tq < which reflects that the scattering repre- 
sented by Eq. (1.7) is resonant. Our analysis gives a and 
ro in terms of the channel coupling Qs and detuning coq 



ro 



Stt 



(1.8) 



In the p-wave case, the scattering amplitude is found 
to be 



fp{k) 



-v~'^ + ^kok'^ — ik^ ' 



(1.9) 



where v is the magnetic field controlled scattering vol- 
ume, and ko is a parameter with dimensions of inverse 
length which controls the width of the resonance appear- 
ing at negative scattering volume, v and ko can in turn 
be further expressed in terms of interchannel coupling gp 
and detuning ojq 



ko = - 




m'^g^ 



(1.10) 
(1.11) 



In contrast to our many-body predictions (that are only 
quantitatively accurate in a narrow resonance limit), 
above two-body results are exact in the low-energy limit, 
with corrections vanishing as 0{p/A), where A ~ l/d 
is the ultra-violet cutoff set by the inverse size d of 
the closed-channel molecular bound state. We establish 
that at the two-body level this model is identical to the 
extensively studied Fano- Anderson [25] of a continuum 
of states interacting through (scattering on) a localized 
level (see Appendix B). For completeness and to put the 
two-channel model in perspective, we also calculate the 
two-body scattering amplitude for two other models that 
are often studied in the literature, one corresponding to 
a purely local, (5-function two-body interaction and an- 
other in which both a local and resonant interactions are 
included. By computing the exact scattering amplitudes 
of these two models we show that the low-energy scat- 
tering of the former corresponds to ro ^ limit of the 
two-channel model. More importantly, we demonstrate 
that including a local interaction in addition to a reso- 
nant one, as so often done in the literature [9, 32, 33] 
is superfluous, as it can be cast into a purely resonant 
model [19] with redefined parameters, that, after all are 
experimentally determined. 

For the s-wave resonance we predict the zero- 
temperature molecular condensate density, ni, = 
\B{wo)\'^- In the BCS regime of uio ^ ^^f + 'Js^f we 
find 



nb{uJo) 



48n 



exp 



^ a;o - 2eF 

Is^F 



(1.12) 
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computed in Rcf . [6] . The crossover in the speed of sound 
as function of detuning loq should in principle be observ- 
able through Bragg spectroscopy. Extending our analysis 
to finite T, we predict the detuning-dependent transition 
temperature Tc(wo) to the s-wave resonant superfluid. In 
the BCS regime 



/ o^o^A a;o»2e^, (1.16) 



ep exp I — 



where C is the Euler constant, InC « 0.577. In the 
BEC regime Tc{uiq) quickly approachc^s the standard 
BEC transition temperature for a Bose gas of density 
n/2 and of particle mass 2m 



FIG. 6: Normalized condensate order parameter B = 
y2n6/n and normalized chemical potential jl = iijtp as a 
function of normalized detuning wo = oJo/ep- 




FIG. 7: A sketch of the critical temperature Tc{iOo) as a func- 
tion of detuning wo, displaying a maximum at intermediate 
ojQ. Tbec denotes the asymptotics of Tc at large negative uio. 



and in the BEC regime oiwo ^ —ep 



nfc(wo) 



1- 



4V2 



(1.13) 



where n is the total density of the original fermions. The 
full form of nb is plotted in Fig. 6. 

Following Ref. [50] we also compute the zeroth sound 
velocity and find that it interpolates between the deep 
BCS value of 



,,BCS 



(1.14) 



where vp 
value of 



sj2tFfm is the Fermi velocity, and the BEC 



„BEC 



5/4 r- 1 
7s£f 1 

2V4^|a;oP/2" 



(1.15) 



The BEC speed of sound quoted here should not be con- 
fused with the BEC speed of sound of the s-wave con- 
densate undergoing wide resonance crossover, which was 




Wo 4; —ep- 



(1.17) 



Taking into account bosonic fluctuations reviewed for 
a Bose gas in Ref. [51], we also observe that Tc is 
approached from above, decreased. The full 

curve is plotted in Fig. 7. In the broad-resonance limit 
of 7^ — > oo this coincides with earlier predictions of 
Refs. [4, 5, 9, 33, 48]. 

For a single-species p-wavc resonance we determine the 
nature of the p-wavc; superfluid ground state. Since the 
p-wave resonance is observed in a system of effectively 
spinless fermions (all atoms are in the same hyperfine 
state), two distinct phases of a condensate are available: 
Px + iPy phase which is characterized by the molecular 
angular momentum m = ±1 and a Px whose molecular 
angular momentum is equal to m = 0. 

We show that in the idealized case of isotropic res- 
onance, the ground state is always a. p^ + ipy super- 
fluid regardless of whether the condensate is in BCS or 
BEC regime. In the BCS limit of large positive detun- 
ing this reproduces the seminal result of Anderson and 
Morel [52] for pairing in a spin-polarized (by strong mag- 
netic field) triplet pairing in ^He, the so-called A\ phase. 
Deep in the BCS regime we predict that the ratio of 
the condensation energy Ep^+ip^ of this m = 1 state to 
the Ep^ of the competitive m = Pz state is given by 
R = Ep^+ipjEp^ = e/2, exactly. 

A much more interesting, new and experimentally rel- 
evant are our predictions for a Feshbach resonance split 
into a doublet of m = ±1 and m = resonances by 
dipolar anisotropy 6 [37]. Our predictions in this case 
strongly depend on the strength of the dipolar splitting, 
S and the resonance detuning, wq- The three regimes of 
small, intermediate and large value of splitting (to be de- 
fined more precisely below) are summarized respectively 
by phase diagrams in Figs. 9, 10, and 11. 

Consistent with above result of vanishing splitting, for 
weak dipolar splitting, < 5 < Sf^'^ we find that the 
p-wave m = 1 superfiuid ground state is stable, but 
slightly deformed to Px + ictPy, with function a{5,u)o) 
that we compute. For an intermediate dipolar splitting, 
^BEC < (5 < dc'"^ ground is a Pa, -|- ipy-superfluid 
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Px + ^Py \v 




gapless 
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FIG. 8: Temperature vs detuning phase diagrani of a p-wave 
resonant Fermi gas, for the case of no resonance sphtting, 5 — 
0, i.e., isotropic system. This phase diagram is also expected 
to describe a resonance with a splitting much larger than the 
Fermi energy for wo tuned to the m = ±1 resonance doublet. 



T 




2ep Wo 



FIG. 11: Temperature vs detuning phase diagrani of the p- 
wave resonant Fermi gas for the case of a high resonance split- 
ting, 6 > 5c'~'^ ■ This phase diagram is also expected to de- 
scribe a resonance with a splitting much larger than the Fermi 
energy for wo tuned to the m = resonance. 




(£b(H=0) 2Ep cOo 



FIG. 9: Temperature vs detuning phase diagram of a p-wave 
resonant Fermi gas, for the case of a small resonance splitting, 
Q<5< <5f 

(to = 1) in the BCS regime and is a p2:-superfluid (to, = 0) 
in the BEC regime. We therefore predict a quantum 
phase transition at Wqc between these two p-wave super- 
fluids for intermediate range of dipolar splitting [42, 53]. 
For a large Feshbach-resonance splitting, 5 > 5^^'^ the 
ground state is a stable p^-superfluid for all detuning. 
We show that in all these anisotropic cases the Px-Bxis of 
the p-wave condensate order parameter is aligned along 
the external magnetic field. Finally, wc expect that for 
an extremely large dipolar splitting, much bigger than 
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FIG. 10: Temperature vs detuning phase diagram of a p- 
wave resonant Fermi gas, for the case of an intermediate reso- 
nance splitting, S^^'^ < 5 < Sl^^^ . The critical temperature 
Tci{ciJo) vanishes in a universal way at the quantum critical 
point woc, according to Eq. (1.19). 



ep'yp (which could quite well be the current experimen- 
tal situation), the system can be independently tuned 
into TO = and to = ±1 resonances, and may therefore 
display the Px + ipy and Px states separately, depending 
on to which of the TO = OorTO = ±l resonances the 
system is tuned. Thus even in the case of an extremely 
large dipolar splitting, phase diagrams in Fig. 8 and in 
Fig. 11 will be separately observed for tuning near the 
to = 1 and to = resonances, respectively. 

As illustrated in the phase diagrams above, we have 
also extended these results to a finite temperature, us- 
ing a combination of detailed microscopic calculation of 
the free energy with more general Landau-like symmetry 
arguments. Wc show quite generally that for a dipolar- 
split (anisotropic) resonant gas, the normal to a p-wave 
superfluid transition at Tc2(a;o,^) is always into a Px- 
superfluid, that, for an intermediate dipolar splitting is 
followed by a Pa;-superfluid to Px + ipy superfluid transi- 
tion at Tci{u}o, 6). The ratio of these critical temperatures 
is set by 

^~e*/«S (1.18) 

-tcl 

where Oi, given in Eq. (7.22), is an energy scale that we 
derive. As seen from the corresponding phase diagram. 
Fig. 10, we predict that Tci{u}o) vanishes in a universal 
way according to 

Tci(a;o) ~ |wo -woc|^/^ (1-19) 

at a quantum critical point wqo that denotes a T = 
quantum phase transition between px and Px + ipy su- 
perfluids. 

In addition to these conventional quantum and clas- 
sical phase transitions, we predict that a p-wave reso- 
nant superfluid can exhibit as a function of detuning, 
luq quite unconventional (non-Landau type) phase tran- 
sitions between a weakly-paired (BCS regime of /x > 0) 
and a strongly-paired (BEC regime of < 0) versions of 
the px and px + ipy supcrfluids [10, 12, 46, 54]. In three 
dimensions these are clearly distinguished by a gapless 
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(for /i > 0) and a gapped (for /i < 0) quasiparticle spec- 
tra, and also, in the case of a + ipy superfluid via a 
topological invariant that we explicitly calculate. 

While the existence of such transitions at /i = have 
been previously noted in the literature [10, 12, 46, 54] 
our analysis demonstrates that these (previously purely 
theoretical models) can be straightforwardly realized by 
a p-wave resonant Fermi gas by varying the Feshbach 
resonance detuning, wq- 

Moreover, if the condensate is confined to two dimen- 
sions, at a positive chemical potential this state is a Pfaf- 
fian, isomorphic to the Moore-Read ground state of a 
fraction quantum Hall ground state believed to describe 
the ground state of the plateau at the filling fraction 
z/ = 5/2. This state has been shown to exhibit topo- 
logical order [11, 12], guaranteeing a 4- fold ground state 
degeneracy on the torus and vortex excitations that ex- 
hibit non-Abelian statistics. 

As was shown by Read and Green [11], despite the fact 
that both weakly- and strongly-paired p-wave superfluid 
states are gapped in the case of a p^; + iPy- (but not px- 
) superfluid the topological order classiflcation and the 
associated phase transition at jj. — remains. Consistent 
with the existence of such order, we also show [55] (via 
an explicit construction) that for /i > 0, an odd vorticity 
vortex in a Px + ipy-superfluid will gcnerically exhibit a 
single zero mode localized on it. In an even vorticity 
vortex such zero-energy solutions are absent. 

In the presence of far separated vortices, these zero- 
modes will persist (up to exponential accuracy), lead- 
ing to a degenerate many-particle ground state, and are 
responsible for the non-Abelian statistics of associated 
vortices [11, 12, 54, 55]. This new concrete realization 
of a topological ground state with non-Abelian excita- 
tions, may be important (beyond the basic physics in- 
terest) in light of a recent observation that non-Abelian 
excitations can form the building blocks of a "topolog- 
ical quantum computer", free of decoherence [47]. We 
thus propose a Feshbach resonant Fermi gas, tuned to a 
Px + ipy-supcrfluid ground state as a potential system to 
realize a topological quantum computer [42, 56]. 



II. RESONANT SCATTERING THEORY: 
PHENOMENOLOGY 

A discussion of a two-body scattering physics, that de- 
flnes our system in a dilute limit, is a prerequisite to a 
formulation of a proper model and a study of its many- 
body phenomenology. We therefore flrst focus on a two- 
particle quantum mechanics, that, for short-range inter- 
action is fully characterized by a scattering amplitude 
/(k, k'), where ±k and ±k' are scattering momenta be- 
fore and after the collision, respectively, measured in the 
center of mass frame. In the case of a centrally symmet- 
ric interaction potential U{r), the scattering amplitude 
f{k,9) only depends on the magnitude of the relative 



momentum, namely energy 

zrrir 

(with rrir = mi 7712/ (mi +m2) the reduced mass) and the 
scattering angle 9 (through k • k' = k^cosO), and there- 
fore can be expanded in Legendre polynomials, P^(cos^) 

00 

/(fc,^) = ^(2^ + l)/^(fc)P^(cos^). (2.1) 

i=0 

The scattering amplitude is related to the differential 
scattering cross section, the probability density of scat- 
tering into a solid angle Q, by a standard relation 
da/dO. = l/l^. The ^-th partial- wave scattering ampli- 
tude ft{k) measures the scattering in the angular mo- 
mentum channel I, conserved by the spherically symmet- 
ric potential U{r). For later convenience, when we focus 
on s- and p-wavc channels, we denote £ ^ and £ = 1 
quantities with subscripts s and p, respectively, as in 

/. = fe=o (2.2) 
fp = h=i. (2.3) 

In terms of the scattering matrix Sg, = e^"^^^ in channel 
deflned by a phase shift 5i, the scattering amplitude is 

given by fi = {e'^^'- - l)/{2ik). 

Analyticity and unitarity of the scattering matrix, 

[5^1 = 1, then restrict the scattering amplitude to a 

generic form 

= k-^^F,lk^)-ik ^ ^'-'^ 

where Ff {k'^) is a real function Taylor expandable in pow- 
ers of its argument [49] . It is directly related to the scat- 
tering phase shifts Se{k) through the scattering matrix 
Se = e*^'^^ via k'^^+^ cot5e{k) = Fe{k'^). Notice that at 
small k, 

h{k) - (2.5) 

Important information is contained in the poles £^poic 
of scattering amplitude (deflned by /^"^(ii^poie) = 0), 
when it is studied as a function of complex energy E. 
Poles in fi{E) correspond to discrete cigcnstates with dif- 
ferent boundary conditions that can be obtained without 
explicitly solving the corresponding Schrodinger equa- 
tion. However, because k = \/2mrE, the scattering 
amplitude, while a single-valued function of the momen- 
tum is a multi-valued function of the energy, and one 
must be careful to specify the branch on which a pole is 
located in identifying it with a particular eigenstate of 
a Schrodinger equation. Starting with a branch where 
E > Q and /c > 0, negative energy E < Q can be ap- 
proached from the positive real axis either via the upper 
or lower half complex plane. A pole which lies on the 
negative real axis, approached via the upper half plane is 
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equivalent to k = +i^2'mr\E\, i.e., Im(A:) > 0, and there- 
fore corresponds to a true bound state of the potential 
U{r), with a wavefunction ^p{r) ^ e^l'^l'" that properly 
decays at long distances. On the other hand, a pole on 
the negative real axis, approached via the lower half plane 
is not associated with a bound state, since it corresponds 
to A; = —iiy2mr\E\, i.e., Im(A;) < and therefore to 
an unphysical wavefunction that grows at largo distances 
as ijj{r) ^ el'^l''. Although it reflects a real low-energy 
feature of a scattering amplitude fe{E), the so-called vir- 
tual bound state [49] does not correspond to any physical 
bound state solution of a Schrodinger equation as it does 
not satisfy decaying boundary conditions demanded of a 
physical bound state. 

On the other hand a pole 



-2'poie = Ej. — ir/2, 



(2.6) 



of ME), with Re Epoie = Er>0, Im Ep^ie = -r/2 < 
is a resonance, that corresponds to a long-lived state 
with a positive energy E,, — Re ii^poio a-nd width T = 
—2 Im ii^poie, latter characterizing the lifetime r = l/F 
for this state to decay into a continuum. A complex con- 
jugate pole that always appears along with this resonance 
pole, corresponds to an eigenstate that is time reversal 
of the resonance solution [29, 49]. 

Coming back to the scattering amplitude Eq. (2.4), a 
low-energy scattering (small k) is characterized by a first 
few low-order Taylor expansion coefficients of F([k^), and 
therefore only weakly depends on details of the interac- 
tion potential U{r). This observation is at the heart of 
our ability to capture with a simple model Hamiltonian 
(see Sec V, below) the experimentally determined two- 
body phenomenology, governed by a complicated atomic 
interaction potential U{r) or even multi-channel model 
as in the case of a Fcshbach resonant systems. To do this 
we next specialize our discussion to a particular angular 
momentum channel. 



A. Low energy s-wave scattering 

We first concentrate on s-wave {£ = 0) scattering that, 
by virtue of Eq. (2.5), is the channel, that, for two fermion 
species dominates at low energies. 



1. Scattering in the asymptotically low energy limit 

Scattering at low energies can be analyzed by expand- 
ing the amplitude FgihF') in powers of its argument, that 
to lowest order leads to a simple form 



fs{k) 



+ ik' 



(2.7) 



with a = —1/Fg{0), where a is called the s-wave scat- 
tering length. The latter can be identified with particle 
effective interaction (in Born approximation proportional 



to a Fourier transform of the potential), with a > 
(a < 0) generally (but not always) corresponding to a 
repulsive (attractive) potential. We observe that at zero 
momentum the scattering amplitude is simply equal to 
the scattering length, /(O) = —a, leading to tr = Atto,^ 
scattering cross section. 

Wo can now give a physical interpretation to the only 
pole of Eq. (2.7) located at 



Apolc 
-^pole 



2mrd 



2 ■ 



(2.8) 
(2.9) 



The key observation at this stage is that by virtue of 
Eq. (2.8) and the fact that to be a physical bound state 



ifepolo'' 



must decay at large r, the pole Eq. (2.. 



corresponds to the true bound state with energy, E'poie 
only if a > 0. In contrast, for a < the scattering am- 
plitude pole corresponds to a wavefunction that grows 
exponential with r and therefore, despite having a neg- 
ative Spoie, is not a physical bound state or a resonance 
solution of a Schrodinger equation, but is what is called 
a virtual bound state. [49] Hence, a physical bound state 
characterized by a binding energy l/{2mra?), that van- 
ishes with 0, only exists for a positive scattering 
length a and disappears for a negative a. 

Thus, lacking any other poles at this lowest order of 
approximation, the scattering amplitude (2.7), while cap- 
turing the asymptotic low-energy bound states of the po- 
tential U{r), does not exhibit any resonances [29], i.e., 
states with a positive energy and a finite lifetime. fs{k) 
in Eq. (2.7) corresponds to a scattering from a relatively 
featureless potential of the form illustrated in Fig. 12, 
where for a sufficiently deep well, there is a bound state 
and a > 0, but only a continuum of states with a < 
and no resonance for well more shallow than a critical 
depth. This is despite the existence of an (unphysi- 
cal) virtual bound state for a < 0, with a negative en- 
ergy E = —l/{2mra?) identical to that of a true bound 
state (only present for a > 0). This point is, unfortu- 
nately often missed in the discussions of Eq. (2.7) that 
have appeared in the literature. As we discuss in de- 
tail below, this scattering phenomenology is captured 
by a featureless short-ranged attractive 2-body interac- 
tion (pseudo-potential) such as the commonly used delta- 
function four-Fermi many-body interaction. 

We notice that in order to be able to trust Eqs. (2.8) 
and (2.9), all higher order terms in the expansion of 
Fs{k'^) calculated at this value of energy have to be neg- 
ligible when compared with |fc|poic = . In other 
words, a has to be sufficiently large, and |ii^|poic suffi- 
ciently small, with precise criteria determined by the de- 
tails of the scattering potential and the corresponding 
coefficients of higher order terms in the Taylor expansion 
ofF,(A;2). 
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FIG. 12: A weakly attractive potential is progressively made 
more and more attractive, until a bound state appears, indi- 
cated in the figure by a dEished line. Although such potential 
leads to strong resonant scattering, when a bound state is 
close to or has just appeared, in contrast to a potential in 
Fig.5 it does not exhibit a resonance in a sense of a long-lived 
state with a positive energy and finite width [29] . 



amplitude is given by 
fsik) = - 



— irok'^ + + ik' 



(2.11) 



Equivalcntly, in terms of energy E = fc^/2mr (in a slight 
abuse of notation) takes the form 



fs{E) = - 



1 



r, 



1/2 



in which 



Wo 



1 Ip ro 



rrirroa 2 a 



(2.12) 



(2.13) 



and, as discussed in the Introduction, Sec. IB, a charac- 
teristic energy scale 



2. Intermediate energy resonant scattering 

In order to capture the resonant states (absent in the 
approximation Eq. (2.7)), which could be present in the 
potential U{r), Fs{k^) in fs{k) must be expanded to the 
next order in k'^, 



Fs{k^) 



1 



rok^ 



(2.10) 



with parameter ro usually called the effective range of 
the interaction potential. For a generic, everywhere at- 
tractive potential, U (r) < 0, ro can be shown to be pos- 
itive [49], and moreover, to roughly coincide with the 
spatial extent of U{r), hence the name "effective range". 
However, as is clear from physical considerations and an 
analysis of pole structure of /^(fc^), a potential which is 
attractive everywhere cannot support a resonance. In 
order to be able to capture a positive energy quantum 
particle for a significant amount of time, the potential 
must be attractive at short scales and exhibit a positive 
energy barrier at intermediate scales, of a generic form 
illustrated in Fig. 5. It can be shown that for such a 
potential, ro is in fact negative, with its magnitude hav- 
ing nothing to do with the range of U{r). Instead for 
such resonant U{r) as shown on Fig. 5, |ro|^^ reflects the 
barrier transmission coefficient, with the liiglicr barrier 
corresponding to a longer resonance lifetime and larger 
|ro|. Therefore, focusing on resonant potentials, we will 
take ro < 0, keeping in mind that |ro| can be much longer 
than the actual microscopic range of the scattering po- 
tential, d = 2n / A. In short, to leave open the possibility 
for the scattering to go in the presence of low-energy res- 
onances, in addition to bound states and virtual bound 
states, ro must be negative and "anomalously" large, a 
condition that will be assumed throughout the rest of 
this paper. 

At this higher level of approximation, the scattering 



rrirr. 



(2.14) 



r' 



is made explicit, with ro = —^y2/mrTo. It marks a 

crossover energy scale between a low- and intermediate- 
energy behaviors of fs{E). Also, as we will see below, 
Fq defines an energy scale for the low-energy pole above 
(below) which, l/m^a^ > Tq (l/m^a^ < Fq) a resonant 
state appears (disappears). 

The poles of the scattering amplitude are given by 



k^ 

'^pole 



-^pole — 



j_ _|_ V'2ai\t - 
ro aro 



(2.15) 



1 



m^rg 



a 



= too- ^Fo (l - Vl-4u;o/ro) , (2.16) 

where in -Epoic; Eq.(2.16) we only kept the "minus" pole 
(with the minus sign in front of the square-root of) fcpgig, 
as the other pole fep^j^ (with a plus sign) corresponds to 
an unphysical virtual bound state (regardless of the sign 
of the scattering length a) , and therefore will be ignored 
in all further discussions. 

The real part of the energy ii^poic, Eq- (2.16), as a func- 
tion of —a~^ (with ro < 0) is illustrated in Fig. 13. As 
—a~^ is changed from — oo to -|-oo, the pole first repre- 
sents a bound state, then a virtual boimd state (plotted 
as dotted curve), and finally a resonance. This is fur- 
ther illustrated in Fig. 14, where the position of the pole 
i^poie is shown in a complex plane of energy E, with ar- 
rows on the figure indicating its motion with increasing 
—a~^. The bound state and virtual bound state corre- 
spond to Im -Epoic = 0''', respectively, with the former 
(latter) approaching negative real axis from above (be- 
low) the branch cut. The resonance, on the other hand, 
corresponds to Im i^poie < and a positive real part of 
the energy, Re Epo\e > 0. 

We note that for l/|a| < l/|ro|, the Eq. (2.16), lying 
close to zero, approximately coincides with Eq. (2.9), as 
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bound state 



ReE 



resonance 




virtual 
bound state 



FIG. 13: The polo of the scattering amplitude fs{E), 
Eq. (2.16) as a function of — 1/a for ro < 0. As discussed 
in the text, only a bound state and a resonance correspond to 
physical solutions of the Schrodinger's equation with proper 
boundary conditions. The thin dotted line indicates asymp- 
totic linear behavior of the bound state for small positive a. 



bound state 



a = 2 r 




resonance 



FIG. 14: The pole of the scattering amplitude fs{k), 
Eq. (2.16), shown in a complex plane of E. The arrows indi- 
cate pole's motion as —1/a is increased. 



expected, since the higher order term ^rok^^i^ in 
is subdominant to ifcpoie- In other words, at a suf- 
ficiently large scattering length, the scattering is well- 
approximated by the asymptotic low-energy (scattering- 
length) approximation of the previous section. For such 
large positive a ^ |ro| it gives bound state energy, 
Eq. (2.9), that grows quadratically with 1/a. Further 
away from the resonance, on the positive a side, where 
the scattering length drops significantly below the effec- 
tive range, a <C |ro|, the bound state energy crosses over 



to a linear dependence on 1/a, as illustrated in Fig. 13 
and summarized by 



bound 



(a) 



1 

m^|ro|a ' 



for |ro| < a > 
for |ro| > 0, > 0. 



(2.17) 



More importantly, however, on the other side of the 
resonance, where the scattering length is negative, unlike 
Eq. (2.9), Eq. (2.16) also describes a resonant state, that 
appears for a < and shorter than the effective range ro, 
i.e., for — i > when the real part of the energy Epo\e 
becomes positive. The resonant state is characterized by 
a peak at energy Er = £^resonance and a width given 

by 



resonance 



rrirr. 



n 



2ro 



(2.18) 



(2.19) 



1. 



Hence, we find that in the s-wave resonant case, generi- 
cally, even potentials that exhibit a resonance for small 
\a\ (high energy), lose that resonance and therefore re- 
duce to a nonresonant case for sufficiently large and neg- 
ative a (low energy). 

The transition from a bound state to a resonance as a 
function of a is exhibited by scattering via a generic res- 
onant potential illustrated in Fig. 5. A sufficiently deep 
well will exhibit a true bound state, whose energy will 
vanish with decreasing depth and correspondingly in- 
creasing a, according to -Ebound = — l/(2TOra^), Eq.(2.9). 
We note, however, that as the potential is made even 
more shallow, 1/a crosses and the true bound state dis- 
appears (turning into an unphysical virtual bound state), 
the resonant state (positive energy and finite lifetime) 
does not appear until a later point at which scattering 
length becomes shorter than the effective range, i.e., un- 
til \a\ < |ro|. 

This somewhat counterintuitive observation can be un- 
derstood by noting that the lifetime of an s-wave reso- 
nant state is finite, given by the inverse probability of 
tunneling through a finite barrier, which only weakly de- 
pends on the energy of the state as long as the potential 
in Fig. 5 is not too long-ranged. Thus, even when the 
energy of the resonance goes to zero, its width remains 
finite. Hence, since the bound state's width is exactly 
zero, small deepening of the potential cannot immedi- 
ately change a resonance into a bound state, simply by 
reasons of continuity. There has to be some further deep- 
ening of the potential U{r) (range of a), over which the 
resonance has already disappeared, but the bound state 
has not yet appeared. During this intermediate range 
of potential depth corresponding to < l/|a| < l/|ro|, 
when the potential is not deep enough to support a true 



15 



bound state but not yet shallow enough to exhibit a res- 
onance, the scattering is dominated by a virtual bound 
state pole, as illustrated in Figs. 13, 14. 



B. P-wave scattering 

As rcmarkcid ciarlior, in a low-energy scattering of a 
particle off a potential U{r), the s-wave {£ = 0) channel 
dominates over higher angular momentum I ^ Q contri- 
butions, that by virtue of the generic form of the scat- 
tering amplitude, Eq. (2.4) vanish as fc^^, Eq. (2.5). This 
suppression for £ ^ arises due to a long-ranged centrifu- 
gal barrier, that at low energies prevents a particle from 
approaching the origin where the short range scattering 
potential U{r) resides. 

Hence, in the case of a Feshbach resonance of two hy- 
perfine species Fermi gas, where the scattered particles 
are distinguishable (by their hyperfine state), at low en- 
ergies, indeed, the interaction is dominated by the s- 
wave resonance, with higher angular momentum chan- 
nels safely ignored. However, an exception to this is the 
scattering of identical fermions, corresponding to atoms 
in the same hyperfine state in the present context. Be- 
cause Pauli exclusion principle forbids fermion scattering 
in the ,s-wavc channel, the next higher angular momen- 
tum channel, namely p-wave {£ = 1) scattering domi- 
nates, with s-wave and £ > 1 channels vanishing at low 
energies [44]. Thus wv svc that p-wavc Feshbach reso- 
nance is quite special, being the dominant interaction 
channel for a single species Fermi gas. With this motiva- 
tion for our focus on a j>wave Feshbach resonant super- 
fluidity and in preparation for its study, we next analyze 
a p-wavc scattering amplitude. 

Starting with Eq. (2.4) and expanding Fp{k'^) similarly 
to Eq. (2.10) we find 



(2.20) 



Here v is the so-called scattering volume analogous to 
the scattering length a of the s-wave case, diverging and 
changing sign when the system is taken through a p-wave 
Feshbach resonance. A characteristic wavevector fco is ev- 
erywhere negative and plays a role similar to that of the 
effective range ro in the s-wave channel, but has dimen- 
sions of an inverse length. 

Hence at low energies the p-wave scattering amplitude 
takes the form 



fp{k) 



fc2 



—v 



ik^ 



(2.21) 



Although the poles of the scattering amplitude Eq. (2.21) 
can be found by solving a qubic equation, their exact po- 
sitions are not very illuminating and will not be pursued 
here. Instead, it will be sufficient for our purpose to only 
consider an important low-energy limit \v~^\ ^ |fcoPi in 
which the relevant pole of Eq. (2.21) is close to zero and 



ReE 



bound state 




FIG. 15: The pole of a p-wave scattering amplitude Eq. (2.21) 
as a function of — 1/w for ko < 0. 



its position can be found by neglecting (actually treating 
perturbatively in powers of |ffcg|) ik^ term in the scat- 
tering amplitude. To lowest order the pole is then simply 
given by 



1 



'pole 



rrirvko 



(2.22) 



This corresponds to a real bound state for -Epoio < 
(when V > and a resonance for i^poie > (when t; < 0), 
with a width easily estimated to be 



rrirko 

-^pole? 



pole 



(2.23) 
(2.24) 



near a resonance, where \vk^\ oo. Thus, in contrast 
to the s-wave case, where at sufficiently low energies 
To = 2/nirrQ) the width F^ 



y/ToE > E, here, 

a p-wave resonance becomes arbi- 



(E <To = 2/m,r^ 

because 1 „ ^ 

trarily narrow at low energies. Consequently, as the in- 
verse scattering volume is tuned through zero and 
the relevant two-body energy iJpoie = 1/(2 rrirvko) van- 
ishes, the real bound state immediately turns into a res- 
onance without going through an intermediate virtual 
bound state (as it did in the s-wave case). This is illus- 
trated on Fig. 15. This resonant pole behavior extends 
to all finite angular momentum {£ > 0) channels. 

The physical reason behind such a drastic difference 
between s-wave and p-wave (and higher ^ > channels) 
resonances stems from the centrifugal barrier that adds 
a long-ranged l/r^ tail to the effective scattering poten- 



tial C/efF(r) = U{r) 



The width of a low lying 



resonant state in siich potential can be estimated by com- 
puting the decay rate through C/eff(?'), dominated by the 
long-ranged centrifugal barrier 1/r^. Employing the 
WKB approximation, at low energy E the decay rate is 
well approximated by 



(2.25) 
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g-2VW) dr/T 
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(2.26) 
(2.27) 



In above d and are the classical turning points of the 
— {7eff(r), where d can be taken as the microscopic range 
of the potential (closed-channel molecular size), and more 
importantly is determined by 



E. 



Combining this with Eq. (2.27) gives 



(2.28) 



(2.29) 



Although WKB approximation does not recover the cor- 
rect exponent of ^ + 1/2, Eq. (2.4) (required by unitar- 
ity and analyticity) except for the expected large £ limit 
(consistent with the fact that for small € the semiclassi- 
cal criterion on which it is based fails), it does correctly 
predict a narrowing of the resonance at low energies and 
with increasing angular momentum £. 

Of course, the expansion Eq. (2.20) is only a good 
approximation for small k. But in this regime it cap- 
tures both low energy real boimd state (for 1/v > 0) 
and narrow resonant state (for 1/v < 0) . Experimentally 
this regime is guaranteed to be accessible by tuning the 
bound state and resonance energy ii^poio = l/(mrvk()) 
sufficiently close to zero so that Ivk^l ^ 1. In this range 
the scattering amplitude Eq. (2.21) correctly captures the 
physics of a resonant scattering potential and the related 
Feshbach resonance without the need for higher order 
terms in the expansion of Fp{k'^). 



III. RESONANT SCATTERING THEORY: 
MICROSCOPICS 

A. Potential scattering 

The next step in our program is to develop a model 
of a gas of fermions interacting via a resonant pairwise 
potential U{r) of the type illustrated in Fig. 5, that ex- 
hibits a real bound state or a resonance, controlled by 
tuning its shape (e.g., well depth). It is of course pos- 
sible to simply use a many-body theory with a pairwise 
interactions literally taken to be U{r) of Fig. 5, with a 
(normal-ordered) Hamiltonian given by 

= /A.i;*i(-|l)^, (") 

where ^ai^) (^/'^^(r)) is an annihilation (creation) field 
operator of a fcrmion of flavor a at a point r. We 
would like first to discuss how a problem defined by the 
Hamiltonian, (3.1) leads directly to scattering amplitudes 
Eq. (2.4). 

Motivated by experiments where studies are confined 
to gases of no more than two fermion flavors (correspond- 
ing to a mixture of two distinct hyperfine states) we will 
refer to a as simply spin, designating a projection (a) of 
the corresponding two-flavor pseudo-spin along a quan- 
tization axis as a spin up, t, and down, J,. In an equiva- 
lently and sometimes more convenient momentum basis 
above Hamiltonian becomes 



^ = Y.Y.^''U<'^+2V^ ^ C^(|k-k'|)4,^E,^«lk'+E,<.'«-k+f,c7'ak+f,a, (3-2) 

IT k CTjct' k,k',p 



where ak,<7 (ajc „) annihilation (creation) operator of 
a fermionic atom of flavor a with momentum k, satisfying 
canonical anticommutation relations and related to the 
field operator by ipai^) = V~^/^ J2k^K<^^^^ ■ With our 
choice of momentum variables above the relative center 
of mass momenta before (after) the collision are ±k (±k') 
and p is the conserved momentum of the center of mass 
of the pair of scattering particles. 

In the rest of this section, we would like to calculate 

the scattering amplitudes ff given in Eq. (2.1) in terms 
of the interaction potential [/(|k — k'|). With this goal in 
mind, it is convenient to make the symmetry properties 
of the fermion interaction Hint explicit, by taking advan- 
tage of the rotational invariance of the two-body poten- 



tial [/(|r — r'l) and the anticommutation of the fermion 
operators. To this end we decompose the angular depen- 
dence (arising through k • k', where k is a unit vector 
parallel to k) of the Fourier transform of the two-body 
potential, C/(|k — k'|) into spherical harmonics via 

00 

C/(|k - k'l) ^ [/k,k' = J2{2£ + l)ui%,P,{k . k'). (3.3) 

1=0 

The £-th orbital angular momentum channel interaction 

amplitude j,., can be straightforwardly shown to be 
given by 

POO 

4^fe'=47r/ dr r''U{r)je{kr)M{k'r), (3.4) 
Jo 
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where ji{x) is the £-th spherical Besscl function. into the singlet and triplet channels by introducing the 

Using anticommutativity of the fermion operators, it two-body interaction vertex f/^i^^l (k, k') defined by 
is possible to decompose the interaction term in Eq. (3.2) 



1 



int - 2y 13 



with 



(3.5) 



(3.6) 



that automatically reflects the antisymmetric (under ex- 
change) property of fermions, namely 



(k,k') 



" Ualai ( k, k ) , 



(3.7) 
(3.8) 



The vertex can be furthermore decomposed into spin sin- 
glet (s) and triplet (t) channel eigenstates of the two- 
particle spin angular momentum, 



t/;^.1(k,k') = f/^f. 



.,<.^(k,k') + ^!;l„.;.,(k,k'). 

(3.9) 

The singlet and triplet vertices 



U^'^ , ,(k,k') = ?7(^)(k,k')x^'^ ' / 



(3.10) 



C.2,.;.^(k.k') = t/(°^(k,k')xl 



are expressed in terms of an orthonormal set of singlet 
and triplet projection operators 



(s) 

with coefficients 



'^ai(72<5<7'<ri), (3.11) 



2 ('^0-icri^O-2Cr2 "0-i<72"0-2 0-l 



i7(«)(k,k') = i(c/(|k-k'|) + i/(|k + k'|)), (3.12) 
[/W(k,k') = i(i7(|k-k'|)-C/(|k + k'|)), 

that, by virtue of decomposition, Eq. (3.3) and symmetry 
of Legendre polynomials, Pti—k ■ k') = (— l)^P|(k • k') 
are vertices for even and odd orbital angular momentum 
I channels, respectively, as required by the Pauli exclu- 
sion principle. Physically, these irreducible even and odd 
verticies make explicit the constructive and distructive 
interference between scattering by angle Q and tt — ^ of 
two fermions. 

The two-body scattering amplitude /(k, k') is propor- 
tional to the T-matrix, 



27r 



(3.13) 



X - 



+... 



FIG. 16: The renormalized 4-point vertex for potential scat- 
tering, determining the T-matrix T\^y . 



where = m/2 is the reduced mass of two fermions. 
The T-matrix can be computed via standard methods. 
As illustrated on Fig. 16, it equals to a renormalized 4- 
point vertex (IPI) F^^) (k+p/2, -k-^p/2, k'-hp/2, -k'-F 
p/2;£) for particles scattering with initial (final) mo- 
menta ±k -I- p/2 (±k' -I- p/2), and at a total energy in 
the center of mass frame given by 



(k + p/2)2 ^ (-k + p/2)2 



2m 



2m 



V 

4m ' 



(3.14) 



m 



m 



2mr ' 



with the last relation valid due to energy conservation by 
a time independent interaction. 

Given the retarded Green's functions of fermions, 



G(k,a;) 



w — 



2m 



(3.15) 



the main ingredient of the sequence of diagrams from 
Fig. 16 is the polarization operator, denoted by a bubble 

in the figure, and physically corresponding the Green's 
function of the reduced fermion with momentum q and 
mass m^-, 



n(q,e) 



G 



(f-q,-a;) 



1 



(3.16) 



Although perhaps not immediately obvious, n(q, e) as 
defined above is independent of p, the center of mass 
momentum of a pair of fermions. 



18 



The sequence of diagrams in Fig. 16 then generates a 
series for a T-niatrix given by 

Tk,k' = f/k,k' + «7k,qn(q, £)i7q,k' + . . . , (3.17) 



that can formally be resummed into an integral equation 

Tk.k' = [(l-C/n)-iC/]k,k', (3.18) 

where a martix product over wavevectors inside the 
square brackets is implied. 

Utilizing the channel decomposition, Eq. (3.12) of the 
vertex t/k,k' together with the closure-orthogonality re- 
lation 



- - 47r 

^F,(k.q)P,,(q-k')df^q = ^^< 



'Pe{k-k'), (3.19) 



the T-matrix series separates into a partial waves sum 



Tk,k'=E(2^+lTO'^^(k-k'), 



with 



'■k,k' 



'■k,k' 



(3.20) 



(3.21) 



a T-matrix for scattering in an angular momentum chan- 
nel i, conserved by the spherical symmetry of the two- 
body interaction potential. This demonstrates explicitly 
that the interaction vertices in different £ channels do not 
mix, each contributing only to the corresponding scatter- 
ing amplitude channel fe{k) in Eq. (2.4). 

Without specifying the interaction potential U{r), a 
more explicit expression for the T-matrix can only be 
obtained for the so-called separable potential, discussed 
in detail in Ref . [5] . Such separable interaction is a model 
that captures well a low-energy behavior of a scattering 
amplitude of a more generic short-range potential. To 
see this, we observe that a generic short-range potential, 
with a range d, leads to a vertex in the i-th channel, 
which at long scales, kd<^l, separates into 



with 



(3.22) 



(3.23) 



[(2^-M)!!]2' 

Uodl'+^ = / drr^'+^U{r). (3.24) 
Jo 

Assuming that this separation holds at all k (a definition 
of a separable potential), we use this asymptotics inside 
Eq. (3.21). This reduces the T-matrix to a geometric 
series that resums to 



'-k,k' 



(3.25) 



where 11^^^ (e) is the trace over momentum of the atom po- 
larization "bubble" corresponding to the molecular self- 
energy at energy e. 



nW(.) = l^,^5Wn(,,£),^5, 



q ' 



(27r)3£_£+,o' 



(3.26) 
(3.27) 



= -mK^'+^R [j^j - ^m'+^l^ e'+^'\ (3.28) 

In above R{x) is a Taylor-expandable function of its di- 
mensionless argument, the momentum cutoff A « 27r/d 
is set by the potential range d, and, as before, fc^/m = e. 
Putting this together inside the T-matrix, we find the 
low-energy ^-channel scattering amplitude 



1 



^(i+mA2mi?(A^))+^fc 



-, fc«:A. 



(3.29) 

This coincides with the general form, Eq. (2.4) arising 
from the requirement of analyticity and unitarity of the 
scattering matrix. However, we observe that in the s- 
wave case, for the full range of accessible wavevectors 
up to ultraviolet cutoff, fc < A the scattering ampli- 
tude Eq. (3.29) is well approximated by the non-resonant, 
scattering-length dominated form (2.7), with the scatter- 
ing length given by = An / {m\) + Atx KR{Q) . The "ef- 
fective range" ro extracted from Eq. (3.29) is ro ~ 1/A, 
namely microscopic, positive, and is of the order of the 
spatial range of the potential d. Yet, as we saw in 
Sec. II A, in order to capture possible resonances, ro must 
be negative and much longer than the actual spatial range 
of the potential. The fact that our calculation does not 
capture possible resonances is an artifact of our choice of 
a separable potential. 

Although a more physical (nonseparablc) potential 
C/(r), of a resonant form depicted in Fig. 5, indeed ex- 
hibits scattering via a resonant state (not just a bound 
and virtual bound states), calculating the scattering am- 
plitude ft,{k) (beyond Eq. (3.22) approximation) is not 
really practical within the second-quantized many-body 
approach formulated in Eq. (3.1). In fact, the only way 
to derive the scattering amplitude in that case is to go 
back to the Schrodinger equation of a pair of fermions, 
reducing the problem to an effective single-particle quan- 
tum mechanics. However, because we arc ultimately in- 
terested in condensed states of a finite density interact- 
ing atomic gas, this two-particle simplification is of little 
value to our goals. 

However, as we will show in Sees. IV and V, a signifi- 
cant progress can be made by formulating a much simpler 
pseudo-potential model, that, on one hand reproduces 
the low-energy two-atom scattering of the microscopic 
model (3.1) in a vacuum (thereby determining its pa- 
rameters by dilute gas experiments), and on the other 
hand is amenable to a standard many-body treatment 
even at finite density. 
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Furthermore, as will see below, in cases of finite angu- 
lar momentum scattering, Eq. (3.29) can in principle de- 
scribe scattering via resonances as well as in the presence 
of bound states. Thus the assumption of separability is 
no longer as restrictive as it is in the s-wave case. 

B. Feshbach- resonant scattering 

As discussed in the Introduction, in fact, the phys- 
ically most relevant resonant scattering arising in the 
context of cold atoms is microscopically due to a Fes- 
hbach resonance [17]. Generically it can be described as 
a scattering, where the two-body potential, J7a,a'(|r — r'|) 
depends on internal quantum numbers characterizing the 
two-atom state. These states, referred to as channels, are 
not eigenstates of the interacting Hamiltonian and there- 



fore two atoms coming in one channel a in the process 
of scattering will generically undergo a transition into a 
different channel a'. 

The simplest and experimentally most relevant case 
is well approximated by two channels a = o, c (of- 
ten referred to as "open" and "closed"), that approxi- 
mately correspond to electron spin-triplet and electron 
spin-singlet states of two scattering atoms; this is not to 
be confused with the hyperfine singlet and triplet states 
discussed in the previous subsection. Such system ad- 
mits an accessible Feshbach resonance when one of the 
channels (usually the electron spin-singlet) admits a two- 
body bound state. Furthermore, because pair of atoms in 
the two channels have very different magnetic moments, 
their Zecman splitting can be effectively controlled with 
an external magnetic field. The corresponding micro- 
scopic Hamiltonian is given by 



(|r-r'|) V'^(r)V'i,,,,(r')V'.',..(r')V'.,.,(r), (3.30) 



where s labels the channel. The interaction Utlt2{V ~ 
r' I ) can be more conveniently reexpressed in terms of the 
two-atom electron spin-singlct and spin-triplet channels 
basis, ?7a,a'(|r-r'|), where f/o,o(|r-r'|), J7c,c(|r-r'|) are 
the interaction for two atoms in the open (triplet) and 
closed (singlet) channels, respectively, and C/o,c(|r ~ r'|), 
characterizes the interchanncl transition amplitudes, i.e., 
the strength of o-c hybridization due to the hyperfine 
interactions. 

The corresponding scattering problem would clearly 
be even more involved than a single-channel model stud- 
ied the previous subsection. Yet, as the analysis of Sec- 
tion V will show, at low energies, the scattering ampli- 
tude of two atoms, governed by Eq. (3.30), is still of the 
same form, (2.4), as that of a far simpler pseudo-potential 
two-channel model. Indeed, the form of a scattering am- 
plitude is controlled by unitarity and analyticity, not by 
precise details of realistic Hamiltonians. Thus, to capture 
either a microscopically potential- or a Feshbach resonant 



scattering wc will replace a realistic Hamiltonian. such 
as Eq. (3.30) with a simpler model, which, nevertheless 
exhibits a low-energy scattering amplitude of the same 
form. To this end, in the next two sections we examine 
two such effective models and work out their scattering 
amplitudes. We will thereby determine and justify our 
subsequent choice of a many-body model with the correct 
low-energy two-body physics. 

IV. ONE-CHANNEL MODEL 

A. 5-wave scattering 

The most drastic simplification of a resonant Fermi 
gas is to model the two-body interaction by a featureless 
and short-ranged single- channel pseudo-potential, that at 
long scales and low energies is most commonly taken to 
simply be U{r) = \S^^\r), with the corresponding many- 
body Hamiltonian 



H^-'"^ = y dV[^^t (^_|!^ + A ^|(r)^|(r)^T WVii(r)" • 

(4.1) 

I 

In analyzing the Hamiltonian like this one, one has to function potential is known to have a vanishing scattering 
exercise a certain amount of caution, as the repulsive S- amplitude in three dimensions, and therefore does not 
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make sense if understood literally [57] . 

Hence (5-function potential must be supplemented with 
a short-scale cutoff 1/A (i.e., given a finite spatial extent), 
that we will take to be much smaller than the wavelength 
of a scattering particle, i.e., k/K <^ 1. Furthermore, 
for calculational convenience, but without modifying the 
properties on scales longer than the cutoff, we will im- 
pose the cutoff A on each of the momenta k and k' in- 
dependently, modeling the interaction in Eq. (4.1) by a 
featureless separable potential 



C/k,i 



A ^(A^ - k^) ^(A^ - k'2), (4.2) 



a 
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with Q{x) the usual step function, and interactions in 

all finite angular momentum channels vanishing by con- 
truction. We note that this separability of the potential 
is consistent with the general long wavelength form of 
a generic short-scale potential foimd in Eq. (3.22), al- 
though it does lead to some minor unphysical features 
such as only a single bound state, independent of how 
strongly attractive the potential (how negative A) is [5]. 

As discussed in Sec. Ill A, the Dyson equation (3.21) 
can be easily resummed into Eq. (3.25), with the s-wave 
polarization bubble H^e) = W=^'^{e) (cf Eq. (3.26)) eas- 
ily computed to give 



n«(£) = 



d^q 0{A - q) 

(2^)3£_2i+,0 

m 



(4.3) 



where we used s <C h? /{2m). This then directly leads to 
the s-wave scattering amplitude (vanishing in all other 
angular momentum channels) 



/.(k,k') 



iZL _|_ 2A _|_ ,-1.' 



(4.4) 



which coincides with Eq. (2.7), where the scattering 
length is given by 



, 47r 2A\ m ^ 



m 



(4.5) 
(4.6) 



47rl-A/Ac' 

where A^ can be called the renormalized coupling and 

27r2 



Ac = - 



Am 



(4.7) 



is a critical value of coupling A at which the scattering 
length diverges. 

Hence we find that scattering off of a featureless po- 
tential of a microscopic range 1/A (modeled by the cut- 
off 5-function separable potential), is indeed given by 
Eq. (2.7), with this form exact for k/A <C 1, i.e., for the 



FIG. 17: The scattering length a as a function of strength A 
of the separable (5-function potential. 



particle wavelength 1/k longer than the range of the po- 
tential. We also note that in the limit A — > oo, the scat- 
tering amplitude vanishes, in agreement with the afore- 
mentioned fact that the ideal 5-function potential does 
not scatter quantum particles [57]. 

For finite cutoff, the scattering length a as a function 
of A is shown on Fig. 17. We note that in the "hard 
ball" limit of a strongly repulsive potential, A ^ 
the scattering length is given simply by its spatial ex- 
tent, a = 7r/(2A) ~ d. For an attractive potential the 
behavior is more interesting. For weak attraction, the 
scattering length a is negative. However, for sufficiently 
strong attractive potential, i.e., sufficiently negative A, 
the scattering length a(A) changes sign, diverging hyper- 
bolically at the critical value of Ac, and becoming positive 
for A < Ac. The critical value of A at which this takes 
place corresponds to the threshold when the potential 
becomes sufficiently attractive to admit a bound state. 
There is no more than one bound state in a separable 5- 
function potential, regardless of how strongly attractive 
it is [5]. 

Finally, as above discussion (particularly, Eq. (4.4)) 
indicates, although a one-channel s-wave model can suc- 
cessfully reproduce the very low energy limit of the 
generic s-wave scattering amplitude, such ultra-short 
range pseudo-potential models cannot capture scattering 
via a resonance. The actual Feshbach resonance experi- 
ments may or may not involve energies high enough (large 
enough atom density) for the scattering to proceed via a 
resonance (most do not, with the criteria for this derived 
in Subsection C, below). However, our above findings 
show that the ones that do probe the regime of scatter- 
ing via a resonance must be described by a model that 
goes beyond the one-channel ^-function pseudo-potential 
model. We will explore the simplest such two-channel 
model in Sec. VAl. 



21 



B. Finite angular momentum scattering 

Unlike their s-wave counterpart, one-channel models 
for higher angular momentum scattering can describe 
scattering via resonances. This is already clear from the 
analysis after Eq. (3.29). Let us analyze this in more 
detail. 

Above s-wave model (4.1) can be straightforwardly 
generalized to a pseudo-potential model at a finite angu- 
lar momentum. This is most easily formulated directly 
in momentum space by replacing the two-body ^-wave 
interaction in the microscopic model (3.2) by a separable 
model potential 



(4.8) 

k^) 6I(A2 - k'2). (4.9) 



that simply extends the long wavelength asymptotics of a 
microscopic interaction Eq. (3.22) down to a microscopic 
length scale 27r/A. 

Using results of the previous section, this model then 
immediately leads to the scattering amplitude Eq. (2.4) 
with 



47r 



2^2^+1 



2A2^- 



+ (4.10) 



m\ 7r(2^+l) 7r(2£-l) 
The corresponding scattering amplitude is given by 



fm = 



1-) rvn rv Crv 



(4.11) 
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with the analogs of the scattering volume (of dimensions 
2^+1) and effective range parameters given respectively 
by 



k 



ve 



2e-i 



^ 'k{2£ + 1) 

4A2^-i 
"7r(2£- 1)' 
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Af, (4.12) 



(4.13) 



We note that diverges (hyperbollically) for a suffi- 
ciently attractive interaction coupling, reaching a critical 
value 



27r2 {2i + 1) 



From the structure of fi{k) it is clear that at low ener- 
gies (length scales longer than k^^), the imaginary term 
if^'^W jg subdominant to the second k"^ term in the de- 
nominator. Consequently, the pole is well-approximated 
by 

£^poie « (4.14) 



where we defined 



'''pole 



-2k, 



-2^+l„ -1 



(4.15) 
(4.16) 



For a positive detuning, < 0, leading to the first term 
of -Bpoie real and positive, while the second one —iT/2 
negative, imaginary and at low energies {k < ko) much 
smaller than Re[Ep]. Thus, (in contrast to the s-wave 
case) for finite angular momentum scattering, even a 
single-channel model with a separable potential exhibits 
a resonance that is narrow for large, negative vg. For 
Vi > 0, the term —iT/2 becomes real and this resonance 
directly turns into a true bound state, characterized by 
a pole i^poie) tliat is real and negative for i;^ > 0. 



C. Model at finite density: small parameter 

Having established a model for two-particle scattering 
in a vacuum, a generalization to a model at finite density 
n, that is of interest to us, is straightforward. As usual 
this is easiest done by working within a grand-canonical 
ensemble by introducing a chemical potential ji that cou- 
ples to a total number of particles operator N via 



(4.17) 



One thereby controls the average atom number and den- 
sity by adjusting /it. 

The single-channel models of the type Eq. (4.1) and 
its corresponding finite angular momentum channel ex- 
tensions have been widely studied in many problems of 
condensed matter physics. Although (as most interacting 
many-body models) it cannot be solved exactly, for suf- 
ficiently small renormalized couphng Xr{X), (4. 5), (4. 13), 
(whether positive or negative), we expect that one can 
analyze the system in a controlled perturbative expan- 
sion about a mean-field solution in a dimensional mea- 
sure of Xr, namely in the ratio of the interaction energy 
to a typical kinetic energy e^. 



1. Small parameter in an s-wave model 

In the s-wave case this dimensionless ratio is just the 
z,as parameter 



(4.18) 



For weak repulsive s-wave interaction, A > 0, an^^^ 1, 
and the perturbation theory generically leads to a Fermi 
liquid [581. For weak attractive interaction A < and 



1/3 



<C 1, it predicts a weak-coupling BCS supercon- 



\a\n 
duct or. 

However, as A is made more negative (increasing the 
strength of the attractive interaction) \a\ increases ac- 
cording to Eq. (4.5), as illustrated in Fig. 17 and eventu- 
ally goes to infinity when A reaches the critical value of 
Ac. Near this (so-called) unitary point, the gas parame- 
ter is clearly large, precluding a perturbative expansion 
within a one-channel model. 
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On the other (BEG) side of the unitary point, a molec- 
ular bound state appears and the phenomenology is that 
of interacting bosonic molecules with a molecular scat- 
tering length proportional to that of fermionic atoms, 
am ~ 0.6a [6, 35, 59]. Since on the BEG side a also di- 
verges (this time from a positive side), the bosonic gas 
of these molecular dimers is strongly interacting near the 
unitary point and the situation is as hopeless for quanti- 
tative analysis as it was on the BCS side of the unitary 
point. 

Yet, at large negative A, the bound state drops to a 
large negative energy and a becomes small again (this 
time positive). In this deep-BEG regime, the resulting 
dilute repulsive gas of tightly bound molecules then also 
exhibits the same small gas parameter as that deep in the 
BGS regime. Hence its ground state is a weakly inter- 
acting superfluid Bose-condensate [4, 5], with properties 
that can be computed perturbatively in a small parame- 
ter av}/^, although careful analysis of this sort was only 
done recently [6]. 

We note in passing, that, at a finite atom density the 
effective measure of the strength of interaction is actu- 
ally a dispersive coupling Xkp, given by the T-matrix 
Tk,,k, = (47r/m)|/,(fcF)| 



A 



kp 



kp\fs{kp)\, 

kp 

\a~^ + ikp\ ' 

kp\a\, for kplal <C 1, 
1, for kplal 1, ■ 



(4.19) 
(4.20) 
(4.21) 

(4.22) 



Thus, in contrast to a two-body case, at finite density the 

growth of this effective dimcnsionlcss coupling, ^ fci?|a|, 
actually saturates at 1 (i.e., at a large, nonperturbative, 
but noninfinite value), due to a cutoff of the growing 
scattering length a by atom separation kp^. 

Hence, despite its many successes to predict qualitative 
behavior, the Hamiltonian Eq. (4.1) has a limited abil- 
ity to describe a resonant interacting Fermi gas. First 
of all, as we just saw, its two-body scattering amplitude, 
as given by Eq. (2.7), does not describe scattering via a 
resonant state, capturing only a true bound state (for a 
sufficiently attractive A and positive a), but not a reso- 
nance [29] (possible for negative a). Thus, if resonances 
(states at positive energy and finite lifetime) are present, 
the model given by Eq. (4.1) is insufficient. Even in the 
absence of such resonant states, the perturbation the- 
ory about the mean-field state commonly used to ana- 
lyze Eq. (4.1) breaks down in the course of the BGS to 
BEG crossover, where the scattering length surpasses the 
inter-particle spacing and |o|n^/"^ > 1 is no longer small. 

However, it is quite common in literature to ignore 
these issues and simply extend the mean-field analy- 
sis of Eq. (4.1) into the nonperturbative unitary regime 
near Ac. Given the absence of a phase transition in 
the s-wave case, the prediction of such mean-field the- 



ory is undoubtedly qualitatively correct even in the 
strong coupling regime that smoothly interpolates be- 
tween Pauli-principlc stabilized large Gooper pairs and 
a BEG of tightly bound molecules. However, as we 
just discussed, such approach (all the perturbative em- 
bellishments notwithstanding) cannot make any quan- 
titatively trustworthy predictions for A w Ac, a regime 
where a bound state is about to, or just appeared and 
jajn"'^/^ ':§> 1. Since the question of the s-wavc BGS- 
BEG crossover is intrinsically a quantitative one, quan- 
titatively uncontrolled studies performed within above 
nonresonant model provide little information about the 
details of such crossover, particularly near the so-called 
unitary regime. 



ar momentum model 



2. Small parameter in a finite 



As can be seen from the form of the scattering am- 
plitude (4.11) and its parameters and ko, the case of 
a gas resonant at a finite angular momentum is qualita- 
tively quite different from that of the ,s-wave model just 
considered. The reason is that, as discussed in Sec. IV B, 
on length scale longer than the spatial range of the po- 
tential A"-'^ (i.e., on effectively all accessible scales) the 
j^2^4-i -j^ /^(fc) is subdominant and a one-channel finite 
angular momentum model exibits a resonant state that 
continuously transforms into a bound state. As discussed 
in Sec. II B, physically this stems from the existence of a 
finite £ centrifugal barrier that strongly suppresses the 
molecular decay rate at low positive energies. 

Analogously (but distinctly) to the s-wave case, a 
dimensionless parameter that measures the relative 
strength of interaction and kinetic energy in the ^-wave 
case is given by 



Afe^(^) 



l-^kp I 



kf n 



ep 



\\v7'\+kr'kj.\' 

' kf+^\vel for kl'^'^kHvel < 1, 
for kl^-^k%\vi\ » 1, 



(4.23) 
(4.24) 



(4.25) 



Since, (as found above) ko ~ A, we find that for a fi- 
nite angular momentum resonance, although the effective 
coupling A^^^ grows with V£ as a resonance is approached, 
it saturates at a value <C 1, cutoff by a finite density. 
Thus, this heuristic argument suggests that in princi- 
ple a controlled perturbative treatment of full BGS-BEG 
crossover is possible for finite angular-momentum Fesh- 
bach resonances, even within a one-channel model. Such 
analysis has not yet been done, and it is an interesting 
research problem for future work. 
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V. TWO-CHANNEL MODEL 

As we have seen, there are considerable shortcomings 
of a local onc-channcl model, particularly for the s-wavc 
case, as it does not exhibit a resonant state, nor does 
it have a dimensionless parameter that can be taken to 
be small throughout the BCS-BEC crossover. Thus, wc 
now consider a more involved fermion-boson two-channel 
model that is free of these deficiencies. Furthermore, the 



appeal of this two-channel model is that it is inspired by 
and more acciiratcly reflects the microscopies of the Fesh- 
bach resonance physics discussed in the Introduction and 
above, but applies more universally to any system where 
a resonant interaction (e.g., a shape resonance of the type 
illustrated in Fig. 5) is at work. A general two-channel 
model Hamiltonian, that in cold-atom context for the 
special s-wave case was first introduced by Timmermans 
[19], is given by 




In above model Hamiltonian nk.tj (oJc a) is a fcrmionic an- 
nihilation (creation) operator of an atom of flavor a with 
momentum k, representing atoms in the open-channel 
(typically corresponding to the electron [physical, as op- 
posed to flavor] spin-triplet state of two atoms) con- 
tinuum. The annihilation operator h p.e (k) destroys a 

<t,<t' 

bosonic diatomic molecule of mass 2m, with a center of 
mass momentum p, internal (atoms') momenta ±k. It 
is a cartesian spin-tensor that transforms as a tensor- 
product of two spin-1/2 representations and an orbital 
angular momentum £ representation. It is convenient to 
decompose it into 2^-|- 1 components bp,t,m correspond- 

ing to its projections along an orbital quantization axis, 
according to: 



symmetric matrices, with the relations 



&p..^(k) = ^ 6p,.,™r,:„(k) 

a, a a, a 

■m=-e 
= bp,. -YKk), 



(5.2) 



where l£_m(k) are the spherical harmonics, k is a unit 
vector along k, and in the last line the scalar product is 
over 2£ + 1 components labeled by — £ < m, < L These 
bosonic orbital components can be further decomposed 
into a singlet = s) and a triplet (with three spin pro- 
jection components = 0,±1 linear combinations of 
/i = {x, y, z) cartesian components) spinor representa- 
tions according to: 



(5.3) 



with a a vector of Pauli spin matrices; notice that iay 
is a fully antisymmetric (and thus a singlet) spin tensor 

and the components of ayff are linear combinations of the 
spin-triplet projections Sz = 0, ±1, represented by 2 x 2 



:(o) 



.(±1) _ 



(5.4) 
(5.5) 



Within a Feshbach resonant system context the molecule 
bp ^ represents a (quasi-) bound state of two atoms in 
a closed channel (usually electronic spin-singlet state of 

two atoms) , a true bound state in the limit of a vanishing 

coupling g'^^ (proportional to o-c channels hybridization 
energy Uoc) for the decay of a closed- channel molecule 
into an open-channel pair of atoms. As discussed in the 
Introduction, in this case the 'bare' molecular rest en- 
ergy eQ^'™-* (the detuning relative to the bottom of the 
open-channel contimmm) corresponds to the Zeeman en- 
ergies that can be readily tuned with an external mag- 
netic field. For generality we allowed this detuning to 
have a nontrivial m dependence, encoding an explicit 
breaking of orbital rotational invariance seen in the ex- 
perimental systems [37]. This ingredient will be central 
to our analysis in Sec. VII, for a determination of the cor- 
rect ground state of a p-wave paired superfluid. Focusing 
on the closed-channel bound state, the model clearly ig- 
nores the continuum (with respect to relative coordinate) 
of closed-channel states. Because for the experimentally 
interesting regime of a resonance tuned to low energies, 
these states are at a finite energy, they can be adiabat- 
ically eliminated (thereby only slightly modifying model 
parameters) and can therefore be safely omitted. In the 
context of a shape resonance the molecule bp ^ represents 
a resonance that is long-lived in the limit of large poten- 
tial barrier. 

We would like to emphasize that for a nonzero Fes- 
hbach resonance coupling g'^^ it would be incorrect to 

consider 6-particles to be the true bound states (diatomic 
molecules) of a-particles (atoms). Indeed, freely propa- 
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gating 6-particlcs arc not even cigcnstatcs of the Hamil- 
tonian Eq. (5.1). The true, physical molecule is a linear 
combination of 6-particles and a surrounding cloud of a- 
particles. They can be found by studying the scatter- 
ing problem posed by Eq. (5.1). In particular, the true 
bound states of Eq. (5.1) can be spatially quite large 
with their spatial extent set by a scattering length a and 
at finite atom density can easily overlap. In contrast, 
the ^-particles (related to the true bound states only in 

the limit of vanishing Feshbach resonance couplings 5^^^) 
are point-like, with their size set by a microscopic length 
scale corresponding to the range d = 27r/A of the inter- 
atomic atomic potential, Utlt^{\r — r'|). 

With this in mind, we now turn to the analysis of 
the two-channel model, considering separately the s-wave 



and p-wave cases. 



A. S-wave 

1. Two-atom scattering 

As discussed in the Introduction, in the case of a two- 
flavor atomic gas, at low energies it is appropriate to 
focus on the dominant s-wave channel, which by virtue of 
Pauli principle automatically also selects the singlet two- 
atom states. Ignoring all other scattering channels in the 
model (5.1), the s-wave two-channel model Hamiltonian 
reduces to 
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where to simplify notation we defined 



"p — "p,0,0> 
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(5.7) 
(5.8) 



incorporating the short-scale (shorter than the atomic 
interaction range = 27r/A) falloff of the Feshbach reso- 
nant coupling g^^'' as an implicit sharp cutoff at A on the 
momentum sums. 

Within this model, the fermions of the same spin do 
not interact at all, and the scattering amplitude of two 
fermions of opposite spin can be calculated exactly. In 
addition to the free fermion Green's function Eq. (3.15), 
the Green's function for a free boson is given by 



-Do(p,w) 



1 



(5.9) 



The T-matrix is then given by a geometric series depicted 
in Fig. 18, and written algebraically as 



r(k, k') = gsDogs+gsDogsU.s9sDogs-\-. . . 
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Above, 



(5.10) 



= 7^o( P,£+|^ ) =I?o(0,/c7m), (5.11) 



(5.12) 



latter defined in Eq. (4.3) and e = k^ /m, is the two-atom 
center of mass energy. The scattering amplitude is then 



given exactly by the low-to-intermediate energy form in 
Eq. (1.7), with the s-wave scattering length and effective 
range given by 



a = — 



m9s 



ro = 
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(5.13) 



(5.14) 



We note that as expected from general considerations 
of Sec. II A 2, the effective range parameter ro is indeed 
negative. Proportional to l/g^, it controls the lifetime of 
the 6-bosons to decay into two atoms, corresponding to 
the inverse width of the Feshbach resonance, ro therefore 
becomes arbitrarily long (compared to the microscopic 
range 1/A) of the potential with decreasing gg. We recall 
from Sec.IIA2 that these two conditions are precisely 
those required for the resonance to exhibit a positive en- 
ergy, finite lifetime resonant state. 

From Eqs.(5.13,5.14) we also identify the characteristic 
crossover energy scale Fq = 4/(mrQ) and the parameter 
Wo = 2/(mroa) appearing in the s-wave scattering am- 
plitude (2.12), given by 



To 
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In terms of these derived quantities, all the scattering 
phenomenology discussed in Sec. II A 2 follows immedi- 
ately. 

We observe that the scattering length diverges at a crit- 



ical value of the bare detuning eg 



-A, correspond- 
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FIG. 18: The diagrams contributing to the T-rnatrix of the 
two-channel model Eq. (5.6). The straight and wavy hnes rep- 
resent fermionic and bosonic Green's functions, respectively. 

ing to the point at which the bound state appears. This 
should be contrasted with the naive expectation that the 
bound state, being a 6-particle, appears when eo goes 
through zero. We refer to this shift as a renormahza- 
tion of detuning ("mass rcnormaUzation" of the closed- 
channel boson b in the field theory parlance). The origin 
of the shift from eo to uq lies in the fact that the 6-particle 
is, of course, not the bound state (physical molecule) of 
the two a-particles (atoms). Rather, an actual bound 
state is a superposition of a ^-particle and a cloud of a- 
particles. The 6-particlc only corresponds to the part of 
the physical bound state (molecule) which lies within the 
closed-channel. Wc emphasize that while the 6-particle 
can be safely treated as a point particle, whose size is 
related to the detuning- independent cutoff 1/A, the size 
of the actual bound state (physical molecule) can get ar- 
bitrarily large, with its size diverging with a ^ oo, as is 
further discussed in Appendix A. 

Since a diverges where the parameter cjq vanishes, we 
identify this additively-renormalized detuning parameter 
Wo with the physical detuning corresponding to Zeeman 
energy splitting between closed- and open channels, con- 
trolled by the magnetic field H and vanishing at field Hq. 
Comparing the prediction (5.13) for a with its empirical 
form, Eq. (1.1), [24], allows us to identiiy parameters of 
the two-channel model with the experimental parameters 
according to 

uo « 2hb{H-Ho), (5.17) 
To « AmixlalgHl/h\ (5.18) 

where we estimated the magnetic moment responsible 
for the Zeeman splitting between the open- and closed- 
channels (corresponding respectively to electron spin 
triplet and singlet, respectively) to be 2/1^. 

Hence, the conclusion is that indeed the two-channel 
model faithfully describes a scattering in the presence of 
a resonant state, as well as a bound and virtual bound 
states, depending on the value of the detuning parame- 
ter. It is thus a sufficient model to capture all the generic 
features of a Feshbach-resonant atomic gas, without re- 
sorting to a fully microscopic (and therefore typically in- 
tractable) description. This should be contrasted with 
the one-channel model Eq. (4.1), which is only able to 
capture scattering in the absence of a resonant state, i.e., 
only in presence of either bound or virtual-bound states 
and as such insufficient to capture an intermediate energy 
behavior of a Feshbach-resonant atomic gas. 

We close this section with a comment. In the litera- 
ture it is common to study models that in addition to 



the two-channel Feshbach resonant interaction consid- 
ered above, a featureless nonresonant four-Fermi atomic 
interaction is also included. It is simple to show that in 
three-dimensions, doing so does not add any new physics 
to the pure two-channel model considered here. Instead 
it just amomits to redefining the relation between model's 
parameters (eo, 5s, A) and the experimentally determined 
parameters ojq and ro- Please see Appendix B2 for de- 
tails. 



2. Model at finite density: small parameter 

As already discussed in the Introduction, a finite den- 
sity ,s-wavc resonant Fermi gas, described by a two- 
channel model, Eq. (5.6) is characterized by an average 
atom spacing n~^/^ ~ l/fci? in addition to the scattering 
length a and the effective range ro, derived above and 
discussed phenomenologically in Sec. II. Hence, in addi- 
tion to the dimensionless gas parameter kpa discussed in 
Sec. IV C 1, a two-channel model admits another key di- 
mensionless parameter 7^ oc l/(fci?|ro|) that is the ratio 
of the average atom spacing fc^^ to the effective range 
length |ro|. Equivalently, 7^ is related to the square-root 
of the ratio of the Feshbach resonance width To (con- 
trolled by the Feshbach resonance coupling and defined 
by Eq. (2.12)) to the Fermi energy, and equivalently to 
the ratio of the resonance width (at the Fermi energy) 
s/Tqcp to the Fermi energy: 

»^^v5 /r^^^^s 1 

ep TT V ep ve? 7rfcF|ro| 

The two-channel model, Eq. (5.6) is described by an in- 
teracting Hamiltonian, whose interaction strength is con- 
trolled by a coupling gg. The corresponding dimension- 
less parameter 7^ oc g1 controls a perturbative expansion 
in the Feshbach resonant interaction (about an exactly 
solvable non-interacting = limit) of any physical 
quantity. The key observation is that 7^ is independent 
of the scattering length a and detuning wq, and as such, 
if indeed small, remains small throughout the crossover, 
even for a Feshbach resonance tuned through zero en- 
ergy, where the scattering length a and the gas parameter 
|a|n^/^ diverge. 

Hence, we arrive an important conclusion: the two- 
channel model predictions for a narrow Feshbach reso- 
nance, (defined by 7^ <C 1, i.e., width of the resonance 
To much smaller than the Fermi energy, or equivalently 
effective range ro much longer than atom spacing n~^/^) 
are quantitatively accurate throughout the BCS-BEC 
crossover, no matter how large the value of the gas pa- 
rameter |a|n^/^ gets. 

As discussed in Sec. IV C 1, this availability of small pa- 
rameter in the two-channel model contrasts strongly with 
the one-channel model, characterized by a dimensionless 
gas parameter lajn^/"^ that diverges for a Feshbach reso- 
nance tuned to zero (i.e., "on resonance"), that therefore 
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docs not admit a small perturbative expansion parame- 
ter throughout the entire crossover (with the exception 
of deep BCS and deep BEC regimes). 

In contrast, for the broad-resonance 7s ;» 1 system, 
the two-channel model is no more solvable than the one- 
channel model; in fact, as we will show in the next sub- 
section, in this limit the two models become identical. 
The perturbativcly accessible and nonpcrturbative re- 
gions of the two-channel model in the kp and a"^ param- 
eter space are illustrated in Fig. 4. In terms of the Fig. 13, 
the broad and narrow resonance limits respectively corre- 
spond to fci? falling inside and outside the virtual bound 
state regime, defined by l/\a\ < l/(2|ro|). 

The dimensionless parameter 7s naturally emerges in a 
perturbative expansion in atom-molecule coupling. More 
physically, it can also be deduced by estimating the 
ratio of the energy associated with the atom-molecule 
Feshbach-resonance interaction to the kinetic energy, i.e, 
the non-interacting part of the Hamiltonian Eq. (5.6). 
To see this, note that the atom-molecule coupling energy 
^FR P^^ atom scales like 



kF\fs{kF) 



(5.25) 
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where we estimated the value of 6(r) by 5(r) ~ y/n. 
This interaction energy is to be compared to the non- 
interacting part of the Hamiltonian, i.e., the kinetic en- 
ergy per atom 



Eq ~ ep, 
with the square of the ratio 

7s ~ {E].n/Eo) 
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giving the scale of the dimensionless parameter 7s in 
Eq. (5.19). 

In the spirit of the discussion in Sec. IV CI, another 
instructive way to estimate the interaction strength and 
to derive the dimensionless coupling that controls the 
two-channel model's perturbation theory is to integrate 
out (in a c;olicrcnt-statc path-integral sense) the closed- 
channel molecular field 6(r) from the action. As 6(r) 
couples to atoms only linearly this can be done exactly 
by a simple Gaussian integration. The resulting action 
only involves fermionic atoms that interact via an effec- 
tive four-Fermi dispersive vertex. After incorporating 
fcrmion-bubble self-energy corrections of the T-matrix 
the latter is given by Tkp « {Att /m)fs{kp), with a key 
factor that is the finite-density analog of the scattering 
amplitude, fs{k), Eq. (2.11). To gauge the strength of the 
molecule-mediated interaction energy we compare the in- 
teraction per atom {Att /m)fs{kp)n to the kinetic energy 
per atom ep. Hence, dropping numerical prefactors, the 
dimensionless coupling that is a measure of the atomic 
interaction, is given by 
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(5.24) 



At large detuning (i.e., deep in the BCS regime) A|p 

kp\a\ <^ 1 and the theory is perturbative in A|^ . How- 
ever, as detuning is reduced |a(wo)| and A^ (wq) grow, 
and close to the resonance may be neglected in the 
denominator of Eq. (2.11). In this regime, the coupling 



saturates at A^ : 



■\ oo 
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whose magnitude crucially depends on the dimensionless 
ratio 7s oc l/(fcF|''o|), with 



\ oo 



for \ro\kp » 1, 
for |ro|A;F <C 1. 



rok 
1 



(5.27) 



Hence, in contrast to a two-particle vacuum scattering, in 
which the cross section diverges when the Feshbach reso- 
nance is tuned to zero energy, at finite density, for suffi- 
ciently large a, the effective coupling A|^ ceases to grow 

and saturates at A^ , with the saturation value depend- 
ing on whether this growth is cut off by the atom spacing 
1/kp or the effective range tq. The former case corre- 
sponds to a narrow resonance [7s oc {\rQ\kp)~^ <C l],with 
the interaction remaining weak (and therefore perturba- 
tive) throughout the BCS-BEC crossover, right through 
the strong-scattering l/(fcp'|a|) = point. In contrast, in 
the latter wide-resonance case [7s oc (|ro|fcF)~^ 1], dis- 
cussed in Sec. IV C 1, sufficiently close to the unitary point 
1/a = the effective coupling A^, Eq. (5.27), grows to 
0{1) precluding a perturbative expansion in atom inter- 
action near the unitary point. 



3. Relation to one-channel model 

In this section we would like to demonstrate that 
in the broad-resonance limit, of relevance to most 
experimentally-realized Feshbach resonances to date, the 
rok"^ contribution to the dispersion (arising from the 
molecular kinetic energy) of the effective coupling A| can 
be neglected and one obtains an effective single (open-) 
channel description. Thus the one and two channel mod- 
els are equivalent in the limit 7s ^ 00. 

The reduction to a single-channel model in the broad 
resonance limit can be executed in an operator formalism, 
with the derivation becoming exact in the infinite Fesh- 
bach resonance width (7s 00) limit. (For this same 
reduction in the functional integral formalism, see Ap- 
pendix A of Ref. [6].) The expression for the scattering 
length, Eq. (5.13) 
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dictates that a proper transition to the broad resonance 
limit corresponds to gs ^ oo while adjusting the bare 
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detuning according to 



(5.29) 



such that the physical scattering length a remains fixed. 
This allows us to trade the bare detuning ep and coupling 
Qs for a new coupling A that physically corresponds to 
a non-resonant attractive interaction depth, that can be 
used to tune the scattering length. The Heisenberg equa- 
tion of motion governing the molecular field bp dynamics 
under Hamiltonian (5.6), with condition Eq. (5.29), is 
given by: 















LV 4m 





(5.30) 
(.5.31) 



Now, in the large — > oo limit (keeping A fixed) the 
molecular kinetic energy term cx /Am on the right and 

the 6p term on the left are clearly subdominant, reducing 
the Heisenberg equation to a simple constraint relation 



-k+Ex"k+E|- 



(5.32) 



Hence, we see that in the extreme broad-resonance limit 
the molecular field's dynamics is "slaved" to that of the 
pair of atoms, according to Eq. (5.32). Substituting 
this constraint into the Hamiltonian, (5.6) allows us to 
eliminate the closed-channel molecular field in favor of a 
purely open-channel atomic model with the Hamiltonian 



Hi 
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(5.33) 



k,k',p 



a momentum space version of the one-channel model, 
Eq. (4.1) discussed in Sec. (IV). 

A clear advantage of the one-channel model is that, as 
shown above, it naturally emerges as the correct Hamil- 
tonian in the experimentally-relevant case of a wide res- 
onance, 7s » 1. However, as discussed in Sec. IV CI, 
a notable disadvantage is that, in the most interesting 
regime of a Feshbach resonance tuned to zero energy, its 
dimensionless gas parameter kF\a\ oo precluding a 
controlled perturbative calculation throughout the BCS- 
BEC crossover. 



B. P-wave 

1. Two-atom scattering 

As discussed in the Introduction, for a single compo- 
nent Fermi gas Pauli principle forbids interaction in the 



s-wave channel, and, consequently the dominant interac- 
tion is in the p-wave channel. In addition to this motiva- 
tion, a study of a p-wave resonance is attractive because, 
as we will see below, (and is already clear from scatter- 
ing phenomenology discussion in Sec. Ill A) they can in 
principle be made arbitrarily narrow by simply decreas- 
ing the particle density (as opposed to increasing n in 
the s-wave case) and therefore arc amanablc to a quanti- 
tatively accurate description possibly in experimentally 
accessible regimes. Finally, as we will see, p-wave super- 
fluids exhibit richc^r set of possibilities and thereby allow 
genuine phase transitions (some quite exotic), not just 
crossover as a function of detuning. 

With this motivation in mind, in this section we focus 
on the dominant p-wave, Sz = +1 triplet channel (a gas 
of atoms in a single hyperfinc state |) in the model (5.1), 
described by the following p-wave Hamiltonian 
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(5.34) 



Here, as before, we defined the p-wave coupling to be gp, confused with momentum. To simplify notation, we de- 
where the subscript p refers to the "p-wave" , not to be fined three (cartesian tensor components) p-wave bosonic 



28 



operators {a = x,y,z) in terms of the three bosonic 
(closed-channel) operators &p*i^^' with definite projec- 
tions of orbital angular momentum, m 
Sz = +1), defined in Sec.V 
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,(+1) , 
,(+1) 

'p,l: 
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(±1,0) (and 

(5.35) 

(5.36) 
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and we have dropped the hyperfine subscript t on these 
molecular operators. We also defined the corresponding 
Feshbach resonance coupling and bare detunings 
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incorporating the short-scale (shorter than the atomic 

interaction range d = 2tt/ A) falloff of the Feshbach reso- 
nant coupling g^^ as an implicit sharp cutoff at A on the 
momentum sums. The coupling gp is the amplitude for 



the transition between a pair of identical fermionic atoms 
with one unit of orbital (relative) angular momentum into 
a closed-channel molecule with an internal angular mo- 
mentum £ ~ 1. 

In Hp, Eq. (5.34), we have specialized to the ex- 
perimentally relevant time-reversal invariant Hamilto- 
nian [37] for degenerate m = ±1 resonances and thereby 
omitted a contribution 

Ht-break = ^(4^'"^^^ " 4^'"^^) Xl^^P.y^P,^ " ^p,x^P,!/)' 

(5.41) 

that vanishes in the case Co^'^^' = Cq^' ^'^ of interest to 
us here. 

By construction, the fermionic atoms (a-particles) 
scatter only in the p-wave channel. The scattering am- 
plitude can be easily calculated in the same T-matrix 
formalism, as in the s-wave case, Eq. (5.10). 

The propagator of the fea-particles is given by 



Da/sip, = 



D„{p,uj)Sc,0, (5.42) 



Graphically, the T-matrix is represented by the geometric 
series in Fig. 18, with vertices proportional to fc^. It is 
given by 



Tk,k' = 2gl ^ kaDak'^ + 2g^V-^ ^ kaD^qa 2U qnDpk'p + . . . = ^ 
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where stands for ^0, , 11 stands for the po- 
larization bubble n ^q, defined in Eq. (3.16), a = 
X, y, z, and overall factor of 2 comes from the definition of 
the T-matrix in this many-body language (see factor of 
1 /2 in the definition of the interaction term in Eq. (3.2)). 
A related symmetry factor of 2 appearing in front of 11 
in Eq. (5.43) is also a consequence of identical fermions 
appearing the diagrams in Fig. 18, that allows two possi- 
ble contractions of atomic lines inside 11, which contrasts 
to one such contraction for s-wave scattering of atoms 
distinguished by (hyperfine-) spin. 

Calculating the momentum q sum in the p-wave po- 
larization bubble in the denominator of Eq. (5.43), we 
find 
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where as before e = k'^/m is the molecule's internal en- 
ergy in the center of mass frame. Just as in Eq. (4.3) in 
lip{e) we have cut off the (otherwise ultra- violently di- 
vergent) integral at high momentum A corresponding to 
the inverse (closed-channel) molecular size, with the cal- 
culation (and the whole approach of treating 6 as a point 
particle) valid only as long as e <^ A^/m. However, in 
contrast to the ,s-wavc, Eq. (4.3), here the integral for the 
p-wave case scales as at large momenta. As a resiilt, 
in addition to the constant contribution (first A'' term 
that is analogous to hnear A term, Eqs.(4.3), (5.16)) that 
leads to the detuning shift, the polarization bubble shows 
a second A-depcndent contribution that multiplic;atively 
renormalizes the molecular dispersion. For a future refer- 
ence, we introduce two cutoff-dependent parameters re- 
lated to these two terms 
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where ci is a constant with dimensions of energy, C2 is an 
important dimensionless constant and we defined a new 
momentum scale 



kg — 



(5.47) 



The two-body scattering ampUtude is obtained 
through its relation /(k,k') = -^Tk,k', Eq. (3.13) to 
the T-matrix. Combining this with Eq. (5.43), Eq. (5.44) 
and Eqs. (5.45), (5.46), we thereby obtain 



/(k,k') = E 



(5.48) 

For an isotropic interaction, m = 0, ±1, Feshbach reso- 
nances are degenerate, = eo) and the scattering ampli- 
tude is (not surprisingly) entirely in the p-wave channel 



/(k,k') = 3/p(fc)cos(^), 



(5.49) 



where 9 is the angle between momenta k and k' before 
and after the scattering event. The partial wave scat- 
tering amplitude fp{k) in the p-wave channel, as follows 
from Eq. (2.1) and (3.20), is given by 



fp{k) 
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Therefore, as argued on general grounds, the scattering 
of identical fermionic atoms is indeed exactly of the form 
Eq. (2.21), with 



ko = - 



2 (^0 - Cl) 

mgp 

Yin 



2^2 



m _ 

kn 



(I + C2), 



1 + 



(5.51) 
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and this result essentially exact, valid on all momentum 
scales up to the cutoff A. As required for a resonant state, 
indeed fco < is negative definite. 

We note that fco is the characteristic momentum scale 
beyond which the width of the resonance becomes 
larger than its energy k^k"^ /2, i.e., a crossover scale be- 
yond which the resonant state disappears. It is clear from 
its form, Eq. (5.52) that fco is given by the following limits 



fco = X 
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for kg > A, 

for kg < A, ^^-^^l 



depending on the ratio A/fcg, but with A:o > A for all kg, 
set by Feshbach resonance coupling gp. 

It is useful to introduce the physical detuning wq 
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that corresponds to the energy of the pole in fp(k), 
Eq. (5.50) when this pole is tuned to low energy. In 
terms of the detuning u)q, the p-wave scattering ampli- 
tude is given by 



fp{k) = 



{l + C2){uJo-^)-ik^' 



(5.55) 



Adjusting luq from negative to positive, turns the scat- 
tering in the presence of a low- lying bound state at — |a;o| 
into the scattering in the presence of a resonance at loq. 

Thus, the p-wavc two-channel model Eq. (5.34) cap- 
tures the most general low-energy scattering in almost 
exactly the same way as its s-wave counterpart does. 
The most obvious difference from the ,s-wavc model lies 
in how the cutoff A enters the scattering amplitude. In 
the s-wave case A could be eliminated via a redefinition 
of the detiming energy from eo to ivq, as in Eq. (5.16), and 
thereby disappears from all other computations. In the 
p-wave case, however, A enters the scattering amplitude 
not only additively but also multiplicatively, and there- 
fore explicitly appears in the scattering amplitude fp{k), 
Eq. (5.55) through the dimensionless parameter C2, even 
after the shift to physical detuning ojq- 

Interestingly, appearance of C2 persists when we calcu- 
late the phases of the p-wave condensate later in this pa- 
per. While the parameter ci drops out of all predictions 
when written in terms of physical parameters, the dimen- 
sionless parameter C2, controlled by the closed-channel 
cutoff A and coupling gp continues to appear explicitly. 
Unfortunately, it is not easy to extract C2 from experi- 
mental measurements of the scattering amplitude, as it 
enters the amplitude in the combination {1 + C2) / {rn^ g^) . 
We note that if C2 ^ 1 then it and the uv-cutoff A in- 
deed drop out from all physical quantities with fco ~ ~kg. 
However, if C2 » 1, then fco ~ — (1 -|- C2)/(m^g'p) —A 
reduces to a quantity that is completely independent of 
g. Indeed in this limit, the bare dispersion of the closed- 
channel 6-field can be ignored (in comparison to the po- 
larization bubble A-dependent corrections) and the field 
b can be integrated out, just like in the strong coupling 
,s-wavc two-channel model, and leads to a p-wave single- 
channel model analog of (4.1). 



2. Model at finite density: small parameter 

A finite density p-wave resonant Fermi gas is charac- 
terized by the following three length scales: the average 
atom spacing (Fermi wavelength) l/kp, the analog of 
effective range (characterizing resonance intrinsic width 

1/3 

Fo) 1/fco, and the scattering length \v\ . Consequently, 
we can form two dimensionless constants. One is the 
p-wave gas parameter fci?|w|^^'^, that, although small for 
large positive and negative detuning, diverges near the 

resonance (tuned to low energy), and thereby precludes 

1/3 

a controlled perturbativc expansion in fci? I u I (or equiv- 
alently in nv) throughout the phase diagram. However, 
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in the two-channel p-wave model the second dimension- 
less parameter, defined in Eq. (1.6) and approximately 
given by 



7p 



(5.56) 



( that is approximately independent of Wq and v), offers 
such a controlled expansion even in the region where v 
diverges. 

As in the s-wave case, we can get a better physical 
sense of the dimensionless parameter that controls per- 
turbation theory by looking at the ratio of the p-wave 
Feshbach resonance interaction energy 



gpTi^^'^kF, 



to the typical kinetic energy per atom 
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where we estimated the value of 6(r) by 6(r) ~ i/n. and 
took the value of a typical internal momentum k (appear- 
ing in the j>wave vertex) to be fcir. We find 
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We first note that, in principle, at high densities (ener- 
gies) kp > ko, the kokp is subdominant and the dimen- 
sionless coupling is given by the gas parameter vkp vn 
that saturates at 1 as v^^^ grows beyond atom spacing 
1/kp. However, given the ko asymptotics, Eq. (5.53), we 
have kp ko > A and in contrast to the s-wave case this 
nonperturbative regime is never accessible in the p-wave 
case. Namely, at all physically accessible densities the 
width term kp is subdominant and the effective coupling 
is given by 
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kpM = ^{kokl\v\), forfcofc|k|«l, 
for kokllvl » 1. 
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Since, as emphasized above kp <^ ko for all densities 
of physical interest (see Eq. (5.53)), we conclude that a 
p-wave resonant gas is always in a pcrturbative regime 
with the perturbation theory controlled via a small di- 
mensionless coupling 7p given in Eq. (5.56). It thereby 
allows quantitatively accurate predictions given in pow- 
ers of 7p (However, if C2 ^ 1, this description might go 
beyond mean field theory [45]. We are grateful to Y. 
Castin [60] for pointing this out to us.) 



that, as expected is indeed controlled by ■jp in the limit 
of large kg, i.e., small gp. 

As in the s-wave case above we can more carefully 
gauge the interaction strength and the corresponding di- 
mensionless coupling of the p-wave two-channel model 
by formally integrating out the closed-channel molecular 
field b(r). This leads to an effective one-channel model, 
with atoms that interact via an effective four-Fermi dis- 
persive vertex. After incorporating the fcrmion-bubble 
self-energy corrections of the T-matrix the latter is given 
by Tkp « {4:Tr/m)fp{kp), with a key factor that is the 
finite-density analog of the scattering amplitude, fp{k), 
Eq. (2.21). Following the finite angular momentum one- 
channel model analysis of Sec. IV C 2 we can gauge the 
strength of the molecule-mediated interaction energy by 
comparing the interaction per atom {Air /m)fp{kp)n to 
the kinetic energy per atom ep- Hence, dropping numer- 
ical prefactors, the dimensionless p-wave coupling that is 
a measure of the atomic interaction, is given by 



47rn 



\fp{kp)\, 



kplfpikp) 



kp 



\koki/2-ik%y 



(5.60) 



(5.61) 



VI. 5- WAVE BCS-BEC CROSSOVER 

In Sec. V A we developed and justified the proper two- 
channel model of an s-wavc resonantly interacting atomic 
gas and related its parameters to a two-body scattering 
experiment. We now turn to the main goal of our work, 
namely a study of this model at a fixed chemical poten- 
tial, with the aim to establish the thermodynamics of an 
s-wave resonant Fermi gas as a function of temperature 
T, density n, and detuning loo- 

The thermodynamics is encoded in the partition func- 
tion Z = Tr e~^^ and the corresponding free energy F = 
—T\nZ. The partition fimction can be conveniently for- 
mulated in terms the imaginary-time path-integral over 
coherent states labeled by commuting closed-channel 
fields (/>(r), (^(r) (bosonic molecules) and anticommut- 
ing open-channel fields tp{Y), ^p{r) (fermionic atoms), and 
their complex conjugates 
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with the action [4>, "0] corresponding to the Hamilto- 
nian H^'"^, Eq. (5.6), given by 
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Eq. (6.1), they can be formally integrated out, giving 
an effective purely bosonic action 



where (3 = 1/T is the inverse temperature. In above, 
because fermionic atoms and bosonic molcciilcs are in 
equilibrium, able intcrconvert into each other via the Fes- 
hbach resonant coupling Qg, the chemical potential cou- 
ples to the total conserved particle number, that is the 
sum of the number of free atoms (fermions) and twice 
the total number of bosons. 

Given that fermions ■0 appear quadratically in where the effective action Sg [(f)] is given by 
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The bosonic action S's[0] completely characterizes our 
system. However, it is nonlinear in 0, describing effective 
bosonic interactions controlled by gs and therefore can- 
not be solved exactly. Nevertheless it can be studied via 
standard many-body methods as we will describe below. 



A. Infinitely-narrow resonance limit 

It is enlightening to first consider the limit of a vanish- 
ingly narrow resonance, 7^ ~ f/gm'^/^ep,^^^ — » 0. As can 
be seen most clearly from the original action, Eq. (6.2) or 
the corresponding Hamiltonian, Eq. (5.6), in this — > 
limit the system breaks into two non-interacting parts; 
fermionic atoms of mass m and bosonic molecules of mass 
2m. Despite of a vanishing interaction, we emphasize an 
implicit order of limits here. Namely, the vanishing in- 
teraction is still sufficiently finite so that on experimental 
times scales the resulting fermion-boson mixture is never- 
theless in equilibrium, with only total number of particles 
(but not the separate fermion and boson number) 



renormalized physical detuning Eq. (5.16), disappears 
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that is conserved, with n/, Uh the atom and molecule 
densities, respectively. This key feature is captured by a 
common chemical potential /U even in the gig ^ limit. 
In this limit the boson lifetime, even for eo > when 
the boson is actually a resonance [29], becomes infinitely 
long (but still short enough to establish equilibrium with 
fermions), and therefore the bosons can be considered as 
stable particles on equal footing with fermions. Also in 
this limit the distinction between parameter eo and the 
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At T = 0, the condition 2/i < loq holds, since the lowest 
energy level of bosons cannot be negative. We therefore 
arrive at the following picture illustrated in Fig. 19. For 
a large positive detuning loq > 2ei?, molecules are too 
energetically costly to be produced in equilibrium, and 
all particles are fermionic atoms forming a Fermi sea, 
with the chemical potential locked to the Fermi energy 
ji = ep, with 
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set completely by the total particle density n. However, 
as illustrated in Fig. 19, for an intermediate range of de- 
tuning wq < 2eF, it becomes energetically advantageous 
to convert a fraction of fermions in the Fermi sea between 
ujQ and 2eF into Bose-condensed molecules, thereby keep- 
ing the effective bosonic chemical potential 2/i — at its 
lowest value of zero. This atom-to-molecule conversion 
regime continues as detuning is reduced, with /x locked 
to ci;o/2. It terminates when wq reaches 0, at which point 
atom-to-molecule conversion is complete and the system 
enters into the BEG regime of a pure molecular conden- 
sate for Wo < 0. The full range of behavior can be sum- 
marized by the evolution of the molecular boson density, 
nb(wo) with detuning luq, that can be easily found to be 



rib 



0, 



t(l-(e)'''' 



for luq > 2ep 
for < luq < 2eF 
for Wo < 0. 
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FIG. 19: An illustration of the BCS-BEC crossover in the 
limit of a vanishing resonance width 7s —» 0. The evolution 
with detuning tJo is illustrated, with (a) the BCS regime of 
ujo > 2eF, where particles are all free atoms forming a Fermi 
sea, (b) the crossover regime of < ojo < 2eF, where a frac- 
tion of atoms between u>o and ep have converted into EEC of 
bosonic molecules, and (c) the EEC regime of luq < 0, where 
only Eose-condensed molecules are present. 



and is displayed in Fig. 20. 

For finite temperature the chemical potential is no 
longer locked to the detuning and is determined by the 
particle number equation Eq. (6.5), together with the 
noninteracting expressions for the fermionic density 
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FIG. 20: The normalized density of bosonic molecules hb = 
2nb /n vs the normalized detuning lio = ioo/ep in the limit of 
a vanishing resonance width, 73 ^ . 
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(6.11) 

that uniquely gives Tc{ujq). The numerical solution of 
Eq. (6.11) is presented in Fig. 21. 

The limiting behavior of Tciujo) is easy to deduce. 
Deep in the BEC regime, for lo^ <C — e_F, the first in- 
tegral is exponentially small, reflecting the fact that in 
this regime the fermion chemical potential /i is large and 
negative and a number of thermally created fermionic 
atoms is strongly suppressed. The second integral then 
gives the critical temperature, that in this regime coin- 
cides with the BEC transition temperature 



Tc(wo < -e_F) ~ Tbec = — 



TO V2C(3/2) 



2/3 



(6.12) 



and the bosonic density 



that is indeed on the order oi ep- As o^o is increased 
through the BEC and crossover regimes, Tc{uj[j) de- 
creases, as the contribution of thermally created free 
atoms from the first integral increases. When detuning 
reaches wq = 2eF, the solution of Eq. (6.11) drops down 
to Tcilep) — 0. Beyond this point, for loq > 2ep in the 
BCS regime, the bosons are completely converted into 
free fermions forming a Fermi sea and T(,{loo) sticks at 0. 



/d k 1 
(^ ,^~2^_i +"°' ('-^'^ 

where no = \B\'^ is, as usual, the density of the bosonic 
condensate. This total number constraint must be sup- 
plemented by the free-energy minimization rule that 
no > only if 2/i = loq and vanishes otherwise. 

These equations can then be used to determine the 
normal-superfluid transition temperature Tciuj^), defined 
as a temperature at a given detuning wo at which no first 
vanishes. Setting 2/z = and no = 0, we find an implicit 



B. Narrow-resonance limit 

We extend our study of the s-wave resonant Fermi gas, 
described by the Bose-Fermi mixture action to the limit 
where 7s is small but nonzero. The overall qualitative 
picture is quite similar to the —t Q limit discussed in 
the previous subsection and summarized in Fig. 19, with 
only a few new features. 

Because of the nonlinearities in [0], Eq. (6.4), the 
functional integral Eq. (6.3) in general cannot be evalu- 
ated exactly. However, as discussed in the Introduction 



33 



T 



c 



-3 



-1 




FIG. 21: The normalized critical temperature Tc = Tc/Tbec 
as a function of the normalized detuning uiq = wo/ef in the 
limit of a vanishing resonance width 7^ — > 0. 



and in Sec. VA, for small gs {'js ^ 1) the theory can 
be analyzed by a controlled perturbative expansion in 
powers of 7s around the saddle-point (mean-field) ap- 
proximation of Zg. To this end, we look for the spatially 
uniform field configuration (j){r) — B that minimizes the 
action Sg [(j)] ■ We find the following saddle-point equation 



1 SSg [0] 
B 66 



2m 



= 0. 



g'sBB 



(6.13) 

where w„ = 'KT{2n + 1) are the fermion Matsubara fre- 
quencies. The sum over the frequencies can be done in a 



closed form, leading to the so-called BCS-BEC gap equa- 
tion for the mean field B{T, jj, eg) (and the corresponding 
condensate density |-B|^) 



where is given by 



2 r (fk tanhf^ 



(6.14) 



(6.15) 



The integral on the right hand side of the BCS-BEC gap 
equation is formally divergent, scaling linearly with the 
uv momentum cutoff A. However, expressing the bare 
detuning parameter in terms of the physical, renor- 
malized detuning using Eq. (5.16), we can completely 
eliminate the appearance of the microscopic uv scale 
A in all physical quantities, and thereby obtain a uv- 
convergent form of the BCS-BEC gap equation 



wq — 2fi = 



(Pk 
(2^ 



tanh^ 
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(6.16) 



To calculate Zg in the saddle-point approximation, we 
write the bosonic field 



cl>{r)=B + ip{r), 



(6.17) 



in terms of a fluctuation field i^(r) about the saddle point 
B. Expanding the action, we find 



Ss[B,^]^So[B]+S(iuct[B,ip], 



(6.18) 



where 



So[B] = Ss 



^(eo-2M)BB-y^ J 



dPk 



In 



^-H +9IBB, 



(6.19) 



and 



S20 (P, ^n) 
^Sb{p, -^n) + Sii (p, -U.n) 



'f-p,-n„ 



(6.20) 



Here Sb{p, i^n) = i^n + ^0 — + ^ represents the free 
part of the bosonic action, Sn {p, lu) and S02 {p, i-^) arc the 
normal and anomalous fermion polarization operators, 
and 0„ = 27rTn are the bosonic Matsubara frequencies. 
We used here a fact that S20 is symmetric under the 
sign change of its arguments. The explicit expressions 
for the polarization operators are not very illuminating 
at this stage and will be discussed later (see Eqs. (6.71) 
and (6.72)). Within this saddle-point approximation, the 



partition function Zg is given by 




VifVifexpi-Sfiuct)- (6.21) 



In the applications of Eq. (6.14) to atomic gases, it is 
the total particle number N, rather than the chemical 
potential that is controlled experimentally. Of course, as 
usual, in the thermodynamic limit there is no distinction 
between the two ensembles and it is sufficient to work in 
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the grand-canonical ensemble (as we have done above), 
and then eliminate /x in favor of N through the particle 
number equation 



N = T^lnZ,. 

O/J, 



(6.22) 



Solving (6.22) simultaneously with the BCS-BEC gap 

equation (6.16) determines the condensate density B and 
chemical potential /i as a function of experimentally con- 
trolled parameters, the detuning Uq, temperature T, and 
particle number N. 



1. Zero temperature BCS-BEC crossover 

At zero temperature the BCS-BEC gap equation re- 
duces to 



wo - 2/1 = ^ 



f (Pk 


■ 1 2m" 
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(6.23) 



The particle number equation can also be evaluated not- 
ing that at T = and small 7s, most of the weakly- 
interacting bosons remain in the condensate B (with only 
a small interaction-driven depletion set by 7^) and there- 
fore the fluctuations (p{y) are small and can be safely 
neglected. Omitting S'fluct from Zg above and using 
Eq. (6.19), the particle number-density (n = N/V) equa- 
tion is then given by 



-I 



TdSo 
V djji 
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where Ek is still given by Eq. (6.15), and 



(6.24) 



(6.25) 



Before solving this equation together with the BCS- 
BEC gap equation, (6.23) for the condensate B and the 
chemical potential /z as a function of uiq, let us comment 
on the nature of the ground state. The bosons, within 
the narrow resonance approximation, where the conden- 
sate depletion has been neglected, are all located in the 
condensate. The fermions, on the other hand, form a 
paired superfluid described by the BCS-like Hamiltonian 

H - n Nf = ^^k a^^akCT + 

k,<T 

9s [b al^a^_j^^ + B a-kiakt) , (6.26) 

k 

with the condensate B appearing as pairing parameter 

to be self-consistently determined. In practice, B can 
always be chosen to be real due to the symmetry B 



e}fB, so that B = B. The role of the BCS gap is played 
by the bosonic density, via 



(6.27) 



The ground state of this Hamiltonian is the BCS wave- 
function [1] 

|BCS) = n ("fc + 4T«-fci) |0) . (6.28) 



with Ufe, Vk given by the standard expressions 




-^). (6.29) 



When Wo is large and positive, wq > 2eF, the closed- 
channel molecules (bosons hp) are energetically sup- 
pressed and Eq. (6.26) predicts phenomenology of a BCS 
superconductor [1]. Namely, most of the particles will 
be fermionic atoms with a weak attraction due to ex- 
change of virtual bosons (resonances) with energy de- 
tuned much higher than the chemical potential. Such 
degenerate fermions will therefore form a BCS ground 
state (6.28). Due to fermionic pairing, the bosons will 
still be present, albeit in the exponentially small num- 
bers since the BCS gap Eq. (6.27) and the associated 
Cooper-pair density are exponentially small. 

In a qualitative picture similar to the (7s ^ limit 
(see Fig. 19), as ujq is decreased the number of bo- 
son condensate will grow, while the number of fermions 
will decrease, with the substantial change approximately 
(within a window 7s cf) taking place when cjo drops be- 
low 2eF- As Wo crosses 0, the chemical potential that 
tracks it for wq < 2eF will also change sign and become 
negative. At this point the remaining fermions will form 
a "strongly-coupled superconductor" (in the notation of 
Ref. [11]) that exhibits pairing in the absence of a Fermi 
surface, driven by the existence of a true two-body bound 
state. Such situation is not typically encountered nor 
experimentally accessible in ordinary, solid-state super- 
conductors. As discussed in Ref. [11], a strongly coupled 
s-wave superconductor is not qualitatively different from 
the more standard one with /i > 0. This contrasts with 
the p-wave case, where /z > and jj < regimes are sep- 
arated by a (possibly topological) quantum phase transi- 
tion, that we will discussed in Sec. VII devoted to p-wave 
resonances. Because throughout the entire range of de- 
tuning wq (particularly for < wq < 26^) the system 
will be a superposition of Bose-condensed bosons and 
Cooper-paired fermions, we refer to it as a BCS-BEC 
condensate. 

It is useful to contrast this small 7^ picture of BCS- 
BEC condensate with the earlier studies of wide reso- 
nances (corresponding to the large 7s limit in our set- 
ting). When the resonance is wide, most of the particles 
are fermions regardless of the value of wq (open-channel 
fermions in the atomic physics parlance) and only a small 
fraction of the condensate will be bosons (closed-channel 
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fermions). As a result, for a wide resonance no sharp 
features exist in the BCS-BEC crossover. 

Now we are in the position to solve the BCS-BEC gap 
and the particle number equations. The former can be 
written as 



Wo — 2/x = 



47r2 



I{u), 



(6.30) 



where u = QsB / ^ (with B = B), and the integral I{u) is 
given by 



i{u) = r 

Jo 



dx 



\J (x^ — sign uf' +11? 



- 1 



(6.31) 



with sign u = sign /i. Similarly, the particle number 
equation reduces to 



where K{u) is defined as 



K{u) + = N, 



(6.32) 



K{u) 



dx 







1 



sign u 



\/ [x"^ — sign u)2 + 

(6.33) 

K{u) essentially measures the deviation of the particle 
number from the usual Fermi distribution in the absence 
of pairing, with K{u — > 0+) — > 1. 

In the BCS regime we expect a small condensate with 
w <C 1, for which 



I{u) w In 



(6.34) 



The logarithmic divergence of /(w) as u goes to zero is the 
standard Fermi surface contribution to the gap equation. 
In the same regime we can replace K{u) by K{0'^) = 1. 
We thus find 



B ^8e ^/uexp 



Att'^ (wq - 2/i) 
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and 
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+ 2B^ = n. 



(6.35) 
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Since B is exponentially small, Eq. (6.36) gives n ~ ep- 
Combined, these two equations give us B{luq) as a func- 
tion of detuning. The quantitative validity of this regime, 
u <C 1 is given by 



(6.37) 



where 7^ is the small parameter given by Eq. (5.19) char- 
acterizing an s-wave resonant gas at density n. 

As wq is decreased below 2ep^ a crossover regime is 
entered as fermions begin to be converted into bosons. 



and the chemical potential /i tracks wo/2 to accuracy 
0(75). As is clear from the infinitely narrow resonance 
analysis, Eq. (6.8) and Fig. 19, most of the atoms will 
pair up into bosons that become true bound states inside 
the EEC regime of wo < 0. In this BEC regime, for a 
sufficiently negative wo, so that u < and \u\ <C 1, we 
can use the following simple asymptotics /(— 0) = — f , 
and K{u) w ff w^, to reduce the gap and particle number 
equations deep in the BEC regime to be 



Wo — 2/i 
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These give a BCS-BEC condensate, in which most parti- 
cles are molecular (true bound states) bosons, and only a 
small fraction of the total number are the Cooper-paired 
fermions. 

The full solution of Eqs. (6.30) and (6.32) can only be 
found numerically and is displayed in Fig. 6 for 7 = 0.1. 
We can see that in contrast of the infinitely narrow res- 
onance case presented earlier on Fig. 20, the molecular 
condc;nsate density rib = B^ extends to eo > 2eF, al- 
though only as an exponentially small tail. This repre- 
sents the BCS condensate absent in the limit of an in- 
finitely narrow resonance, 7^ — > 0. 



2. Ground state energy across BCS-BEC crossover 

It is instructive to also calculate the zero tem- 
perature grand canonical ground-state energy density 
£gs{I-^, Wo, -B), that is given by the T ^ limit oiTSo/V, 
where Sq is given by Eq. (6.19). Calculating it at arbi- 
trary B and then minimizing it with respect to B will 
be of course be equivalent to solving the gap equation 
Eq. (6.23). As in the previous subsection, to lowest or- 
der in 7s, ignoring quantum fluctuations of (I), we find 

B^ 

£GS = (wo - 2ii)B^ + V(a -Ek+ gl — ). (6.40) 
k 

Here we again chose B to be real, traded eo for wo, and 
used the notations Eqs. (6.25), (6.15) with the additional 
notation 



2m' 



(6.41) 



Setting B = Q the normal state energy is easily com- 
puted as 



eGs{B = Q) = ^(^,-1^,1), 



(6.42) 



= -^cn^'Ml^), (6.43) 
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where we converted the sum to an integral and used the 
three-dimensional density of states N{E) = c\fE with 



3/2 



(6.44) 



Combining this with Eq. (6.40) then gives: 



£GS = (wo - 2/z)S^ 
where 

J(/i,B) = 



ciJ.^/^e{n) + J{n,B), (6.45) 



15 



{2n) 



(6.46) 



and we have converted the momentum sum to an integral. 

The gap equation discussed previously, Eq. (6.23), and 
the particle number equation Eq. (6.24) obviously follow 
from 







and 
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For a narrow Feshbach resonance (7^ <SC 1), we can find 
an accurate analytic approximations to eas in Eq. (6.45) 
in all relevant regimes. The first step is to find an ap- 
propriate approximation to Eq. (6.46), which has drasti- 
cally different properties depending on whether /i > (so 
that the low-energy states are near the Fermi surface) or 

< (so that there is no Fermi surface and excitations 
are gapped with energy bounded from below by |/i|). We 
proceed by first evaluating the derivative ^ and then 
integrating the expression with a constant of integration 
J(Ai,0) = 0. 



dB = -^^^ 



(fk,l 1 . 
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, M > 0;)i:^ gsB, (6.51) 
y ,fx«-gsB, (6.52) 



This calculation proceeds through the evaluation of the 
integral for I{u), Eq. (6.31). Integrating back up with 
respect to B, we thus have 
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a. BCS regime: The BCS regime is defined by ujq 3> 
2ep, where ggB <^ ji and /x ~ ep > 0, with pairing tak- 
ing place in a thin shell around the well-formed Fermi 
surface. In this regime, egs is given by 

^2 



^ > 0; ^ > gsB, 
IJ.<0; \n\ > gsB, 
(6.53) 



Having computed £gs{i^,B) in the regimes of interest, 
the phase diagram is easily deduced by finding B that 
minimizes SGsip-^B), subject to the total atom number 
constraint Eq. (6.48). 



£Gs -c^A^ + ^{S-2ii) + Cy/JIA'^ In ■ 
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(6.54) 



with A = gsB, see Eq. (6.27). 

It is convenient to work with the dimensionless vari- 
ables defined by 

fi= — , A = — , uJo = — . (6.55) 

ep ep ep 

The normalized ground-state energy ecs in the BCS 
regime is then given by 



£GS 



ce 



5/2 



-|-^/AA2 In 



2 

8e-2/t 15^^ ' 



(6.56) 



where, 7^, defined in Eq. (5.19), is a dimensionless mea- 
sure of the Feshbach resonance width Fq to the Fermi 
energy. With this, Eqs. (6.47) and (6.48) become 
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(6.58) 
(6.59) 

, (6.60) 



that admits the normal state (A = 0, /t = 1) and the 
BCS SF state 



A 
4 
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" 3^ 



Se-^^e-''" '('^o-2A)/V7I^ (6.61) 
3/"^ + 2A'^'y-\ (6.62) 



Abcs(A) 
4 



where in the second line we approximately neglected the 
first term on the right side of Eq. (6.60), valid since 
Abcs ^ 1 (and 7^ ^ 1). It is easy to show that the 
BCS solution is always a minimum of £gs{B). 

The meaning of the two terms on the right side of 
Eq. (6.62) is clear once we recall its form in terms of 
dimensionful quantities: 



n : 



lcix'/^ + 2\B\ 



(6.63) 



i.e., the first term simply represents the total unpaired 
atom density, reduced below n since /i < ep, while the 
second term represents the density of atoms bound into 
molecules, i.e., twice the molecular condensate density 
|i?|2. Qualitatively, we see that at large a)o, A <C 1, 
implying from the number equation that /t « 1. 
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b. BEC regime: We next consider the BEC regime 
defined by < 0. As we shall see, in this regime /it < 
and ^ A, so that Eq. (6.53), J(/i, A), applies. This 
yields, for the normalized ground-state energy. 



ec ^ (a>o - 2/<)A 7, 



2„,-l 



A2 



TT .A^ 
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, (6.64) 



and, for the gap and number equations (dividing by an 
overall factor of A in the former) 



~ 2771(^0 -2A) + ^4^ 

- ~ 2'Y~^A^ H . 
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(6.65) 
(6.66) 



In the BEC regime the roles of the two cqiiations are 
reversed, with jl approximately determined by the gap 
equation and A approximately determined by the num- 
ber equation. Thus, jl is well-approximated by neglecting 
the term proportional to A^ in Eq. (6.65), giving 




(6.67) 



At large negative detuning, 3> 1, in other words in 
the BEC regime, Eq. (6.67) reduces to /i « with 
the chemical potential tracking the detuning. 
Inserting Eq. (6.67) into Eq. (6.66) yields 
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Using A = A/ep and the relation A^ = gliib between A 
and the molecular density, we have 



rift 



-lA^n, 



(6.69) 
(6.70) 



which, as expected (given the fermions are nearly absent 
for /i < 0) simply yields rih « n/2 in the asymptotic 
(large |a)o|) BEC regime. 

These results of coiirse match those derived purely on 
the basis of the gap and particle number equations in 
Sec. VI Bl. 



Zero temperature collective excitations and condensate 
depletion 



We would now like to calculate the spectrum of col- 
lective excitations of the BEC-BCS condensate, which is 
contained in the Sfluct part of the effective action. In or- 
der to do that, we need expressions for the self-energies 
Ell and S20 appearing in Eq. (6.20). At zero tempera- 
ture these are given by 



En (g,l^) = f 
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and 



5^20 (q,^,) 



g'^B^ f duj d^k 



I 



27r (27r)3 



[LU- 



(6.72) 



where 

The self-energy En involves an IR divergent integral over 
k. This divergence can be regularized if one notices that 
Ell enters the effective action 5'fluct in the combination 
eo — 2/x -h 2Eii. It is straightforward to check, however, 
that 

eo-2n + 2Eii(0, 0) = 2E2o(0, 0), (6.74) 



by virtue of the saddle-point equation Eq. (6.23). This 
situation is typical in the interacting Bose gas, and 
Eq. (6.74) is nothing but the Goldstone theorem ensur- 
ing that the collective excitations remain massless (also 
referred to as Hugcnholtz-Pines relation [61] in the in- 
teracting Bose gas literature). Therefore, we are really 
interested not in Eii(g, 51), but rather in the linear com- 
bination Eii(q, — Eii(0, 0), which remains finite. 

The spectrum of collective excitations is given by the 

condition that the propagator computed with the help 
of Eq. (6.20) has a pole. To simply the calculations, we 
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will only compute the spectrum at low momentum and 
energy. Following [50] , in anticipation that the collective 
excitations are sound waves and so, Qq q, we expand 
the self-energies in powers of energy and momentum ac- 
cording to 

eo -2/x + 2Sii(g,f]) f« 2Eo + Ein + Sifi2 + nig2 



S2o(g,0) « So + 220^+0252 (6.75) 

where we have used Eq. (6.74). Therefore, the spectrum 
is given by the condition that the determinant of the 
matrix in Eq. (6.20) vanishes. That matrix in our case 
takes the form 



det 



i (ifi (1 + S) + ^ + 2So + Si^^ + nig2) 
So + S2O2 + n2f22 



So + S2O2 + 02^2 
I (-iQ (1 + S) + £^ + 2So + Sifi2 + Hig^) 



0. (6.76) 



r 



This gives for the spectrum fig of excitations 

" (l + S)'+4SoS 
where we introduced the notation 

n = Hi - 2n2, S = Si - 2S2. 



(6.77) 
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In other words, the excitations are indeed sound modes, 
with the speed of sound 



4S, 
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We now evaluate So, S, 11, and S. Doing the frequency 
integral in S2,o(0, 0) gives 
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Differentiating Sii(0, fl) with respect to f2 at f2 = gives 
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Finally, we also find after some algebra 
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In general, evaluation of the integrals in Eqs. (6.80), 
(6.81), (6.82), and (6.83) is straightforward but cumber- 
some. Wc will present results only in the deep BEC and 
BCS regimes (that is, wq < —ep and ^ ep)- 



First, consider the BEC side of the crossover. There 
Wo and /x are negative, and it is clear that at small gs, 

S ^ 1, SqS ^ 1, and mil <^ 1, and thus they can be 
neglected. The speed of sound is then simply given by 



cbec = 



(6.84) 



Since g^B"^ in this BEC regime, we can neglect 

QgB^ in the denominator of Eq. (6.80) to find 
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where we have used 2/z = wo in this BEC regime. There- 
fore, the square of the speed of sound is simply 



2 _ atB^m'/^ 
Cbec — 



(6.86) 



where w n/2. We compare this expression for c^^q 
with a standard expression for a BEC condensates of 
point bosons (see, for example, Ref. [36]) 



AnabB"^ 
(2m)2 
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where at is the boson scattering length (not to be con- 
fused with the scattering length of fermions, given by a in 
Eq. (5.13)). By inspection, we therefore conclude that on 
the BEC side of the crossover, the BCS-BEC paired con- 
densates behaves as an effective gas of weakly repulsive 
bosons with a scattering length 



ab 



gW" 
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(6.88) 



We note that the scattering length ab, together with the 
speed of sound c, decreases as the detuning wo is made 
more negative, deeper into the BEC regime. This is of 
course to be expected, as paired-bosons interaction arises 
due to their polarization into their constituent fermions. 
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followed by a fermions exchange - the process of the order 
of g'^. Since this virtual fermion creation process costs a 
molecular binding energy, Eq. (2.16), deep in the BEC 
regime approximately given by \coo\, it is suppressed with 
increasing |a;o|, as is the effective bosonic interaction and 
at. 

Here an important remark is in order. Our narrow res- 
onance result for Ub, Eq. (6.88), contrasts sharply with 
the well-known (g^ and m independent) result for the 
molecular scattering length deep in the BEC regime, 
namely ttf, w 0.6 a, where a is the fermion scattering 
length [35]. The short answer explaining this difference 
is that the ah f« 0.6 a prediction is for the wide reso- 
nance BEC regime, corresponding to 7s ^ oo instead of 
the limit of narrow resonance considered here. In more 
detail, the results of Ref. [35] apply only in the regime 
where a ^ \ro\- In our narrow resonance problem, this 
regime is realized only in a very narrow range of wq, sat- 
isfying — Fo < Wo < (Fo is the resonance width given 
in Eq. (5.15)). This is not what one should call the BEC 
regime of the narrow resonance crossover, which should 
be defined as wq ^ — ef, with e^, in turn, being much 
bigger than Fq. 

Moreover, if one does tune loq to this narrow window, 
the Fermi energy of the gas under study here will be much 
bigger than the binding energy of the bosonic molecules, 
and the condensate cannot be treated at all as weakly in- 
teracting bosons. So even though the scattering length of 
bosons within this window of loq is indeed 0.6 a, this will 
not get reflected in the speed of sound in the condensate. 

Contrast this with the BEC regime of the broad res- 
onance BCS-BEC superfluid, where a :» ro for a wide 
range of the detuning wq < 0, and where Fermi energy is 
small compared to Fq. For further details, including the 
calculation of the speed of sound in the broad resonance 
BEC regime of Eq. (5.6), see Ref. [6]. 

Thus we conclude that Eq. (6.88) is the correct scatter- 
ing length of molecules in the narrow resonance problem. 
In fact, Eq. (6.88) can also be derived independently by 
studying the scattering of bosons in vacuum perturba- 
tively. We will not do it here. 

Let us now turn to the BCS regime luq > 2ep, where 
/X « ep- We evaluate Eq, S, S, and 11, from Eqs. (6.80), 
(6.81), and (6.82), (6.83). These integrals are easiest to 
compute if we change variables k'^/(2m) — /i = ^, and 
notice that only small ^ essentially contribute to the in- 
tegrals. We find 

_ gs'(2m)3/2 0^0 -2m , . 
^0 ~ — 77-^ — , S ~ — : , (6.89) 
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that reassuringly recovers the well-known result for the 
speed of sound in a neutral BCS superconductor [62, 63]. 

In the intermediate crossover regime between BEC and 
BCS, where < ujq 2ei?, the integrals in Sq, S, S, 
and n should be evaluated numerically to give the speed 
of sound which interpolates between its BEC and BCS 
values. 

Using our understanding of the collective excitations, 
wc can now compute the interaction-driven depletion of 
the condensate, namely the number of bosons that are 
not Bose-condensed into a single-particle k = state, 
even at zero temperature. As wc will show, the depletion 
number turns out to be much smaller than the number of 
particles in the condensate (with the ratio controlled by 



the smallness of 7^ 



5jm 



3/2,-1/2 



), which justifies our 



neglecting it in the analysis of the crossover, above. We 
also note that smallness of depletion justifies the expan- 
sion in powers of fluctuations (controlled by 7^) across 
the whole range of the BCS-BEC crossover in a narrow 
resonance atomic system. 

The number of excited bosons can be simply computed 
from the Green's function of the fluctuations, 
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Evaluating the Green's function gives 



SB(a,-f^) + 2Sn(g,-fl) 



det M 

(6.93) 

where M is the matrix in Eq. (6.20). Evaluating this 
expression analytically in general is difficult, so we con- 
centrate on limiting BEC regime. 

In the BEC regime, where we can approximate E20 ~ 
So, eo — 2fj. + Ell w 2Eo and sum over frequency, we find 
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where 51^ is the spectrum, VLq = 2 

(Eq. (6.77) is the small momentum version of it). Doing 
the integral gives 



^(r.Eo)^/^ = -n,/:^, 



8 



(6.95) 



which coincides with the standard expressions of the con- 
densate depletion in a weakly interacting Bose gas. As 
advertised above, the depletion is small and vanishes in 
the limit of a vanishingly narrow resonance, 7s <g; 1. 

As Wo is increased from negative towards positive val- 
ues, Eq. (6.95) is no longer applicable and analysis is best 
performed numerically. 



The speed of sound, Eq. (6.79), is now dominated by S, 
n, and give 



4.. Critical temperature 



Cbcs 




(6.91) 



We expect the condensate to be reduced with in- 
creasing temperature, vanishing at a critical temperature 
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Tc{luo), that wc compute below. In contrast to superflu- 
ids of point bosons, in paired superfluids there are two 
physicahy distinct effects that contribute to the conden- 
sate reduction with temperature [5]. One is the disso- 
ciation of Cooper pairs (and closed-channel molecules 
hybridized with them), and, simultaneously, thermal 
bosonic excitations. One of them is captured by the fi- 
nite temperature gap equation Eq. (6.16) and is respon- 
sible for Tc in the BCS regime, while the other must be 
included in the finite-temperature particle number equa- 
tion and is at work in the BEG regime. 

For simplicity we focus on the Bose-Fermi mixture at 
the critical temperature, where the condensate density 
vanishes. This allows us to take advantage of technique 
of Ref. [5] to find the number of particles. This method 
ignores the interactions between the bosons and concen- 
trates solely on the bosonic propagator modified by the 
presence of fermions. This amounts to approximation of 
Ss[(f>], Eq. (6.4) by a quadratic expansion in (j){r) = (p{r), 
reducing it to 5'fluct[</^], Eq- (6.20). In the case of a broad 
resonance, considered in Ref. [5], this is not quantita- 
tively justified for wq sufficiently close to zero so that 
a ^ n~^^^. In contrast, in the case of a narrow resonance 
system studied here, this expansion is justified, since the 
strength of interactions, governed by 7s, is weak. 



At T = Tc the condensate vanishes, B 
gap equation reduces to 



9s f d^k 



tanh ^ 



2m 
"fc2" 



0, and the 



(6.96) 
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FIG. 22: Tlic contour of integration in Eq. (6.101). The 
crosses depict the positions of the poles of the Bose- Einstein 
distribution 1/(6"/"^ - 1). 



that is a generalization of the vacuum phase-shift, which 
can be deduced from Eq. (1.7) and the relation f{p) = 
(e2"5 - l) /{2ip). We can now transform the sum over 
frequencies in Eq. (6.98) into the integral 

with the integration over done along the contour de- 
picted in Fig. 22. 

Combining this all together, the particle number equa- 
tion (6.22) 



with given by Eq. (6.25). To find the particle number 
equation, we need to evaluate the contribution to the 
partition function due to fluctuations of (i in Eq. (6.4). 
This can be expressed in terms of polarization operators 
Ell and S20. In fact, for B = 0, S20 = 0, and only En 
survives. This gives for Eq. (6.19) 
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(6.97) 

where cj„ = tiT {2n+l) are fermionic Matsubara frequen- 
cies, with fluctuation corrections to 5*0 [B\ given by 



rffluct 
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(6.98) 



where 



Sn{q, 0„) = zf2„ + ^ - 2/x + eo + 2Sn(g, 0„), (6.99) 

and n„ = 2'KTn are bosonic Matsubara frequencies. To 

simplify this expression further, we can use the technique 
discussed in Ref. [5] . To that end, we introduce the many- 
body finite temperature phase-shift 
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4m 



2Y.u{q,iVt) 

(6T00) 



N 



.dSo 
dfj, 



takes the form 
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with the first and second terms giving the number of 
fermion and bosons, respectively. 

To make further progress, we need to know En, which, 
when evaluated at T = Tc is given by 

y ( n) = -i^S^f— - 

11 W i 2 (27r)3[e+-z(a;„-f2)][e-+ic.. 

(6.104) 

where ^_ are given by Eq. (6.73), above. The sum 
over the frequencies is elementary, with the result 



2 r tanh 
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(6.105) 

Unfortunately the remaining integral can only be done 
numerically, and we will not evaluate it in this paper, 

except deep in the BEC regime. Fortunately, however, 
we do not need to know it in the narrow resonance limit 
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of Qs 0. In this limit, En becomes small and the 
phase-shift reduces to 



6{q, O) = Im In ( f2 — h 2/i — wo 



(6.106) 



In this limit, we can use ujq and eo interchangeably since 
they now coincide. 

Substituting into Eq. (6.103), we can transform the 
contour integral over f2 to the form 



1 



1 



2^ e^/T 



(6.107) 



The contour in this integral can now be transformed to 

2 

enclose the pole at O = ^ going in the clockwise direc- 
tion, finally giving for the total particle density 
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(6.108) 



where we used uoq = 2fj, valid the small gg limit. 

This equation coincides with Eq. (6.11) which wc de- 
rived in the gs = limit, as could have been guessed 
from the outset. However, here we are in principle in the 
position to compute corrections to this equation if Sn is 
evaluated and included in Eq. (6.103). 

Let us now use the gap equation Eq. (6.96) and the par- 
ticle number equation Eq. (6.108) to compute the critical 
temperature as a function of detuning luq. 

In the BCS regime, Wq > ^cp, we expect the transition 
temperature to be exponentially small. As a result, the 
particle number equation forces /x to be very close to 
ep (slightly below it). Indeed, the number of excited 
bosons at a low temperature is expected to be small, 
and the particle number is saturated by fermions, whose 
chemical potential must therefore be in the vicinity of ep. 
We recall that for = 0, Eq. (6.108) would be solved 
simply by setting ^ = ep. 

We then need to use Eq. (6.96), with ep substituted for 
/i with sufficient accuracy to determine Tc- The actual 
calculations arc identical to the ones employed by the 
BCS theory. One technique for solving Eq. (6.96) in this 
regime is described in Ref. [36]. Evaluating the integral 
in Eq. (6.96) we find 
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where C is the Euler constant. In C « 0.577. We see that 
indeed, the critical temperature is exponentially small in 
the ratio (wq — 2ep)/g1. This could have been guessed 
without any calculation as this simply coincides with the 
standard BCS result in the same way as Eq. (6.35) co- 
incides with the appropriate BCS result, with Tc/A = 

In the deep BEC regime, where ljq is negative, wc ex- 
pect the chemical potential fi to roughly follow ujq, in the 



way quite similar to the infinitely narrow resonance limit 
described in section VI A. The critical temperature will 
then be given by solving Eq. (6.108) and noting that the 
fermion part of the particle number is going to be very 
small. Therefore, it will reach its asymptotics coinciding 
with the critical temperature of a non-interacting Bose 
gas, given by Eq. (6.12). 

Between the BEC and BCS regime through the 
crossover the temperature will interpolate between the 
BEC Eq. (6.12) and the BCS Eq. (6.109) values, in the 
precise way that can be obtained through a numerical 
solution. 

An interesting question is whether the critical tempera- 
ture decreases monotonously as the detuning is increased 
or perhaps has a maximum at some intermediate value of 
the detuning. Recall that Nozieres and Schmitt-Rink ob- 
served a maximum in the vs wq diagram, see Rcf. [5] , 
and so did subsequent papers which followed thcnr icvh- 
niques. However, as these authors themselves observed, 
their calculations were done in the case of a broad reso- 
nance, where their approach was an uncontrolled approx- 
imation that could not guarantee that the maximum was 
not an unphysical artifact of their approximation. In con- 
trast, in our case of narrow resonance, we can actually 
calculate the entire curve To(a;o) perturbatively in pow- 
ers of gs , and predict the behavior of in a trustworthy 
way, at least for a 7^ < 1 system. 

For our purpose it is sufficient to concentrate on the 
deep BEC regime where loo <^ —2ep. In this regime we 
expect Tc(ijJq) to approach the limiting value (6.12) of or- 
der ep from either above or below. Given that the high 
Wo BCS asymptotics, Eq. (6.109), is exponentially small 
compared to ep, the approach of the asymptotic BEC 
value (at large negative ljq) from above implies unam- 
biguously that the curve Tc{ujq) must have a maximum 
somewhere. 

In the infinitely narrow resonance case we observe that 
the transition temperature decreases with increasing wq, 
since the fermion number, suppressed in the BEC regime 
as e"^"/^^, would start increasing in accordance with 
Eq. (6.11). However, for a narrow but finite width reso- 
nance, fluctuations must also be taken into account. 

Let us first evaluate the contribution of the fluctua- 
tions to the particle number equation Eq. (6.108). First, 
we compute Eq. (6.105), which in the BEC regime can be 
evaluated and leads to a correction to the particle num- 
ber equation which we now discuss. At ujq < —2ep (and 
consequently, /U < — e^), we can safely neglect the hyper- 
bolic tangents in the numerators of Eq. (6.105) to arrive 
at 
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This expression basically coincides with the correspond- 
ing expression for the polarization operator in a vacuum, 
Eq. (3.26). Physically this is expected since deep in the 
BEC regime there are only exponentially small number 
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However, all this ignores interactions between the 
bosons (in other words, higher terms in the expansion 
in powers of (f> in Eq. (6.4)). It turns out that the inter- 
actions change the critical temperature in a way which 
is not exponentially suppressed in the deep BEC regime. 
A weakly interacting Bose gas with a given scattering 
length (for our case given by Eq. (6.88)) has been exten- 
sively studied in the literature. Although the correction 
to the critical temperature due to interactions is still a 
controversial subject, there is a reasonable agreement in 
the literature that this correction is positive and is given 
by a bosonic gas parameter (see Ref. [51]) 



FIG. 23: In the BEC regime we can deform the contour from 
Fig. 22 to this one, which encloses the pole corresponding to 
the bound state and goes around the cut corresponding to the 
scattering states. 



of fermions, so from the point of view of bosons, the sit- 
uation is indistinguishable from a vacuum. 

We now observe that the phase-shift 6{q, fl) has a sin- 
gularity at 



n--i— + 2u-eo- 2Sii(g, iCt) = 0. 
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The value of fl that solves this equation is given by 
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Here T^q represents the critical temperature in the non- 
interacting limit Eq. (6.12) and n/2 is the density of 
bosons. This expression is clearly much bigger than the 
exponentially small corrections due to fermion number 
and Ell and hence those other corrections can be ne- 
glected. 

Therefore, we conclude that as luq is increased from 
large negative values, Tc{(jJo) actually increases, accord- 
ing to 



n„ = 
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(6.112) 
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To sec this we observe that S only depends on fl and q 

2 

through the combination fl — 4— , and therefore the q de- 

2 

pendence of the solution to Eq. (6.111) is simply f^. At 
the same time, at (7 = 0, the solution to Eq. (6.111) must 
be f2 = 0, owing to the Goldstone theorem Eq. (6.74), 
giving the result (6.112). 

In addition to this pole, the phase-shift S{q,Sl) will 
also have a cut along the real axis of fl, corresponding 

to the scattering fermionic states. This cut goes from 

2 

= ^ — 2/U to infinity (notice that /x < 0). Using this 
information, we can transform the contour in the integral 
over fl in Eq. (6.103) to the one depicted in Fig. 23. 
The integral around the pole gives back the atom number 
confined inside thermally-excited bosons, 
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while the remaining integral along the cut gives terms 
suppressed exponentially as e^^o . These terms represent 
corrections to the particle number equation Eq. (6.108). 
They must be combined with properly evaluated fermion 
number in Eq. (6.108) and with the additional terms 
given by the expansion of the hyperbolic tangent in 
Eq. (6.105) to give corrections to the critical tempera- 
ture in the deep BEC regime wq —2eF- The key ob- 
servation is that all these contributions are exponentially 
small as . 
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where a is an unknown constant of the order of 1. At 

the same time, in the BCS regime, at large positive uq, it 
drops off exponentially, according to Eq. (6.109). Thus, 
Tc{coo) must exhibit a maximum somewhere for the in- 
termediate values of ujo. 

Although this conclusion about the existence of a max- 
imum in rc(ci;o) agrees with those appearing in a number 
of papers devoted to broad resonances, beginning from 
Ref. [5], here, in contrast to those studies our arguments 
for a maximum are robust and quite general, being based 
on a quantitatively trustworthy (in a narrow resonance 
case) calculation, rather than on an uncontrolled approx- 
imations. 

A schematic diagram depicting Tc{luo) is illustrated in 
Fig. 7. To set the proper scale on the vertical axis, we 
notice that the ratio oi ep given by Eq. (6.7) to TbeCi 
Eq. (6.12), is approximately 5. Compare this with the 
critical temperature of an infinitely narrow resonance, 
7s 0, shown on Fig. 21. For a finite resonance width, 
the first qualitative difference is that the critical tem- 
perature is nonzero even for wq > 2e_F, representing the 
BCS regime, absent in the limit of infinitely narrow res- 
onances. Secondly, a Tc{ijJo) for a finite- width resonance 
exhibits a maximum at intermediate values of wq, reflect- 
ing the boson-boson interaction correction to the critical 
temperature in the BEC phase. 



43 



VII. P-WAVE BCS-BEC CROSSOVER AND 
PHASE TRANSITIONS 

A. Coherent-state formulation and saddle-point 
approximation 

In Sec. VB wc established the proper two-channel 
model for a p-wave resonantly interacting atomic gas and 
determined its parameters in terms of results of a two- 
body scattering experiment. We now turn to the study 
of this model at a fixed chemical potential, with the aim 
to establish the phases and phase transitions of such a 
Fermi gas at finite density. 

As usual, the thermodynamics is encoded in the par- 



tition function Z = Ti e^" and the corresponding 
free energy F = — TlnZ. The partition function can 
be conveniently formulated in terms the imaginary-time 
path-integral over coherent states labelled by commuting 
closed-channel fields 0(r), i^(r) (bosonic molecules) and 
anticommuting open-channel fields V'(r), '^(r) (fermionic 
atoms), and their complex conjugates 

Zp = j VtpV4>V(j)V$ e"^", (7.1) 

with the action Sp corresponding to the Hamiltonian 
iJ^-c/i, Eq. (5.34), given by 
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As with the s-wave case in Eq. (6.3) the fermionic atoms can be formally integrated out exactly, to give the effective 
bosonic action 



Sp[<P] = -^Tv \n 
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where the trace is over the 2x2 matrix structure, space 
r and the fermionic Matsubara frequencies Un = 7r(2n -|- 
l)//3. 

The field theory Sp[(j}\ is nonlinear in (j) a-nd there- 
fore cannot be solved exactly. However, as discussed 
in Sec. V B 2 for a narrow resonance it is characterized 
by a dimensionless detuning-independent parameter jp, 
Eq. (1.6), and can therefore be systematically analyzed 
as a perturbativc; expansion in 7p 4; 1. 

A lowest order in this expansion in 7^ corresponds to a 
computation of the function integral over (f) via a saddle- 
point method. The dominant saddle-point configuration 
is a constant (/>a(r) = Ba, that is proportional to the 
condensate of the zero-momentum bosonic operator ac- 
cording to 

B5p,o = ^bp=o- (7.4) 

The resulting saddle-point action then becomes quadratic 

in the fermionic fields and can therefore be easily com- 
puted. Within this approximation it gives the free-energy 
density fp = Sp[B]/{pV) 

fp[B] = ^(e„ -2m)B„B„ - ^ E H^n + El), (7.5) 



where 

E. = /(£-M)%4,^|B.kp (7.6) 

and is also the spectrum of the Bogoliubov quasiparti- 
cles discussed below. In above, B is determined by the 
minimum of /p[B], given by the saddle-point equation 

Sfp[B]/5B„ = 0: 

(e«-2M)i?a=E/S^[B]i?^, (7.7) 

where 

r(^) [Rl «2 / d^k fc„fc/?tanh(^) 

U [B] =9pJ ^ ^ , (7.8) 

obtained by simple contour integration over z = icun 

Above expressions can also be equally easily obtained 
by working within the operator (rather than functional 
integral) formalism, approximating the Hamiltonian by 
a fermionic quadratic form with a variational parame- 
ter B, performing a standard Bogoliubov transformation, 
followed by a trace over decoupled Bogoliubov quasipar- 
ticles and minimizing the resulting free energy over B. 

The complex vector "order parameter" B can be 
uniquely and conveniently decomposed according to 

B = u + iv, (7.9) 
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where u and v are two real vectors. 

For latter use it is important to establish a relation 
between vectors u and v (6 real components) and states 
with a definite angular momentum, characterized by 3 
complex wavefunctions ij(™=0'='=i)_ This connection is 
contained in Eqs. (5.35), (5.36), (5.37). Firstly, we note 
that under a global gauge transformation B — > e"^B, 



B • B = (u^ - w^) + i2u ■ V, 



(7.10) 



transforms as a two-dimensional rank-2 tensor, with real 
and imaginary components rotated into each other by an 
angle 2ip, while 



(7.11) 



is a gaiigc-invariant scalar. Using these 
it can bo shown that as long as u and v 
a phase (p can always be chosen to make 
ular. If u and v are parallel, then they 
and B can be made real by a choice of 
state is equivalent to a state with v = 
that a state characterized hy u = v and 
these properties. 

Using Eqs. (5.35), (5.36), (5.37), and (7.9) we find 



transformations, 
are not parallel, 
them perpcndic- 
remain parallel, 
(fi; hence, a u||v 
0. We also note 
u • V = retains 



S(o) = u, + iv„ (7.12) 
B^^^^ = {u:,TVy)+i{v^±Uy), (7.13) 

which shows that the m = p-wave supcrfluid corre- 
sponds to u||v (equivalently v = 0) pointing along the 
m = quantization axis, and m = ±1 superfluids are 
characterized by states with u _L v, u = f with the pro- 
jection of the angular momentum onto u x v equal to 
±1, respectively. All other u, v states are related to a 
linear combination of above three eigenstates by a gauge 
transformation. 

It is also important to summarize symmetries of the 
free energy /p[B], Eq. (7. 5), (7. 6). Firstly, quite clearly 
/p[B] is invariant under gauge transformations. Secondly, 
in a symmetric case of degenerate m = 0, ±1 Feshbach 
resonances with Cq, = eo, the free energy is also rota- 
tionally invariant. Thus, at a quadratic level /p[B] must 
be a function of the only rotationally, gauge-invariant 
quadratic form, (7.11). At a higher order in B, all terms 
can be expressed as powers of this qiiadratic invariant 
and an independent quartic term |B • Bp, a magnitude- 
squared of the quadratic form in (7.10). In the physi- 
cally interesting case where the rotationally symmetry is 
explicitly broken by distinct e^'s, generically /p[B] will 
not exhibit rotational symmetry. However, within the 
saddle-point approximation, it is easy to sec that the 
first, quadratic term in /p[B] is the only one that breaks 
rotational symmetry, with higher order terms a function 
of the two independent gauge- and rotationally-invariant 
combinations B • B and IB • BP. 



B. Zero-temperature: ground state of a p-wave 
resonant Fermi gas 

1. Saddle-point equation and ground-state energy 

We focus on the case of zero temperature, for which the 
free-energy density reduces to the ground-state energy 
density /j=0[B] =eGs[B] 

SGsiB] = ^(e„ - 2m)5„B„ + ^k), 

(7.14) 

with the saddle-point (gap) equation given by (7.7) and 

kakfj 



1^0 [B] =glj 



(27r)3 ■ 



(7.15) 



It is advantageous at this stage to trade the parameter 
6/3 for a physical detuning {jp, according to 



UJfj 



1+C2 ' 



(7.16) 



introduced in Eq. (5.54). This gives a renormalized 
saddle-point equation 



(Wall + C2) - 2/x) Ba 



/3 



{Iaf3 - CiSo^p) Bf3. (7.17) 



To proceed further, we need to calculate Iaf3, that we 
do in detail in Appendix C. Since the integral is for- 
mally divergent as k^, the leading contribution to lajs 
comes from short scales (high energies), cut off by A 
corresponding to the inverse size of the closed-channel 
molecule. This leading A^ contribution is given by 



I. 



(A3) 
a/3 



Cl Sal3, 



(7.18) 



with Cl defined in Eq. (5.45) by the two-atom p-wave 
scattering calculation, Eq. (5.44), that led to the defini- 
tion of u)a- Hence, as in the s-wave case, this leading 
uv-cutoff dependent contribution identically cancels the 
Cl term in Eq. (7.17), and therefore does not contribute 
to any physical quantity expressed in terms of a physical 
detuning w^. 

However, la/s also has a subleading cutoff-dependent 
contributions that scale linearly with A, and are given 

by 



-'a/3 ~ 



(7.19) 



with the dimensionless constant C2 identical to that de- 
fined by the two-atom scattering theory, Eq. (5.46). 

A tensor /^^ also contains uv-cutoff independent 
low-energy contributions coming from momenta around 
Fermi surface. Because these are infrared divergent at 
B = 0, they are nonanalytic in S = |B|, and therefore 
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(as visual) arc in fact dominant at small B, relevant to 
the positive detuning BCS regime. As detailed in Ap- 
pendix C these contributions are easiest to evaluate in 



J 



the gauge where u and v arc perpendicular, and together 
with Eqs.(7.18), (7.19) in the u • v = gauge finally give 
the explicit gap equation 



(1 + C2) (Wa - 2/x)Ba = -7pC2— ^ [daf3 (u^ + V^) + 2UaUp + 2VaV0] B/s 



5n 



-8/3, 



mg^{u + vy 



2UaUf3 2VaVp 



u{u + v) v{u + v) 



B, 



r 



(7.20) 



Integrating these saddle-point equations over Ba we obtain the ground state energy density 
I 



£GS(U,V) 

1 + C2 



y] {ul + vl) [uja - 2/x + ai In {oq {u + v)}] + ai + a2 



(7.21) 



r 



where 



and 



ai = 

02 = 



27p at \ 



5 1 + C2 n 



7p 



2 1/2 5/2 



3^2 



37r2 



kp = kp/kg 



(7.22) 

(7.23) 
(7.24) 

(7.25) 
(7.26) 



is the dimcnsionlcss p-wavc Fcshbach resonance coupling 
discussed previously. It is straightforward to check that 
Sees/SB^ = gives back Eq. (7.20). 

We emphasize that Eqs.(7.20), (7.21) are written in the 
u • V = gauge. However, once obtained we can utilize 
the gauge-invariance of ecs to reexpress it in an arbitrary 
gauge. To this end we note that u^-l-t;^ = B-B is already 
invariant. However, while v? — is not gauge invariant 
(being a real part of B • B, Eq. (7.10)), its square is a 
gauge-invariant operator written in u • v = gauge, i.e., 
in the u • v = gauge (u^ - v^)"^ = |B • Bp. 

Thus, a gauge invariant form of ecsp] is given by 



£Gg(u,v) 
1 + C2 



= + -2iJ, + ai \n{ao {u + v)}] + ai- 



u + v 



a2 



(B.B)%i|B.B|2 



(7.27) 



where 



i=^B.B + |B.B|, 

1 



V2 



B B - IB B| 



(7.28) 
(7.29) 



A global minimum of the energy density function £gs[B], 
Eq. (7.27) then determines the ground state of a p-wave 

paired supcrfluid at fixed chemical potential, and possible 
quantum phase transitions as a function of detuning and 



chemical potential as the nature of the minimum changes. 



2. Particle number equation 

As discussed earlier in the context of an s-wavc super- 
fluid, for atomic gas experiments of interest to us, it is 
more relevant to determine the ground state at a fixed 
total atom number TV, rather than a chemical potential. 
As usual, however, this problem is related to the fixed 
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H result by supplementing a minimization of eas (the 
gap equation, Eq. (7.20)) with the total atom number 
equation. The latter is given 



^(B|A^|B), 

degs 
dfj, ' 



(7.30) 
(7.31) 



where the right-hand side is the expectation value of the 
total atom number computed in the grand-canonical en- 
semble, i.e., at fixed /x, in the ground state |B) (a BCS- 
type variational one, labelled by B in the case of the 
saddle-point approximation). This gives a relation be- 
tween N and /x, thereby allowing one to eliminate the 
latter in favor of the former. We thus turn to the com- 
putation of the atom number equation. 

Within the above saddle-point approximation (that ig- 
nores molecular field fluctuations) valid at a small 7p, the 
atom number density equation is given by 



n ■■ 



2|B|2+n/, 



where the fermion density is given by 



_i r (fi 

~ 27 (2^^ 



1 - 



fc^/2m - n 



(7.32) 



(7.33) 



The coefficient 1/2 in front of the integral, absent in the 

s-wave case, Eq. (6.24) is due to the fact that here there 
is only a single species of fermions ("polarized" isospin). 
Clearly, according to Eq. (7.31), result (7.32), (7.33) can 
be equivalently obtained by differentiating egs with re- 
spect to 11. 

It is essential to note a crucial qualitative difference be- 
tween Eq. (7.33) and its s-wave counterpart Eq. (6.24). 
For qIB'^ < jjb, the s-wave fermion density Eq. (6.24) 
at nonzero B can be estimated by the density p^/(37r^) 
[Pn = v'2m/x) of a degenerate noninteractive fermion gas 
at the same chemical potential. However, because in the 
p-wavc case, for B ^ the occupation number nj{k) (in- 
tegrand in Eq. (7.33)) exhibits a long tail, the integral 
in Eq. (7.33) is formally linearly divergent at large mo- 
menta, cutoff only by the inverse closed-channel molecu- 
lar size A [42, 43]. To compute the fermion number, we 
separate out this large B-dependent short-scale contri- 
bution, finding 



n/ = no/ -I- 2c2|B| 



(7.34) 



where C2 = m^gpA/Sn'^ is the dimensionlcss parameter 
that already appeared in the two-body study of the p- 
wave two-channel model. Sec. VB, see [45], and 



(2^ 



3A:2 



B| 



(7.35) 

is a remaining contribution to nj that is uv-convergent, 
i.e., is not dominated by large momenta, and as a result 



for g(p|Bp <C /U, can be estimated by its Qp = value 



3/2 



67r2 



9{fi), for3^|B|^</z, (7.36) 



where 9{p) is the usual theta-function, equal to 1 for 
positive argument and to for negative argument. 

In the range of detuning where n < g^lBj^, the full 
integral in Eq. (7.35) must be computed more precisely, 
but this is a very narrow range of the chemical potential 
and can (and will) be ignored. 

Thus we find that the atom number-density equation 
is given by 



2(1 + C2)|B| 



no/, 



(7.37) 



to be contrasted with its s-wave analog, Eq. (6.24). As 
noted above the number equation can be directly ob- 
tained from SgSj Eq- (7-21), via Eq. (7.31) and in partic- 
ular the key enhancement factor (I-I-C2) above arises from 
the same factor in sas- Its implication depends on C2. If 
C2 <C 1 , then the number equation is no different than its 
s-wave counterpart and for example in the BEC regime, 
where jj, < the total atom number is "carried" by the 
bosons. If, however, C2 ^ 1,[45] then it shows that even 
deep in the BEC regime, where /i is large and negative 
and correspondingly no / is vanishingly small, the density 
of bosons is given by n/(2c2) and is a small fraction of 
the total atom density, n. In this case the total atom 
number is in a form of free atoms with density given by 
2c2|B|2, the last term in Eq. (7.34). This is a reflection 
of the fact that the p-wave interactions (proportional to 
k'^, due to a centrifugal barrier diverging at short scales) 
are strong at large momenta and therefore for large C2 
lead to a large depletion of the molecular condensate, 
even in the BEC regime where fermions are at a negative 
chemical potential. This is a phenomenon not previously 
discussed in the literature. 



3. Phases and phase transitions of the p-wave BGS-BEC 
superfluid 

Zero-temperature phases, crossover and transitions as 
a function of detuning in a p-wave resonant Fermi gas 
are completely encoded inside the ground state energy 
function eG's[B,/i], (7.21) or, equivalently the associated 
gap and number equations, Eqs.(7.20),(7.37). From our 
earlier analysis of the s-wave BCS-BEC crossover for 
a narrow resonance in Sec. VI A, we can already antici- 
pate some of the qualitative phenomenology associated 
with changing of the detuning. At zero temperature the 
gas will condense into a j>wave superfluid that at large 
positive too will be of a BCS type with weakly-paired, 
strongly-overlapping Cooper pairs and correspondingly 
an exponentially small boson number. As u>o is lowered 
past 2eF, the number of bosons in the condensate will 
grow as a power of 2e^— cjq, while the number of fermions 
will diminish, reflected in the tracking of the chemical 



47 



potential with detuning, fi « luo/2. This intermediate 
crossover regime will thereby consist of a superposition of 
small (of size A^^) closed-channel molecular bosons and 
much larger Cooper pairs. Finally, for loq lowered below 
zero, the tracking chemical potential will change sign to 
/U < and (for small C2) the condensate will transform 
into a purely molecular Bose-Einstein condensate. 

Although very generally this picture remains correct, 
there are a number of qualitatively important differences 
in evolution with the detuning between s-wave and p- 
wave superfluids. Firstly, p-wave superfluid is character- 
ized by a richer complex vector order parameter B, asso- 
ciated with £ = 1 angular momentum of the condensing 
boson, and therefore admits a possibility of a variety of 
distinct p-wave superfluid ground states and associated 
quantum phase transitions between them. Possible su- 
perfluid ground states are distinguished by a projection 
of condensate's angular momentum along a quantization 
axis. This allows for a possibilities of a time-reversal 
breaking to = 1 states (and its rotated and time-reversed 
versions) referred to as apx+iPj/-superfluid with a projec- 
tion of the angular momentum of the condensed bosons 
onto the z-axis equal to +1, or a p^-superfluid (and its 
rotated analog), with a projection of the condensate's an- 
gular momentum onto the 2;-axis equal to 0. As discussed 
in Sec. VII A these two phases are characterized by u and 
V, defined in Eq. (7.9), with the u = v, u -L v state corre- 
sponding to px + ipy phase, and v = state the pz phase 
(which we will often refer to as px phase as well, [64]), 
respectively. It is also possible to have a "superposition" 
phase, where u ± v, but with unequal lengths, u 7^ v, 
corresponding to a time-reversal breaking state in which 
all bosons condense into a linear combination of Pz and 
Px + 'iPy orbitals. 

Secondly, and related to above, a p-wave gas is char- 
acterized by (potentially) three distinct detunings, Wa, 
one for each component of the £ = 1-field As dis- 
cussed by Ticknor, et al. [37], in systems of interest to 
us, this resonance splitting, 6 arises due to the inter- 
atomic dipolar interaction predominately due to electron 
spin. Although it is rotationally invariant in the spin- 
singict closed-channel, the source of anisotropy is a small 
admixture of the spin-triplet channel, with a result that, 
with the quantization axis along the external magnetic 
field H (that we take to be along x), the m = res- 
onance is lower by energy S > than the degenerate 
m = ±1 doublet. Thus we will take 

OJx = Wo, 
^y.z = ^0 + (7.38) 

This feature will be key to a nontrivial phase diagram 
possibilities illustrated in Figs. 8, 9, 10, 11. 

Finally, another important difference that has already 
been noted in the previous subsection is the large C2 limit 
of the p-wave number equation, (7.37), in which even 
for /i < (in what one would normally call the BEC 
regime) the fermion density is large and correspondingly 



the boson density ]B|^ Ri n/(2c2) is vanishingly small for 
C2 > 1 [45]. 

To determine which of the p-wave superfluid phases is 
realized by the BCS-BEC condensate, we minimize the 
ground-state energy ecs, Eq- (7.21) with respect to u 
and V for iVa of interest, while enforcing the total atom 
number-density constraint Eq. (7.37). 

a. Isotropic p-wave Feshbach resonance We first 
consider a simpler isotropic case, where LUa = <-^o for all a. 
Utilizing the rotational invariance of eas, it sufficient to 
minimize it over magnitudes u and v. Analogous to other 
isotropic problems with a vector order parameter (e.g., 
a Heisenberg magnet), the actual global (as opposed to 
their relative) direction of vectors u, v in the ordered 
phase will be chosen spontaneously. 

Although ultimately we need to minimize eGs[u,v] at 
fixed total atom number, i.e., subject to the atom number 
equation contraint 

2{1 + C2){u'^ +v'^)+nof =n. (7.39) 

it is important to first study £(55 [u,?;] at fixed ^. Stan- 
dard analysis of ecs, EQ- (''•21) shows that there are four 
extrema: (i) u = v = (normal state), (ii) u ^ 0, v = 
(paj-superfluid state), (iii) u — 0, v =/= (paj-supcrfluid 
state), and (iv) u = v {px + ipy-superfluid state), 
where clearly (ii) and (iii) correspond to the same su- 
perfluid state. [64] After some standard algebra, one can 
show that the normal (u = v — 0) state is always a max- 
imum with energy £gs[0, 0] = Eq^ = 0. 

The nature and relative stability of the other extrema 
is decided by the parts of sgs[u,v] that do not depend 
on the + combination, namely by terms 

ai I {u^ -h v") ln(« + ^;) + +'^{u^~ v")" . 

\ U-\- V J I 

As illustrated in Fig. 24 and standard analysis shows that 
M ^ 0,u = and u = 0,t; extrema are degenerate 
(guaranteed by w ^ symmetry) saddle-points and u = 
w ^ is a global minimum, independent of the actual 
values of a\ and 02, as long as they are positive. In the 
BEC regime, where a\ = 0, this is clear since the last 
02 term prefers the u ~ v state, but also remains true 
throughout the BCS and the crossover regimes. 

Because there is only one local minimum, this fixed 
chemical potential result automatically applies to the 
minimization of ecs at fixed total atom density n, with 
fi eliminated through Eq. (7.39). Thus we conclude that 
at T = 0, the ground state of a Fermi gas interacting 
with an isotropic Feshbach resonance is a. Px + ipy-wave 
superfluid throughout the BEC-BCS crossover. 

In the BCS regime, this results agrees with the well- 
known prediction by Anderson and Morel [52], who 
showed (in the context of Ai-phasc of •^Hc) that a polar- 
ized p-wave BCS superconductor at T = is always in the 
Px + iPy state. Thus, our above conclusion extends their 
result to the BEC (p < 0) and crossover {0 < fj, < ep) 
regimes of a resonantly-paired superfluid. 



48 



V 








u 



FIG. 24: A contour plot of eas{u,v) in the absence of split- 
ting, S — 0. The global minimum -At u = v ^ 0, saddle-points 
at u = 0,v ^ and u ^ 0, v = 0, and a maximum at 
u = V = can clearly be seen. 



Wc now compute this px + ipy ground state explicitly. 
To this end, we substitute u = v into the ground state 
energy, Eq. (7.21) to obtain s^g^^^[u] = eGs[u,u] 



,+ipy 



GS 



1 + C2 



2u^ [cjo — 2/i + ai ln(2aou)] + aiu^ + 4a2W^ 



(7.40) 

and minimize it with the constraint Eq. (7.39). At fixed 
/U, the saddle-point equation ds^g^''^[u]/du = is given 
by 

u[u!o - 2/i + ai + ai ln(2ao'u)] + 4a2U^ = 0. (7.41) 

As in the s-wave case, once the atom number constraint 
is implemented, the detailed behavior is quite different 
in three regimes, depending on the range of detuning uq- 

BCS regime: 

For LVo > 2eF, closed-channel molecules (b) and the 
corresponding condensate are energetically costly lead- 
ing to a small u. This justifies us to neglect the molecu- 
lar contribution (first term) inside the number equation 
(7.39). Then, with approximation of no/(At) by the nor- 
mal state atom density (i.e., also ignoring the small con- 
densate density there) (7.36) immediately gives ^ « ep- 
Furthermore, similarly neglecting the subdominant quar- 
tic term, a2u'^ inside e^g^''^ {6 = 0), Eq. (7.40), the cor- 
responding saddle-point equation can then be solved an- 
alytically, giving 



5=0 



1 



2aoe 



-{u>o-2€F)/ai 



(7.42) 



= g-ll/6 /f^ g-(a;o-2ei.)(l-F7pA/feF)/(27pejr) 



for coq > 2eF, (5 = 0, 

that is indeed exponentially small in this BCS 
regime. The corresponding condensation energy density 



c+iPy 

'GS 



^P^+iPy 

'GS 



Px+ip, 



, i5] is given by 



45^*"" (0) = -(l + C2)al«+,,J^ (7.43) 

= _4e-lV3 g-(c^o-2ci,)(l+7pA/fei.)/(7p€F) 

for Wo > 2eFj (5 = 0, 

Within the same set of approximations it is also 
straightforward to compute the corresponding quantities 
for the Px-state, obtaining {vP = u^^°) 



< = 



3/2 



-{uo-2eF)/ai 



(7.44) 



V > 

for Luo > 2eF, S = 0, (7.45) 

that gives a ratio n{S = 0) = egs^'^" (0) /eGs(O) = e/2 of 
condensation energies for the two states, consistent with 
the numerical value reported in Ref. [52] and thereby con- 
firms that Px + ipy state is energetically more favorable. 
Crossover and EEC regimes: 

For ojq < 2ep, it becomes favorable (even in Qp — > 
limit) to convert a finite fraction of the Fermi sea 
(between and 2eF) into a BEC of closed-channel 
molecules. Consistent with this, the log contribution in 
Eq. (7.40) is no longer large, with s^g''^^ immediately 
giving a chemical potential that tracks the detuning ac- 
cording to /i w cjo/2 with accuracy of O^-jp). 

As previously noted [28, 65], we observe that the 
roles of number and gap equations interchange in the 
LUo < 2eF regime, with the former determining the molec- 
ular condensate density and the latter giving the chemi- 
cal potential. Consistent with this, the number equation, 
Eq. (7.39) then gives the growth of the bosonic conden- 
sate according to 



2 , 2 
U + V 



2(1 + C2) 



, 3/2 



reaching a maximum value of 



2 , 2 
U + V 



2(1 + C2)- 



(7.46) 



(7.47) 



for a negative detuning. As already noted above, it is 
remarkable that even for a large negative detuning the 

boson density never reaches its (ideal, 7p <C 1) maximum 
value of n/2 corresponding to the total atom density n. 
Instead, due to p-wave interaction that is strong at short 
scales, C2 > 0, and a p-wave molecular condensate is 
depleted into open-channel atoms. We will nevertheless 
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continue to refer to this range of detuning as the BEC 
regime. 

Since as shown above, in the crossover and BEC 
regimes the ground state remains a. Px + ipywave su- 
perfluid the p-wave order parameter is given by 



,1/2 



2(l+C2)l/2 



3/2 



UJo_ 

for uo<2eF. (7.48) 



1/2 



b. Anisotropic p-wave Feshbach resonance We now 
analyze the more experimentally relevant anisotropic 
case [37], where the triplet Feshbach resonance is split 
by dipolar interactions into a m = ±1 degenerate dou- 
blet resonance and an m = resonance, with LOa given 
by Eq. (7.38). With the magnetic field H picking out 
a special direction (that we take to be x), the ground- 
state energy function eas is no longer rotationally invari- 
ant. Within our saddle-point approximation this uniaxial 
anisotropy only enters through the detuning part 



ujq{u^ + v^) + 5{u^y + ul +vl + vl). (7.49) 



With S > 0, this uniaxial single-particle energy is clearly 



minimized by Uy = Uz 



0, i.e., when u 



and V are parallel and point along H = _ffx, correspond- 
ing to the pa;-wave ground state. In our more convenient 
u- V = gauge choice, this Px state is equivalent to either 
u or V pointing along H and with the other vanishing. 
Furthermore, in this transverse gauge for SiPx + ipy state 
[64] (that, as we saw above is preferred by the interac- 
tions) in which neither u nor v vanish, e'Qg°'^' is clearly 
minimized by choosing the longer of the u and v to be 
along H = _ffx, while the shorter one spontaneously se- 
lects a direction anywhere in the (yz-) plane perpendic- 
ular to H. For u = V, their overall orientation is chosen 
spontaneously. 

An explicit minimization over the direction of u-v or- 
thogonal set confirms these arguments, giving 



~-GS 



[u, v] = loq{u^ + v^) + 5 Min[?i^, 



(7.50) 



It is convenient to take advantage of the exchange sym- 
metry u <-^- V, £gs[wj v] = eGs[v^ m], and ior v (with- 
out loss of generality) always choose u to be the longer 
vector, with the other state physically equivalent. With 
this choice and Eq. (7.50) the ground-state energy is min- 
imized by u directed along H. The resulting ground state 
energy as a function of magnitudes u and with u > v 
and u = MX takes the form 



eGs[u,v\ = cD(u2 -hi)^) -h (w^ + i)^)^ + -(u^ --0^)^ +<5i)^ + ai 



(m^ + 1)^) ln(u -h v) + 



u + v 

for u> V, vi = wk, u _L V. 



(7.51) 



where to simplify notation we introduced dimensionless 
variables 



w = 


(wo -2ii + ai ln(ao\/n)) /{a2n), 


(7.52) 


5 = 


6 


(7.53) 


ai = 


ai 

a2n' 


(7.54) 






(7.55) 


(1 -1- 02)02^2 £0 


;/ = 


U ^ V 

v=-=. 


(7.56) 



As in the isotropic case above, the ground state as a 
function of detuning LiJq and dipolar splitting 5 is found by 
minimizing egs [U) v] over magnitudes u and v with the 
constraint of the total atom density equation, Eq. (7.39) 

Standard analysis shows that eGs[u,v\ generically has 

three physically distinct [64] (confined to u > t) > quad- 
rant) extrema: (\) u = i) = (normal state), (ii) u > 0, 



^ = (Px-superfluid state), and (iii) u> v > Q {Px+iPy- 

superfluid state [64]). At zero temperature for a fixed 
atom density the normal state is always a maximum with 
energy £(35 [0,0] = Eq^ = 0. However, in contrast to the 
isotropic case, here the relative stability of the px and 
Px + iPy states crucially depends on the detuning uq and 
dipolar splitting 6. This is summarized by the contour 
plots of eGs{u,v) for (5 = (Fig. 24), S small (Fig. 25) 
and S large (Fig. 26). We now study these in detail. 
BCS regime: 

As a generic property of the BCS regime, for lvq > 
2eF, the molecules are energetically suppressed, and only 
exponentially small condensate is expected. The number 
equation then leads to p, ^ ep. It also allows us to neglect 
the subdominant quartic in u and v contributions inside 
£gs, Eq. (7.51), allowing the corresponding saddle-point 
equation to be solved analytically. For the two candidate 
p-wave superfluid states we find: 

Up.+iPy = (1 + V«i) e-'/'"^<+,p^, (7.58) 

= -^(l+(5/ai) e-('^°-2^^+''/2)/«S 
2aoe 

v,^+,,^ = {l-6/a,)e-'/'^^vl^,^^, 
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FIG. 25: A contour plot of eGs{u,v) in the presence of a 
small splitting S. The global minimum at u ^ v ^ and 
saddle-points aX u = 0,v ^ and u ^ 0, v = 0, can clearly be 
seen. 
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FIG. 26: A contour plot of easiu, v) in the presence of a large 
splitting 5. Only the global minima at w = 0, ^ and u ^ 0, 

V = are present 
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-(wo-2€F+<5/2)/ai 

2aoe^ 

for Wo > 2eir, 



and 



(7.59) 



1 



0, 



-(wo-2eF)/ai 



^^5(0), 

for Wo > 26^?, 



where the corresponding ground-state energies at finite 
splitting b have been expressed in terms of i5 = energies, 
Eqs.(7.43),(7.44). The ratio oi'Px + Wy and "Px ground- 
state energies is then given by 



-GS 



1 + 



-S/ai 



(7.60) 



and reduces to the previously found result of e/2 for ^ = 
0. 

Consistent with the analysis of the ^Hsotropic reso- 
nance" subsection, for low dipolar splitting 6, TZ{6) > 1 
and in the BCS regime the Px + ipy superfluid [64] is the 
ground state, as seen on Fig. 25. However, TZ{d) reaches 
1 at 6^*^^, given by 



1+C2 



€f, 



(7.61) 



(7.62) 



signaling a quantum phase transition from the px + ipy 
to Px ground state for 5 > 6c [64]. This is consistent with 
intrinsically positive quantity Vp^^ip^ ^{S), in Eq. (7.58) 
turning negative (unphysical) for 6 > S^'"^' ■ 
BEC regime: 

We can similarly evaluate the order parameters, 
ground state energies and the Px-Px + iPy quantum phase 
transition boundary in the opposite, BEC regime of a 
large negative detuning luq and < 0, which reduces the 
number equation, Eq. (7.57) to 



U +V = 



1 



2(1 + C2) 



(7.63) 



As discussed for the isotropic resonance case, in the BEC 
regime the condensate is no longer exponentially small 
(given by a finite fraction of total atom density, as seen 
above), and as a result we can neglect the Si terms in 
Eq. (7.51) for small ai. Standard minimization of the re- 
sulting ground-state energy function, together with the 
number equation (7.63) gives for two cxtrema, one cor- 
responding to a.px + ipy superfluid [64] 



•Px+iPy 



l=n'/^(hB+S/2y \ (7.64) 

5(1 +C2) s y/' 



1 nl/2 



= - n 



+ 



Px+iPy 



2 \1 + C2 8C27p CF 
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= 2^ 



1/2 



1 _ 5(1 + C2) S 

1+ C2 8C27p €F 



1/2 



JPx+iPy 

'GS 



{6) = -eo{ni + S'/8), 



2c27p , 5(1 + C2)2 52 



+ 



epn, 



,5(1 + C2)2 ' 8c27p 
and one corresponding to a Px superfluid 

= (7.65) 



Wb, = n ' n 



V2(l + C2) 



3 C27p 



5(1 + C2)^ 



epn. 



As argued earlier for the isotropic case, for low dipo- 
lar splitting 6 the ground state is a Px + ipy superfluid, 
with order parameter and ground-state energy given in 
Eq. (7.64), illustrated on Fig. 25. As wc can sec from 
the form of Vp^+ip^ this minimum and the corresponding 

state disappears for 5 > 2fiB = 5f^^, which gives the 
critical splitting 



^BEC ^ 2a2nnB, 
_ 8 C27p 



(7.66) 



5(1 + C2)^ 



for the quantum phase transition from px + ipy to Px 
superfluid [64]. 

The behavior of the p-wave superfluid order parame- 
ters and ground state energy as a function of splitting 



5 and for full range of detuning luq is best mapped out 
numerically and gives a smooth interpolation between 
above extreme (BCS and BEC) limits derived above. 
However, the phase boundary 5c (wo) for the quantum 
phase transition between px + ipy and Px superfluids can 
in fact be obtained analytically. 

To this end we start at a large dipolar splitting, for 
which the p^^-superfluid (m > v = 0) is a stable ground 
state and therefore the eigenvalues of the curvature ma- 
trix of eGs[u, v] are positive in this state. We then locate 
the critical phase boundary (5c (wg) by a point where the 
eigenvalue along v direction changes sign, becoming neg- 
ative and therefore signaling an instability toward devel- 
opment of a finite value of v characteristic of the Px + ipy- 
superfluid. 

To carry this out, we first minimize eGs[u,v] to im- 
plicitly determine the value of Up^ (with Vp^ = 0), that 
is given by: 



(2(2) -|- 3ai) -|- QUp^ + 2di Inwp^ = 0. 



(7.67) 



Although above saddle-point equation cannot be explic- 
itly solved for Up^ , it can be used to evaluate the eigen- 
values of the curvature matrix at the Px minimum, and 
thereby determine the transition boundary Sc{(-Oo). Com- 
puting the eigenvalues of the curvature matrix of the 
ground-state energy at the Px minimum we find that Px 
superfluid is stable for 



S-ai- 2ul^ > 0, 



(7.68) 



which when combined with the atom number equation 
(7.57) gives (to lowest order in jp) 



S{u)o) « ai{n) + 2ns 1 - 




Sf""' = 4 , for u;o > 2eF 

/ \ / \3/2 (7.69) 



for the (dimensionless) critical boundary illustrated in two values we predict a continuous quantum phase tran- 
Fig. 27, with the system transitioning into the Px + ipy sition at a critical value of detuning, given by (to C(7p)) 
superfluid [64] for 6 < (5c('^o)- In above we used small 
7p (narrow Feshbach resonance) approximation for the 
chemical potential i^{uio) (derived above) appropriate for 

different regimes. As anticipated the phase boundary / (5 S^^^ \'^^^ 

Sci^o) smoothly interpolates as a function of detuning '*^o('^) ~ "^^f ( ~r7=w — rfc^ I • (7.70) 

between the BCS and BEC results found in Eqs.(7.61), V^c - J 

(7.66). Since for all values of the dimensionless coupling 

C2 = ^pA/kp, S^^'^ < 5^'-'^, for 5 falfing between these 
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FIG. 27: The phase diagram of an anisotropic p-wavc super- 
fluid at zero temperature, ilhistrating a phase boundary of 
(dimensionless) dipolar sphtting 5c{(jJo) as function of detun- 
ing, that marks a phase boundary of a continuous quantum 
phase transition between a px- and px + ipy-superfluid. 



C. Finite temperature: phases and transitions in a 
p-wave resonant Fermi gas 

We now extend our study of the phase behavior of a 
p-wave resonant gas to finite temperature. This involves 
a calculation of the frcc-cncrgy density, /p[B], Eq. (7.5), 
and its minimization along the lines similar to the above 
T = analysis of the ground-state energy density ecsfB]. 
The former amounts to a computation of the polarization 
tensor /^'^^[B], Eq. (7.15), details of which we relegate to 
Appendix C. 

The upshot of detailed calculations, presented in the 
Appendix C, is that (as usual) at finite T the low-energy 
singularities arising from Fcrmi-surfacc low-energy con- 
tributions to /a/3[B] are cutoff by T. Consequently, (in 
contrast to the T = case, above), the free energy, /p[B] 
is an analytic function of B, that at high temperatures, 
where |B| is small is Taylor-expandable in powers of the 



J 



gauge- invariant tensor BaBp. Naturally, in the isotropic 
case LVa = t^o a-nd /p[B] only involves rotationally invari- 
ant traces of the powers of a tensor Qa0 = SaB^ and its 
transpose Qpa- 

Within the saddle-point approximation only the 
quadratic contribution is anisotropic, and the resulting 
free-energy density is given by a standard Landau form 

fp[B]= J2 ta|Sa|' + Ai(B-B)' + A2|B-Bp, 

a=x,y,z 

(7.71) 

where, because of the dipolar-anisotropy splitting 
tx{T,LOo) < ty{T,LOo) = tz{T,LOo) = t±{T,LOo), reflecting 
uniaxial symmetry of the system, and these parameters 
vanish linearly at respective T^s, with 



ix(T,wo) 
t±{T,uJo) 



(7.72) 



and 



T^{uo) > T.^icoo). (7.73) 
The parameters Ai,2 are only weakly temperature depen- 
dent. 

Beyond the saddle-point approximation, we expect 
that generically only the gauge-invariance is preserved 
by all the terms in /j,[B] and lack of rotational symme- 
try for finite 6 will be reflected by all terms. However, 
for our purposes it will be sufficient to keep only the 
dominant non-rotational invariant contribution entering 
through the quadratic term as reflected in /p[B] above. 
In terms of u and v parametrization the free-energy den- 
sity is given by 



/,[u,v]= U<+vl) + X,{\u\'' + \v\y + X,{{\u\'-\v\Y+4{u-vf), 



(7.74) 



r 



where the ratio of 4 between the two A2 terms is a 
generic feature that is a reflection of the underlying 
gauge-invariance. 

1. Isotropic 

In the isotropic case {6 = 0), = t± = t{T), and the 
free-energy density is fully rotationally invariant, given 

by 

4^°[B] = f|B|2 + Ai (B • B)' + A2IB • B|2. (7.75) 

For i > (T > Tc), /^"°[B] is minimized by B = and 
the gas is in its normal (nonsuperfluid) phase. Upon low- 



ering T below Tc, a minimum develops at a finite value of 
B. As can be seen from the its expression in terms of u 
and V, Eq. (7.77) the minimum is at m = t; and u • v = 
(or any of its gauge-equivalent states corresponding to 
unequal and nontransverse u and v). Thus, the finite-T 
normal-to-superfluid transition is to a p^; -|- ipy-superfluid 
(SFp^+jpj^), consistent with our earlier finding [52] that 
the ground state is a p^; -I- ipj^-supcrfluid, for all detun- 
ings. At this transition the global U{1) gauge-symmetry 
is spontaneously broken, corresponding to a choice of a 
phase of B or equivalently the relative orientation and 
magnitudes of u and v (as long as they arc not par- 
allel or one of them does not vanish, since this would 
correspond to a Px state that is not connected to the 
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Px + iPy state by a gauge transformation) . In addition, 
an arbitrary choice of an overall orientation of B (i.e., 
of the u — V frame, that by gauge-choice can be taken 
to be orthogonal) spontaneously breaks 0(3) rotational 
symmetry. Clearly time-reversal symmetry is also spon- 
taneously broken in the px + «Pj,-superfluid state. 

This finite-temperature transition is in the complex 
0(3) universality class, which can be thought of as a 
well-explored real 0(6) model[66], explicitly broken by A2 
crystal symmetry-like breaking fields, analogous to 0(3) 
ferromagnet in a crystal-fields due to spin-orbit coupling 
to a lattice. Its critical behavior has been extensively 
explored by Vicari, et al. [67]. 



2. Anisotropic 



We now turn to the more experimentally relevant uni- 
axially anisotropic case, of a Feshbach-resonance triplet 
split by 5 > (as described above) by dipolar interactions 
in the presence of an external magnetic field H = iJx. 
The dipolar splitting considerably enriches the phase di- 
agram, allowing for three possible phase diagram topolo- 
gies, illustrated in Figs 9, 10 and 11. In terms of the 
complex 0(3) model dipolar-splitting leads to an easy- 
axis (Ising) anisotropy, with the free-energy density given 

by 



janisot 



= t,\B,f + tj.|Bi|2 + Al (B • B)' + A2IB • B|2, 

= t,{ul + vl) + t^{ul + vl) + Ai(|up + IvH^ + A2 ((|up - IvH^ + 4(u ■ v)2) 

I 



(7.76) 
(7.77) 



where ± indicates two components in the plane perpen- 
dicular to the external magnetic field H axis that we have 
taken to be x. For tx < t± it is clear that Bx part of B 
will order first, with Bj^ = 0. Namely, since > T^, 
u and V will always both order parallel to the a;-axis, 
showing that for arbitrary small splitting 6 > and ar- 
bitrary detuning ujq. the finite temperature normal to p- 
wavc superfiuid transition is always to the y»a;-superfiuid 
SFp^. We designate this upper-critical temperature by 
Tc2{i-^a) and expect it to be set (up to renormalization 
by fiuctuations) by T^(ti'o)- Clearly from the structure 
of /«"'«°*[B], Eq. (7.76), the noncritical ("massive") B_l 
component can be safely integrated out at the N-SFp^ 
transition, leaving a Landau model of a single complex 
order parameter Bx- Hence the finite-T N-SFp^ classical 
transition is in 3D XY universality class, at which only 
a global U{1) gauge symmetry is broken. 

What follows upon further lower the temperature qual- 
itatively depends on the strength of the dipolar splitting 



6. This follows from the zero-temperature analysis of Sec. 
VII B 3 and is summarized by phase diagrams in Figs 9, 
10 and 11. 

For weak (normalized) Feshbach resonance dipolar 

splitting < S < Sf^'^ , upon further lowering tem- 
perature from a pa;-superfluid phase, the system always 
undergoes a transition to a -I- ipj,-superfluid for all 
detuning loq; we designate this critical temperature by 
Tci{luo)- To see this, we observe that for this low range 
of 6, the parameter t_i_ becomes negatives with reduced 
T and thereby leads to another critical temperature at 
which the B_l component also orders. This ordering takes 
place in the presence of a finite px order parameter Bxo, 
within mean-field theory given by 



Brr 



~tx 



Y 2(Ai+A2) 

The resulting Landau theory for Bj_ is then given by 



(7.78) 



J 



fp^^p^+ip^ = {t^ + 2Ai|i?,on|B^p + (Al + A2)|B^|4 + A2 {BloBs. ■ B^ + B^^B^ ■ B^) 



{i±_ + 2Acos(2(^_L - 2ipxfi)^ 



Bi|2 + A|Bj 



(7.79) 
(7.80) 



r 



where A = Ai -|- A2 and 



A = A2|Bj;o| 



A2 



2(Ai + A2) 
i±_ = t_L + 2\i\Bxo^ = t_L + 



tx I ) 



Al 



Al +\2 



(7.81) 
\tx\- (7.82) 



We note that at this transition the phase tf± of B_l locks 
to the phase (pxfi of Bx^ so that the relative phase is 
±7r/2. This is exactly what is expected upon ordering 
into one of the two degenerate Px ± ipy states. We thus 
find that the SFj,^ to SFp^+jp^^ transition is modified by 
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FIG. 28: Finitc-tcmporaturc phase diagram illustrating con- 
tinuous transitions between normal (N), pa;-supcrfluid (SFp^) 
and Pa; + ipy-superfluid (SFp^+ip^). The parameters txiT) < 
t±{T) are reduced temperatures split by 6. Only tx < t± part 
of the figure is physically relevant. 



the presence of px order and takes place at 



t±{T) 



2(A2-Ai)|B,o| 
A2 — Ai 



X2+X1 



(7.83) 



which then in turn determines Tci{uJo). Since the U{1) 
gauge-symmetry is already broken in the p2;-supcrfluid 
phase and since, as seen above ip± is automatically locked 
to (fix.o, the remaining symmetries that are broken at 
this transition are the 0(2) rotations of Bj^ about the x- 
axis (set by the magnetic field H) and the time-reversal 
symmetry associated with a choice of one of the locking 
angles ±7r/2, corresponding to angular momentum pro- 
jection m = ±1. Thus the SFp^ to SFp^^ip^ transition is 
also in the well-studied 3D XY universality class. Above 
results are summarized by a finite temperature part of 
the phase diagram, illustrated in Fig. 28 

At intermediate dipolar splitting Jf^*^ < ^ < 5^'^^, 
the N-SFp^, transition can also be followed by the SFp^, to 
+ipy transition, but only for dctmiing luq > luoc{T), 
as illustrated in Fig .10. The zero-temperature critical 
frequency wqc (with Hmits woc(<5 S^^'"^') +00 and 
UOc{S S^^^) -00) is given by Eq. (7.70) and ujnc{T) 
is its finitc-T extension. Hence, for this intermediate 
range of 5, we predict on general grounds that this SFp^- 
SFp^_i_ip^ phase boundary Tci(a;o) has a maximum. Thus 
at fixed T the gas for this range of parameters should 



exhibit a reentrant SFr,, 



SFri^ transition 



with detuning u)q. 

Finally, for large dipolar splitting 5f~^^ < (5, p^- 
superfiuid is stable below Tc2(wo) throughout, as illus- 



trated in Fig. 11. We note, however, that for a Feshbach 
resonance splitting much larger than the Fermi energy, 
we expect that on sufficiently short time scales (set by 
time scale for energy relaxation in the system) the two 
(m = and m = ±1) split Feshbach resonances will act 
independently, so that one can come in resonance with 
each of them separately. If so, either Fig. 11 or 9 will 
be experimentally observed, depending on to which of 
the two resonances, m = Oorm = ±l, respectively the 
system has been tuned. 

All of the above discussed transitions are quite conven- 
tional and should be experimentally identifiable through 
(among other signatures) their standard universal ther- 
modynamic singularities (e.g., in heat capacity). Fur- 
thermore, the phases should be distinguishable through 
their spectroscopic properties, with the normal state gap- 
less throughout, and for positive chemical potential, the 
Px-superfluid displaying a gap with an equatorial line of 



nodes, i.e., for kF 



ode 



kp, and p^ 



ipy-superfluid ex- 



hibiting a gap with nodes at the north and south pole. 

In addition to above transitions that are characteriz- 

ablc by an order parameter, wc expect the p-wavc su- 
perfluid to exhibit a number of non-Landau type of (the 
so-called) topological transitions at /x = 0. The sim- 
plest argument for the existence of such transitions is 
the fact that a p-wave supcrfluid exhibits the aforemen- 
tioned gapless excitations around Fermi surface for /z > 0, 
and is gapped to single-particle excitations for ^ < 0, as 
clearly seen in i^k, Eq. (7.6). Thus we predict contin- 
uous ^ ^FSapped ^nd SF^^l^'' SF^'^IZ'^ 

Px Px Px~rtPy Px~rtPy 

transitions at /i = 0. One might expect a clear signature 
of such transitions from the change in the low-T thermo- 
dynamic behavior, e.g., with the heat capacity changing 
from a power-law in T to an activated form, with gap 
for < set by the molecular binding energy. How- 
ever, a collective sound mode, present in any superfluid, 
that contributes power-law in T contributions, might ob- 
scure the distinction between the /it > and < p- 
wave superfluid phases. Local spectroscopic experimen- 
tal probes (some atomic gas analog of tunneling exper- 
iments) should prove useful for detection of these tran- 
sitions. Despite a lack of local Landau order parame- 
ter, these weakly- and strongly-paired p- wave superfluids 
are distinguishable by their topological properties [10- 
12, 46, 54], as we discuss in the next section. 



VIII. TOPOLOGICAL PHASE TRANSITIONS 
AND NON-ABELIAN STATISTICS 

In addition to the rich (but conventional) phenomenol- 
ogy of p-wave resonant gases obtained in previous sub- 
sections, as we will show next, they can also exhibit a 
more subtle (in some cases topological) order and associ- 
ated phase transitions, that cannot be classified by a lo- 
cal order parameter, nor associated Landau theory [68]. 
The existence of such continuous non-Landau type phase 
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transitions have long been appreciated in the literature. 
Examples range from Anderson's metal-insulator tran- 
sition and transitions between different quantum Hall 
ground states to the well-understood Kosterlitz-Thouless 
transition (e.g., in superfluid films) and two-dimensional 
melting, all separating two disordered states undistin- 
guished by any local order parameter or conventional 
symmetry operation. 



A. P-wave superfluid in 3 dimensions 

It is remarkable that p^- and Px + ip^-superfluids 
are examples of a system, that can undergo such non- 
Landau type phase transition when the chemical poten- 
tial changes sign. In three dimensions this can be simply 
seen from the qualitative change in the spectrum E^k, 
Eq. (7.6) of single-particle fermionic excitations, that for 
the Px and + ipy states 

{Bx,By,B,)p^ = A (1,0,0), (8.1) 

{Bx,By,Bz)p^+ipy = (1,^,0), (8.2) 



are respectively given by 




<Px+ip 



— -m) +\mkl+kl), (8.4) 



^'Px~\~'^'Py Q 

ky = 0, and 



with A the maximum gap of each state, related to the cor- 
responding order parameter. Clearly, in a Px-superfluid, 
for n > 0, E^" = (i.e., gapless) for = and ky, kz 
arbitrary, and for < 0, > (i.e., gapped) for all k. 

Similarly, in a,px+ipy-snpei&uid, for > 0, E^ 
(i.e., gapless) for kz = \/2mfj. and kx 

for /X < 0, > (i.e., gapped) for all k. Physi- 

cally these spectral distinctions arise because for > 0, 
a phase that we refer to as SF^'^"-'', the pairing is a col- 
lective Fermi surface phenomenon and finite angular mo- 
mentum forces the gap to vanish on some subspace of 
the Fermi surface. On the other hand, for /x < 0, in the 
gpatrong gg^p jg siuglc-pair of fcrmious phenomenon 
and is simply set by the molecular binding energy, inde- 
pendent of the angular momentum state of the molecule. 

These changes in the spectrum lead to qualitatively 
distinct single-particle correlation functions and there- 
fore require a genuine quantum phase transition (illus- 
trated in Figs. 9, 10, 11) separating two distinct (weakly- 
and strongly-paired) Px- and two distinct (weakly- and 
strongly-paired) px + ipy-superfluids, as /i changes from 
positive to negative, respectively. 

In a classic BCS p-wave paired superfluid, such as He^, 
fj, is always positive and these transitions are not experi- 
mentally accessible. However, in p-wave resonant atomic 



gases they should be easily realizable (if a p-wave super- 
fluid is produced) by changing the detuning parameter 
Wo (controlled by an external magnetic field), that is, as 
we have shown in Sec. VII, is closely tracked (up to terms 
of the order of gp) by /x. 

In addition to the above quasiparticle spectrum and 
correlation function argument for the transition at /x = 0, 
the existence of the SF^^^^^^^ ^^pT+ipy transition can 
be also seen by noting that the two types of px + ipy- 
superfluids can be distinguished by topological order as 
discussed in detail in Refs. [12, 46, 54]. Although, as 
argued above a spectral distinction between /i > and 
A* < Px-superfluids exists, and therefore we expect a 
corresponding SF^J''''' SF^^J'"^^ transition at = 0, 
we are not aware of any topological distinction between 
these two phases similar to the Px+iPy classiflcation. For 
the rest of the section, below, we will focus on the anal- 
ogous transition in two dimensions, where for a Pa, -|- ipy 
superfluid clearly no spectral distinction exists, with both 
SF;"'^^ and S'Ff ™"f gapped, but stiU distinguishable 



by topological order. 



B. P-wave superfluid in 2 dimensions 

The three-dimensional calculations of this paper, argu- 
ing for the existence of a p-wave superfluid can be easily 
extended to two dimensions, with only minor quantita- 
tive distinctions (e.g., the dimensionless parameter C2 in 
2D scales logarithmically with the uv cutoff A). Thus, we 
expect the existence of a fully-gapped two-dimensional 
Px + ipj,-superfluid for > and /i < 0. As we will 
see, a plethora of especially interesting phenomena takes 
place in such a system, that we expect to be realizable 
by confining the degenerate atomic gas to a highly oblate 
magnetic trap. 

Although much of our discussion of this system follows 
an excellent paper by N. Read and D. Green, Ref. [11], 
as well as Refs. [12, 54, 55], we elaborate on details of the 
analysis (particularly on the existence of the zero modes), 
and thereby hope to elucidate a number of points dis- 
cussed there. Furthermore, while above papers are well- 
known and appreciated in the quantum Hall community, 
they are less familiar to the atomic community and thus, 
their main results are worth elaborating on here. 

Following Ref. [11], let us first construct a ground 
state wave function of a two-dimensional px + ipy super- 
fiuid. As discussed in Sec. VB its mean-field Hamilto- 
nian (valid in the narrow-resonance limit) follows directly 
from Eq. (5.34), with the substitution bp^a SpflBa and 
Eq. (8.2), and is given by 



H - jjL Nf = ^ ^fc a]^ak | A (fc^ + iky ) aj^a 

k k 

-I- A (fcj, - iky) S-kSk 



At 



(8.5) 



Here we fix the phase of A by choosing A = A and, as 
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before, is given by Eq. (6.25), 

This Hamiltonian is diagonalized by a unitary transfor- 
mation to the Bogohubov quasiparticles 



7k = ulcLi, + v^aiy., 
7k = t^ka-k + WkSj,, 



taking the form 



^ = £GS + X]^k7k7k. 



(8.6) 



(8.7) 



The ground state of this Hamiltonian is similar to its 
s-wave counterpart Eq. (6.28), and is given by 



(8.8) 



with each pair k, — k in the product taken only once. 

Here ?ik, Wk arc p-wavc analogs of Eq. (6.25), and satisfy 
the Bogoliubov-dc-Gcnnes equations 



-A (fcx - iky) 



-A {kx + iky) \ f u-k \ ^ I Uk 



(8.9) 

The solution of these equations is straightforward with 
the result Eq. (8.4) and with normalized Uk, Wk being 



Uk 



2Ek ' 

iky) , 



V2Sk(£^k + a) 



(8.10) 



We note that unlike the s-wave case Eq. (6.29), the rel- 
ative phase of Uk and Wk is nonzero. Let us construct a 
real space version of Eq. (8.8). It is given by 

*(ri,r2,...) = ^(-l)^fif(rpj - rpj g{rp^ - rpj . . . . 

(8.11) 

Here rj are two-dimensional vectors denoting the posi- 
tion of the i-th fermion, and g{r) is a Cooper-pair (or 
molecular) wavefunction given by 



5(r) 



J (2^)2 



Vk_ 
Uk 



.12) 



P stands for a permutation of numbers 1,2, Nj, 
where Nf is the total number of fermions, and (—1)^ 
is the sign of the permutation, thereby enforcing the an- 
tisymmetrization of the many-atom ground-state wave- 
function. Notice that g{r) = —g{—r) due to the p-wave 
symmetry of the superfluid (since Uk — u_k, i^k ~ ~w-k). 

Now suppose /u > 0. Then at small k <C \/2miJ, or 
equivalently |r| n~^/^, we can estimate the function 
to be integrated in Eq. (8.12) to go as 



Uk 



Ek + ^k' 



kx iky 

fc2 



(8.13) 



It immediately follows that for |r| ::g> Z, 



5(r) 



.14) 



where z is the complex number representing the two- 
dimensional vector r = xx + yy as z = x + iy. Therefore, 
for > the wave function takes the form 



^{zuz2,...)=j2(-^y 



Zp-^ — Zp^ Zp^ — Zp^ 

(8.15) 

This wave function occurs in the context of the quantum 
Hall effect (modulo the Gaussian and Jastrow factors not 
essential for the present discussion) and is called the Pfaf- 
fian or Moore- Read state [69]. 

To understand the connection with the quantum Hall 
effect, we recall that for the last few years attempts have 
been made to realize quantum Hall states [70, 71] out 
of Bose- Einstein condensates by rotating them [72]. In 
the px + ipy condensate, thanks to the relative angular 
momentum £ = I oi each Cooper pair (or closed-channel 
molecule), the (fermionic) condensate already automat- 
ically rotates and therefore does not require any exter- 
nally imposed rotation to be in the quantum Hall ground 
state. 

A key observation is that for /i < Uk/fk ^ kx — iky 
at small k and g{r) no longer has the power- law fall off 
characteristic of the quantum Hall-like ground state in 
Eq. (8.15). Instead, the integral in Eq. (8.12) is then 
dominated by large k, and generally we expect that ^(r) 
will be an exponentially decaying function. The authors 
of Ref. [11] referred to the /i < as the strongly-coupled 
phase. For BCS-BEC condensates studied here, /z < 
corresponds to the BEC regime obtained for negative de- 
tuning Wo- 

As mentioned above, despite the qualitative distinction 
in the ground-state wavcfunctions, in a two dimensions 
Ek > 0, i.e., gapped for all k, for both /x > and fj, < 0. 
The only special point where there are gapless excitations 
is /X = 0. Nevertheless given the qualitative distinction 
between the ground states at /i > and n < 0, they 
much be separated by a quantum phase transition at the 
gapless point /i = 0. The situation is again reminiscent of 
quantum Hall transitions where gapless points separates 
gapped quantum Hall states. 

Although this transition at /x = is not of Landau 
type (not exhibiting any obvious local order parameter) 
the weakly- (/i > 0) and strongly-paired (/i < 0) Px + ipy 
states are topologically distinct and therefore the transi- 
tion is topological. The topological distinction lies in the 
properties of the two complex functions Mk and Wk, con- 
strained by |ukp + |uk|^ = 1- Since their overall phase 
is unimportant, they arc parametrized by two real pa- 
rameters. Thus, ?ik and t'k represent a map from a two- 
dimensional space of k (which can be thought of topolog- 
ically as a sphere 52, if the point A; = oo is added) to the 
two-dimensional space of Uk, Wk- Such maps arc char- 
acterized by the winding numbers, called the homotopy 
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classes corresponding to the homotopy group 7r2 (52) = Z. 
Roughly speaking, these winding numbers are the num- 
ber of times Wk, fk wraps around a sphere as k varies. 
Quite remarkably, one can see that these numbers are 
different for ^ > and /i < 0. 

To see; this explicitly we construct a unit vector n which 
points in the direction of the spinor (uk, i^k)- To this end, 
recall the standard relation between a spinor tpa and a 
vector n^, = a-^^'^*V'/3, which gives 



krrA 



ik'^k + -"k^k 



II. 



i^k 



(8.16) 



It is important to keep in mind that A is in fact a function 
of fc^, being a constant for fc <C A, but quickly dropping 
off to zero at fc A, where A is the ultraviolet cutoff as- 
sociated with the interatomic potential range. The wind- 
ing number associated with 7r2(52) is given by the well 
known topological invariant (discussed in our context in 
Ref. [12], whose notations we borrow here) 



N:- 



Stt/ 



(8.17) 



Substituting Eq. (8.16) into the expression for the 
topological invariant we fin d, af ter an appropriate rescal- 
ing of k and with A = A/\/2to 



1 r 

2 70 



kdk A 



((fc2 - + A;2A2) ' 

(8.18) 

As required by the general form of A^s, Eq. (8.17), this 

expression is a total derivative, and the integral can be 
computed directly with the result 



k—oo 



sign n) . 



fc=0 



(8.19) 

Thus, for /X < 0, A^3 = 0, while for /z > A^3 = 1, and 
indeed Uk, ^^k define a topologically nontrivial map only 
for /X > 0. Hence, a Px + ipj,-superfluid ground state 
exhibits topological order only for /i > 0. 

It is interesting to observe that the topological invari- 
ant N3 for the -state gives N3 = independent of /x, 
since its Uk, fk are real and therefore define a trivial 
map. The same is true for an s-wave condensate. Thus, 
at least based on this topological invariant, neither of 
these states are topological, nor is the transition between 
weakly- (BCS) and strongly-paired (BEC) states in these 
systems. 

Finally, we remark that the topological invariant N3 

constructed here constitutes a particular case of more 
general topological invariants studied in Ref. [12]. 



C. Vortices and zero modes of a two-dimensional 
Px + ipy superfluid: non-Abelian statistics and "index 
theorem" 

We can further elucidate the nature of the Px + ipy 
condensates if wo study the solutions to the Bogoliubov- 
de-Gennes (BdG) equation in the presence of vortices in 
the condensate A. In fact as we will see below, a non- 
trivial topological order exhibited by the wcakly-pairod 
(/Lt > 0) Px + «Py-superfiuid will reflect itself in the nature 
of the spectrum in the presence of vortices. Recall that 
a phase of the condensate wavcfunction changes by an 
integer number times 27r every time one goes around the 
vortex. Thus, in the presence of collection of vortices at 
positions z^, the gap function A(r), proportional to the 
condensate wavefunction can generally be written as 



./2 



D{v), 



(8.20) 



where D{r) is a fmiction of position whose phase is single 
valued. Since its square is the condensate density, and 
D{r) is expected to vanish inside vortex cores. 

Generically. in the presence of vortices, one expects 
solutions localized on them. It has been appreciated for 
some time, based on a variety of arguments [11, 12] (with- 
out an explicit solution of the Bogoliubov-de-Gennes 
equation) , that a,px+ ipy-superfluid is special in that its 
fundamental 27r vortex in thermodynamic limit is guar- 
anteed to carry a state (referred to as "zero mode") at 
exactly zero energy. 

Recently, wc have studied the question of existence and 
robustness of such zero modes for the more general prob- 
lem of a collection of vortices [55] . As we will show below, 
we found that for a macroscopic sample (i.e., ignoring the 
boundary physics), without fine-tuning, strictly speaking 
there is only one or zero Majorana-fermion mode depend- 
ing only on whether the total vorticity of the order pa- 
rameter (in elementary vortex units of 27r) is odd or even, 
respectively. For a collection of well-separated vortices, 
within an exponential accuracy one zero mode per an iso- 
lated odd-vorticity vortex persists. As two of such vor- 
tices arc brought closer together the corresponding pair of 
"zero" modes splits away to finite ±E (vortex-separation 
dependent) energies. Generically, even-vorticity vortices 
do not carry any zero modes. 

Before wc proceed to construct these solutions explic- 
itly, let us discuss in general what we expect from the 
solutions of these BdG equations. A generic Bogoliubov- 
de-Gennes Hamiltonian can always be represented in the 
form 



H = \ a\hijaj - ajhija] + ajAij-a] + a]A*^-aJ 



(8.21) 

Here the indices i, j represent a way to enumerate 

fcrmion creation and annihilation operators, being for ex- 
ample, points in space and/or spin, if the fermions also 
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carry spin, hij is a hermitian operator, while Ay is an 
antisymmetric operator. The study of this Hamiltonian 
is then equivalent to the study of a matrix 



h 

At 



A 



5.22) 



This matrix possesses the following important symmetry 
property 



aiH(j\ = —Ti.*. 



(8.23) 



Here ai is the first Pauli matrix acting in the 2 by 2 space 
of the matrix Eq. (8.22). In the terminology of Ref. [73], 
wc say that this matrix belongs to symmetry class D. 
As a result of this property, if tp is an eigenvector of this 
matrix with the eigenvalue E, then aiip* has to be an 
eigenvector with the eigenvalue —E. Indeed, 



Haii)* = -ain*tl>* = -Eait/j*. 



(8.24) 



As a result, all nonzero eigenvalues of H come in pairs, 

ztE. A special role is played by the zero eigenvectors of 
this matrix, namely the zero modes discussed above. If 
■0 is a zero mode, aiip* is also a zero mode. Taking linear 
combinations tp + ciiip*, i{tp — critjj*) of these modes, we 
can always ensure the relation 



.25) 



for every zero mode. In the absence of other symme- 
tries of H it is quite clear that generically there is noth- 
ing that protects the total number of its zero modes 
under smooth changes of the Hamiltonian matrix that 
preserve! its BdG form, namely retain the properties in 
Eqs. (8.22) and (8.23) . However, since non-zero modes 
have to always appear and disappear in ±E pairs, as 
long as the symmetry property (8.23) is preserved by the 
perturbation, the number of zero modes can only change 
by multiples of 2. Thus, while the number of zero 
modes of the Hamiltonian (8.22) may change, this num- 
ber will always remain either odd or even, with (—1)^^ 
a "topological invariant" [74, 75]. 

The value of this invariant is easy to establish if one 
observes that Ti is an even-sized matrix, with an even 
number of eigenvalues. Since the number of non-zero 
modes must be even, this implies that the number of zero 
modes is also even. Thus (—1)^^ = 0, and generally the 
BdG problem does not have any topologically protected 
zero modes. Furthermore, since, as demonstrated above, 
zero modes must appear in pairs, there can only be an 
even number of accidental zero modes, which will nev- 
ertheless be generally destroyed by any perturbation of 
H (preserving its BdG structure Eq. (8.22)). We believe 
this observation was first made by N. Read [75]. 

The situation should be contrasted with that of the 
Dirac operators V. Those operators, being generally of 
one of chiral classes in the terminology of Ref. [73], obey 
the symmetry 



Thus if Tp is an eigenvector of P with the eigenvalue E, 
a^ip is an eigenvector with the eigenvalue —E. The zero 
modes of T) must obey the relation 

(^3'4'L,R = ±'tpL,R- 

Here "left" zero modes tpL come with the eigenvalue 
while "right" zero modes ipn have the eigenvalue —1 of 
the operator a^. As the operator V is deformed, the 
number of zero modes changes, but the non-zero modes 
always appear in pairs where one of the pair has to be 
"left" and the other "right" . Therefore, while the num- 
ber of zero modes is not an invariant, the difference be- 
tween the number of left and right zero modes has to 
be a topological invariant, determined (through the in- 
dex theorem) by the monopole charge of the background 
gauge- field. 

Contrast this with zero modes of H, which obey the re- 
lation Eq. (8.25). Because of complex conjugation of -0, 
these zero modes cannot be split into "left" and "right" . 
Indeed, even if we tried to impose aiip* = —ip, a sim- 
ple redefinition of tjj iijj brings this relation back to 
Eq. (8.25). Thus, the most an "index theorem" could 
demonstrate in case of the Bogoliubov-de-Gennes prob- 
lem, is whether there is or exactly 1 zero mode. More- 
over, since the Bogoliubov-de-Gennes problem is defined 
by an even dimensional Hamiltonian, generically there 
will not be any topologically protected zero modes [75] . 

Yet it is quite remarkable that in case of an isolated 
vortex of odd vorticity in a macroscopic sample (i.e., ig- 
noring the sample boundaries) of apx + ip^-superfluid of 
spinless fermions, there is exactly one zero mode local- 
ized on this vortex [11, 56, 76, 77]. To be consistent with 
above general property of the BdG Hamiltonian (namely, 
that the total number of zero modes must be even) an- 
other vortex is situated at the boundary of the system 
[11, 75], preserving the overall parity of the number of 
zero modes. Hence, although even in this odd-vorticity 
case the one zero mode is not protected topologically, able 
to hybridize with a vortex at a boundary of the sample, 
it survives (up to exponentially small corrections) only 
by virtue of being far away from the boundary (and from 
other odd-vorticity vortices). 

To sec this explicitly we now consider a + ipy- 
superfluid in the presence of a single rotationally sym- 
metry vortex, characterized by 



A(r) 



Hip f2 



nr), 



(8.26) 



-V. 



where /(r) is a real function of r (vanishing at small 
r), / is the vorticity of the vortex, r, ip are the polar 
coordinates centered on the vortex, and the factor of i/2 
is chosen to simplify subsequent calculations. In this case 
the Bogoliubov-de-Gennes equations take the form 



-^-M^ w(r)-/(r)e"2^^ e'^ f{r)v{v) = Eu{r), 
^-h/x')w(r)-/(r)e-'^^ e'"^ f{r)u{v) = Ev{r), 
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(8.27) 



We remark that once solutions to these equations ?in(r), 
Vn{r), corresponding to energies En, are known, the Bo- 
goliubov quasiparticle creation and annihilation opera- 
tors are given by 



7« 



(fr [<(r)a(r) +<(r)a^(r)] 
t= j dPr[un{v)d\v)+Vn{v)a{v)\. (8.28) 



If the condensate was uniform, then the solutions to 

the Bogoliubov-dc-Gcnncs equations would be plane 
waves, immediately leading to Eq. (8.6). The inverse 
to Eq. (8.28) reads 



=^lnUn{Y) +7l,<(r), 

n 



small deformation of, say, the shape of the order param- 
eter shape /(r) will destroy these modes. We note that 
this argument easily accommodates vortices that are not 
symmetric, as those can be smoothly deformed into sym- 
metric ones without changing the topologically protected 
parity of ■ The conclusion is that generically there are 
no zero modes in the presence of an isolated vortex of 
even vorticity. 

The situation is dractically different if the vorticity of 
the vortex is odd, i.e., ii I = 2n — 1. In this case the 
transformation Eq. (8.30) cannot entirely eliminate the 
vortex from the equations (even with the help of a smooth 
deformation), leaving at least one fundamental unit of 
vorticity. This thereby guarantees at least one one zero 
mode localized on the odd-vorticity vortex. To see this, 
recall that due to the condition Eq. (8.25), the zero mode 
satisfies 



u = v 



(8.29) 



(8.32) 



Next we observe that for the case of a vortex of even 
vorticity, I = 2n, we can eliminate the phase dependence 
of Eq. (8.27) entirely. Indeed, making a transformation 



Combining this with the transformation Eq. (8.30), we 
find the equation for the zero mode 



-.f(r)e- 
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leads to equations 
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- IJ, u- 
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in du s 3 . „ , . , „ 

mr-' dip ^ ' dz 



- f{r)u* 
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mr"^ dip' 



(8.33) 



We look for the solution to this equation in terms of a 
spherically symmetric real function u(r). This gives 



2mr2 



in dv 
mr"^ dip 
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(8.31)A- transformation 
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2mr'^ ) ^ ^ 
(8.34) 



Now we note that these equations are topologically equiv- 
alent to the BdG equation without any vortices. Indeed, 
the only difference between these equations and those 
for a uniform condensate is the presence of the terms 
2in/r^[d/dip], r? jr^, and /(r) that is a constant at large 
r and vanishes in the core of the vortex for v <C Tcore* 

We 

can imagine smoothly deforming these equations to get 
rid of the first two terms (for example, by replacing them 
with a {r? jr^ — 2in/r'^[d/dip\) u and taking a from 1 to 
0), and smoothly deforming /(r) into a constant equal to 
its asymptotic value at large r; in order to be smooth, the 
deformation must preserve the BdG structure Eq. (8.22) 
and the vorticity of the condensate, if there is any. These 
equations then become equivalent to Eq. (8.9) for a con- 
stant, vortex-free order parameter with an exact spec- 
trum Eq. (8.4), that for /x 7^ in two dimensional space 
clearly does not exhibit any zero modes. 

As Eqs. (8.31) are smoothly deformed to get rid of the 
vortex, in principle it is possible that its solutions will 
change and that it will develop zero modes (although, as 
demonstrated above, this can only happen in ±E pairs, 
leading to an even number of these). However, these 
modes will not be topologically protected, and even a 



u{r) = x{r) exp - 



m 



f dr' fir'] 
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(8.35) 



brings this equation to a more familiar form 



X 



X 



2m 2mr 



This is a Schrodinger equation for a particle of mass m 
which moves with angular momentum n in the potential 
m/^(r)/8, that is everywhere positive. We observe that 
this potential vanishes at the origin, and quickly reaches 
its asymptotic bulk value m/o/8 at large r. Then for 
H > m/o/8, there always exist a solution to this equation 
finite at the origin and at infinity. Moreover, if /x < 
™/o/8, then the solution finite at the origin will diverge 
at infinity as 



X 



-2mfj, 



(8.37) 



Combining this with Eq. (8.35), we see that w(r) will still 
be a bounded function at infinity as long as /i > 0. Thus 
the conclusion is, there exist zero mode as long as fi > 0. 
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For a special case of the n = vortex of vorticity —1, the 
small and large r asymptotics of the solution we found 
here was discussed recently in Ref. [56] . 

In the simplest London approximation of a spatially 
uniform condensate when /(r) = /o for all r except inside 
an infinitesimal small core, the zero mode localized on an 
isolated odd-vorticity vortex is simply given by 



u{r) = 



V 



2fj,m ■ 



m?^ 1 



}for, for n > mfy 



m 



f 4 

2 Jo 



2m/x I 



r-^o'', for < /i < m- 
(8.38) 



J 



where Jn{x), In{x) are Bessel and modified Bessel func- 
tions. 

We note that it may seem possible to construct ad- 
ditional zero modes in the following way. Instead of 
the ansatz of a rotationally invariant u{r) just after 
Eq. (8.33), we could have chosen an ansatz 



u{r) = «„(r)e'"'^ + u_„(r)e-'«'^. 



(8.39) 



Then two second order differential equations follow, re- 
lating these two functions. These are 
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(8.40) 
(8.41) 



Generally there are going to be four solutions to these 
equations which go as r'"""' or r'""'""! at small r. 
The other two will diverge as r"'"""' or as 7--l"+"l. 
At infinity the four solutions of these equations go as 



exp 



r ^±y/o ± -\/ /o "^^ - 8/um^ . Obviously, only two 

of these solutions are finite at infinity. However, barring 
a coincidence, none of those solutions finite at r = 
are also finite at infinity. Even if they arc for some 
special value of /x, by the above arguments the addi- 
tional zero modes must appear in topologically unpro- 
tected pairs, that will be split to finite ±E energies by 
a slight generic deformation of the potential (order pa- 
rameter distortion). Hence we conclude that generically 
there will be no additional zero modes (except the one 
found above) for an odd-vorticity vortex. 

Thus we indeed find that the number of zero modes 

in a symmetric odd-vorticity vortex must be one. Since 
a smooth deformations of the order parameter can only 
change the zero mode number by multiples of two, an 
arbitrarily shaped odd-vorticity vortex must also have an 
odd number of zero modes . However, any number of zero 
modes other than one is not generic and will revert to one 
under an arbitrary deformation of the order parameter. 

Now for a collection of well-separated r :» 1/(toA) 
vortices of odd vorticity, each of them will have one zero 
mode localized on it. However, as they are brought closer 
to each other, these zero modes will actually split into a 
band of low lying ±E modes [11]. However, since other 
(^xcitcd modes are separated from the zero modes by a 
gap [76] , the narrow band will only mix very weakly with 
other states of the system. 

It is this band of nearly degenerate zero modes that 
exhibit non-Abelian statistics. Following Ref. [78], we 



briefly describe how it is realized here. Each of the zero 

modes is in fact a Majorana fcrmion, as follows directly 
from Eq. (8.28) and condition (8.32), given by 

7 = j (fr (M*(r)a(r) + u(r)at(r)) . (8.42) 

It is straightforward to check that 



f = 1, 



(8.43) 
(8.44) 



when u{v) is properly normalized. A Majorana fermion 
is essentially half of a real fermion, thus they must always 
come in pairs (in case of an odd number of odd vorticity 
vortices, the last remaining Majorana fermion is located 
at the boundary). Given 2n Majorana fermions, we can 
construct creation and annihilation operators of n real 
fermions, according to 



Cj =72j-l +n2j, 



4 = 72j-i -n2j. 



(8.45) 



with 7j a Majorana annihilation operator of a fcrmion 
localized on a vortex at position Yj. Thus clearly a real 
fermion Cj is actually split between two vortices at T2j-i 
and Y2j. 

In the presence of 2n vortices, there are 2" states cor- 
responding to n pair of vortices being either occupied or 
empty. Now it is possible to show that if two vortices 
are adiabatically exchanged - moved around each other 
- these nearly degenerate zero states mix with each other. 
More precisely, 



(8.46) 
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where t/ is a 2" by 2" unitary matrix, representing the 
unitary transformation of the 2" groimd states ipj. The 
matrix U depends on which two vortices are exchanged 
(and on the direction of the exchange). It does not, how- 
ever, depend on the path along which the vortices are 
moved, and is thus topological. 

The matrix U is not a general unitary matrix. In fact, 
all the 2" states should be split into two subsets of size 
2"/^, one with even, and the other with odd number of 
fermions. U only mixes states within each of these sub- 
sets. The reason for this is that the fermions which oc- 
cupy or vacate the zero modes must come in pairs, since 
they are produced from a Cooper pair which is being 
split into two fermions, or being assembled back from 
two fermions. 

Matrices U can be constructed by considering the 
change in ipj as one vortex is slowly moved around an- 
other, while others are kept fixed. For such adiabatic 
change, the effect on the Bogoliubov-de-Gennes equation 
written in the vicinity of the first vortex is simply through 
A slowly changing its phase. This can be incorporated 
into a change of the phase of u and v by absorbing half 
of the phase into u and the other half (with the oppo- 
site sign) into v. As a result, when one vortex moves 
all the way around another vortex, each of its Majorana 
fermions changes sign. A change in sign of the Majo- 
rana fermions can be translated into the change in the 
states tpj, by constructing an appropriate operator such 
that WjjU — > — 7j. Then the action of U on states ipj 
constitutes a transformation as in Eq. (8.46). For a more 
detailed discussion and an expHcit construction of U for 
a, Px + ipy-supcrfluid we refer the reader to Rcf . [78] . 

The transformation (by U) upon exchange of two vor- 
tices in a + ipjz-superfluid is a generalization of a stan- 
dard quantum statistics of bosons and fermions familiar 
from standard quantum mechanics. This exchange trans- 
formaton also generalizes the two-dimensional anyonic 
quantum statistics (familiar from Abelian quantum Hall 
states), where, upon a two-particle exchange a many- 
particle wavefunction gets multiplied by a phase factor 
e'^ (with a phase 9 not necessarily just or tt). Since 
generically unitary matrices U, corresponding to differ- 
ent pairs of vortex exchanges do not commute, the result- 
ing quantum statistics is termed non-Abclian[69]. Thus 
odd-vorticity vortices in a Px + ipj,-superfluid at positive 
detuning (/x > 0) are excitations ("particles") with non- 
Abelian statistics. 

Now, in addition to a basic interest, recent excitement 
about states that exhibit such non- Abelian statistics is 
the observation that they can form a basis for build- 
ing a fault-tolerant "topological" quantum computer [47] . 
More conventional quantum bit (q-bit) schemes, such as 
the Josephson-junction charge, flux and phase q-bits, ions 
in an electrostatic trap, or spin q-bits suffer from de- 
coherence due to interaction with the environment. In 
contrast, a q-bit based on non-Abelian statics, as e.g., 
a state of 2n vortices is topologically protected because 
to change it requires a global operation on vortices such 



as one encircling another, something that environmental 
noise will not generically do. 

Based on the analysis presented here we propose [42] 
that apx + ipy-superfluid, that is likely to be realized in 
a resonant Fermi gas interacting via a p-wave Feshbach 
resonance is a viable candidate for an implementation 
of such a non-Abelian q-bit and associated topological 
quantum computation. One advantage of the realization 
of such a q-bit in degenerate atomic systems (as opposed 
to solid state superconductors) is the tunability of their 
interaction via an external magnetic field, that allows 
a tuning of the chemical potential closer to the fi = 
transition, while taking care to remain in the topological 
phase /i > 0. This in turn will allow a more energet- 
ically stable BEC superfluid, whose transition temper- 
ature and the size of the gap are set by the Fermi en- 
ergy cf, SiS opposed to a tiny fraction of it as in conven- 
tional superconductors stuck in the exponentially weak 
BCS regime. 

Of course, even if such topological p-wavc superfluid 
state and the associated non-Abelian q-bit are realized in 
atomic resonantly-paired condensates, many challenges 
remain, such as a scheme for addressing the q-bits by 
manipulation of vortices and reading off the state tpj of 
their zero modes [56]. 

IX. COMPARISON WITH EXPERIMENT 
A. S-wave 

An important remaining question which must be ad- 
dressed is whether current experiments are characterized 
by a narrow or a broad resonance. While it is gener- 
ally believed that most s-wave Feshbach resonances real- 
ized in current experiments are wide, we will show below 
that some are indeed narrow in the sense that a relevant 
dimensionless parameter 7^ controlling the quantitative 
validity of our theory is small. 

As discussed throughout the paper, theoretically, an 
absolute characterization of a width of an s-wave reso- 
nance is through the value of a dimensionless parameter 
7s ^ g1m?/^ep^^'^ , Eq. (1.5), that is set by the ratio of the 
resonance energy width to the Fermi energy. Although 
this parameter is never measured directly, it can be re- 
lated to experimentally determined quantities through 
the atomic scattering length a{H) as a function of mag- 
netic field H, that is either measured or calculated (see 
Fig. 2). From that data, magnetic field width can 
be extracted, as the range of the magnetic field change 
between the resonance (where a — > 00) and the point 
where a = 0. Alternatively, one can look at the range 
of H where a deviates significantly (e.g., by a factor 2) 
from the background scattering length ai,g. Both meth- 
ods produce similar definitions of i?^. 

By itself, carries little information about the di- 
mensionless many-body resonance width 7^, which, as 
we recall, is not only a properties of the resonance, but 
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depends also on the particle density. That is, to assess 
how wide the resonance is, its width, controlling strength 
of interactions must be compared to another energy scale, 
which in this case is the typical kinetic energy ep- To es- 
tablish a relation between 7^ and Hyj, we recall Eq. (5.13) 



mroLJo 



where wq is a detuning, that measures the deviation of 
the Zceman splitting (between the open and closed chan- 
nel of Feshbach resonance) from its value at the reso- 
nance, where o — > 00 (see Fig. 1). By matching our 
results for the scattering length with its experimental 
dependence on the magnetic field, we determined that 
oJo « {H — Ho), with Hq the field at which the res- 
onance is tuned to zero energy and /i^ the Bohr magne- 
ton. This allows us to express the effective-range length, 
tq (entering the expression for the energy width of the 
resonance; see below) purely in terms of experimentally 
measured quantities, namely: 



Fo| - 



mabgiiBHw ' 



Here h was restored to facilitate calculations below. Once 
ro is found, we can compute 7^ by using 



7s 



I 



^^0.8^. (9.2) 



From ro we can also estimate the intrinsic (density inde- 
pendent) energy width of an s-wave resonance, 



To 
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(9.3) 

(9.4) 



We note that in contrast to a naive guess, this energy 
width of the resonance is not simply the Zeeman energy 
(converted with a Bohr magneton) associated with the 
width- field ff^ . From Fq the dimensionless parameter 7^ 
is then found to be 



(9.5) 



7s ~ V — 




Equation (9.2) can now be used as a criterion on 
whether a resonance is narrow or wide (which is of course, 
equivalent to the one discussed in Sec. IB) . We also 
remark that the use of a^g to find Hyj is completely ar- 
bitrary; we could have instead define Hyj as a range of 
the magnetic field where a{H) exceeds some given value 
|a| > as a reference point, but this would still lead to 



an exactly the same |ro|, Eq. (9.1) (with a^g replaced by 

aw) 



\ro\ ^ 



ma^jfiBHy 



(9.6) 



and the same criterion for the narrowness of the reso- 
nance; only the product a^H^j enters tq. 

Physically, 7s can be though of as the ratio of the (en- 
ergy) range of loq, where a exceeds inter- particle spacing 
I, to the Fermi energy of the gas. Now, of course, it is in 
principle possible to make a resonance narrow by increas- 
ing the atom density n (reducing spacing I). However, it 
is our understanding that due to experimental limita- 
tions, the Fermi energy is typically in the IfiK or less 
range and cannot be significantly increased above this 
value. 

We now apply this dimensionless criterion to the exper- 
iment reported in Ref. [14]. The s-wave resonance stud- 
ied there is in ^Li at Hq « 543.25 G and is probably the 
most narrow one discussed in the literature. These au- 
thors report the density of their condensate to be 3 x 10^^ 
(9.1) cm""^, corresponding to the inter-particle separation of 
I Ri 7 X 10^^ cm w 1.3 • 10'* au. We also note that the size 
of a closed channel molecule is set by the range of the van 
der Waals interaction, which is about 50 au. This length 
d plays the role of the inverse uv cutoff 1/ A of our theory. 
We note that the ratio d/l « 1/250 thus justifying our 
assumption throughout the paper that A^/ {2m) (A) can 
be treated as the largest energy (momentum) scale of the 
system. 

To estimate 7s for the resonance in Ref. [14], we use 
their Fig. 1 for the scattering length a{H) to extract Hyj 
and a-u;. Although from this figure it is difficult to de- 
duce the range of the magnetic field where a{H) is larger 
than the background length Ubg, we use the arbitrari- 
ness of to pick = 500au. This corresponds to 
\H — Ht)\ < 0.015 G, that can be converted into an en- 
ergy by multiplying by Hb- Using Eq. (9.6) then gives 



and 



7s 



6 • 10* au, 



I 

N 



0.2. 



(9.7) 



(9.8) 



We thus conclude that the 543.25 G s-wavc resonance in 
^Li is in fact quite narrow in the absolute, dimensionless 
sense. It is therefore a good candidate for a quantita- 
tive comparison with our predictions for a narrow s-wave 
BCS-BEC crossover. 

However, we note that, in contrast to above estimate 
of 7s based on Fig. 1, according to Fig. 4 of the same 
Ref. [14] the BCS-BEC crossover occurs over the range 
of magnetic fields of the order of 1 G or so, corresponding 
to the range of the detuning of lOO/xK >> ep ~ 1-5/Lt K. 
Thus based on this Fig. 4, using the narrow-resonance 
theory we would instead conclude that the resonance is 
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wide. The reason for this discrepancy is currently un- 
clear to us. One should also notice that the authors of 
Ref. [14] were unable to convert all the atoms into the 
molecules, so perhaps there were other factors in their 
experiment which made its direct comparison with the 
narrow-resonance theory difficult. 

Although the resonance in ^Li at 543.25 G is unusually 
narrow, more typical resonances have magnetic width 
Hyj ~ 10 G with ~ 50 au [22]. For such a resonance 
(assuming ^Li for the mass m of the atom) , 

[rol^ilO^au, (9.9) 

with 

7, ~ ^ « 10. (9.10) 

Thus a more typical s-wave Feshbach resonance experi- 
ments lie in the class of wide resonances. 



B. P-wave 

Although so far no atomic p-wave BGS-BEG super- 
fluid has been realized, p-wave Feshbach resonances have 
been demonstrated and explored experimentally. To get 
a sense of future p-wave superfluid possibilities, it is use- 
ful to look at the Ref. [37]. Unfortunately, as we will 
sec below from the data reported there it is not possi- 
ble to extract gp. Nevertheless some conclusions can be 
made about which phases of a p-wave condensate may be 
realized with the p-wave Feshbach resonance in "'"K. 

We first look the data concerning the value of param- 
eter fco from Eq. (2.21). The parameter c, as given in 
that paper (Eq. (8) of Ref. [37]) is magnetic field de- 
pendent, but in the relevant range of the magnetic field 
it is roughly fco ~ —0.04 au"'^. As in the previous s- 
wave analysis, we estimate the uv cutoff A to be roughly 
of the inverse size of the closed- channel molecule, i.e., 
A K, 0.02 au^^. This indicates that in the expression for 
fco, Eq. (5.52), 

127r / 2a\ 

most likely the dimensionless uv parameter C2 = 
S'pm^A/(37r^) » 1, which gives 

fco w -1a. (9.11) 

TT 

Otherwise, jfcoj would have been much bigger than A. 
Thus we deduced that C2 1 for the experiment of 
Ref. [37]. 

This implies that this experiment is done in the regime 
where the mean-field theory considered in this paper 
might become quantitatively unreliable [45]. Assuming 
that it does not, large C2 indicates that even in the BEG 
regime of this system, most of the particles will be in the 



form of free atoms, not bosonic molecules, as indicated 
in our analysis above. 

Since the experiment reported in Ref. [37] is likely to 
be in the regime where C2 ^ 1, it is impossible to extract 
the Feshbach resonance coupling gp. We expect that gp 
gets renormalized in the regime of large C2, so that its 
bare value simply drops out. However, since the com- 
plete theory of p-wave superfluids at large C2 is yet to 
be constructed, we cannot tell what this implies for the 
experiment [45]. 

Next, we note that the dipolar- interaction splitting 5 
between the = and Wz = ±1 Feshbach resonances 
is quoted in this paper as approximately 4/LtK. This is 
presumably several times bigger than currently experi- 
mentally achievable ep^ that are typically in the range 
of 0.5 to 1/zK. Thus we conclude that under conditions 
described in Ref. [37], at low temperatures the gas will 
be in the Pa;-superfluid state. However, because the split- 
ting is considerably larger than e_F it might be possible 
to bring = ±1 molecules in resonance with the atoms 
independently of = molecules. If so, a p^ + ipy 
superfluid state should be realizable for tuning near the 
rriz = ±1 resonance. 



X. CONCLUSIONS 

In this paper we presented a study of a degenerate 
Fermi gases interacting through a tunable narrow Fes- 
hbach resonance, as recently demonstrated experimen- 
tally. Starting with an analysis of the two-body scatter- 
ing physics we developed and justified generic models for 
description of such systems. We paid a particular atten- 
tion to regimes of validity for a perturbative analysis of 
such systems at finite density, and showed the existence 
of a small dimensionless parameter, the ratio of the Fes- 
hbach resonance width to the Fermi energy. It allows 
perturbative description throughout the full BEG-BGS 
range of detuning, within the framework of two chan- 
nel model. Focussing on the most interesting cases of 
the s- and p-wave resonances, we analyzed in detail the 
corresponding systems. For the s-wave resonance, we ob- 
tained predictions for the behavior of the system across 
the BEG-BGS crossover, that we expect to be quantita- 
tively accurate for the case of a narrow resonance. For 
the far richer p-wave resonance, dominant for a single hy- 
perfine species of atoms, we showed the existence of and 
analyzed a number of classical, quantum and topological 
phase transitions exhibited by this system as a function 
of temperature, Feshbach resonance detuning and reso- 
nance dipolar splitting, and calculated the correspond- 
ing phase diagrams, illustrated in Figs. 9, 10, and 11. 
Finally, we studied topological properties of the weakly- 
paired px + ipj,-supcrfluid, as well as zero-modes inside 
vortices of such a topologically-ordered superfluid in two 
dimensions. We hope that our analysis will be useful for 
probing the associated non-Abelian quantum statistics of 
such vortices, and more generally, for experimental real- 
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ization and studies of a resonant p-wave superfluidity in 
degenerate atomic gases. 
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APPENDIX A: BOSONIC VACUUM 
PROPAGATOR 

It is instructive to analyze the renormalizcd vacuum 
propagator of the 6-particles of the two-channel model. 
We will do it in the s-wave case. Since the renormal- 
izcd propagator is nothing but the sequence of diagrams 
depicted on Fig. 18 with external legs amputated, its 
calculation parallels that of the two fermion scattering 
amplitude. The answer is given by 



£>(k,w) 



1 



fr2 .„2^3/2 / Z 



(Al) 



4to 



The rcnormalized propagator is of course simply propor- 
tional to the two fermion scattering amplitude. Thus 
the poJes of this propagators, which describe the phys- 
ical bound states (or resonances) of two fcrmions, coin- 
cide with the poles of the scattering amplitude. They are 
given by 



1^2 

W = h Wo + 
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It is now straightforward to calculate the residue of the 
propagator D at its poles. The result is 



Z = 



Stt 



Wo - 
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Although the result is rather cumbersome, it is straight- 
forward to see that the residue goes to 1 if Qs goes to 
zero, and it goes to if gs goes to inflnity. Thus, the con- 
tribution of b to the actual bound state of two fcrmions 
(physical molecule) reduces to zero in the limit of large 
(wide resonance limit). Notice that while the size of the 
closed channel molecules b is of the order oid = 2tt/A, the 
actual size of the molecule, which is a superposition of b 



and a cloud of open channel fcrmions, could be quite large 
(and is in fact of the order of the scattering length a). 
Thus one common perception that since the molecules 
are large b cannot be a point particle is incorrect, b is 
not a molecule, but only its closed channel part, whose 
contribution to the actual molecule may in fact reduce to 
zero in a wide resonance regime. 



APPENDIX B: SCATTERING MATRIX VIA 
SINGLE-BODY HAMILTONIAN 

In this appendix we compute scattering amplitudes 
of a number of models relevant to the problem of res- 
onantly interacting Fermi gases. The results that we find 
here have already been obtained in Sec. Ill, Sec. IV and 
Sec. V of the main text, working directly with many-body 
Hamiltonians. However, a problem of two particles inter- 
acting with a potential U{ri — r2) can always be reduced 
to a decoupled evolution of their center of mass and that 
of their relative coordinate r = ri — r2, whose dynamics is 
governed by a Hamiltonian for a single effective particle 
with a reduced mass rrir = mim2/{mi+m2), moving in a 
single-body potential U{r). Hence, a many-body Hamil- 
tonian, when restricted to act on a two-particle Hilbert 
subspace (as in the computation of two-particle scatter- 
ing amplitude) has an equivalent single-particle Hamil- 
tonian with which scattering physics can be equivalently 
straightforwardly analyzed. Thus the analysis in this ap- 
pendix will complement the main text in that we will 
compute scattering amplitudes for many-body models 
studied there using equivalent relative coordinate single- 
particle Hamiltonian. 



1. Fano- Anderson model 

This is a model of a particle, created by aj^, moving 
freely in space (representing the "open channel"), which 
when it hits the origin can turn into another particle, 
created by b^. The 6-particle (representing the "closed 
channel") cannot move at all and has a fixed energy cq. 
The Hamiltonian of this problem can be written as 

-1 ^ = I]:^4«k + eo6^6 + 5. (6^a(0) + at(0)6) , (Bl) 

k 

where gs is the interconversion rate between a- and 6- 
particles, and rrir is the mass of the o-particle. This is 
called the Fauo- Anderson model [25], or a model of a lo- 
calized state in the continuum [79]. cq is a parameter 
which plays the role of "detuning". It is the energy of 
the fe-particle if it was left alone and did not interact with 
the a-particle. The model represents the two-body ver- 
sion of the s-wave two-channel model given by Eq. (5.6) 
(where to,, corresponds to the reduced mass of fermions 
in Eq. (5.6), hence the notation). 
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The scattering amplitude of the a particles can be 
easily evaluated using the T-matrix formalism. The T- 
matrix is given by 

7k,k' = gsD{E)g, + g,D{E)gsHE)9sD{E)g, + ... 



D-^{E)-glU{Ey 
where D{E) is the Green's function of the 6-particle 

1 



(B2) 



D{E) 



(B3) 



E - eo + iO ' 

and n(£) is the trace of the a-particle Green's function 
d^q 6{\-q) 



(B4) 



and E = fc^/ {2mr). The value of I1{E) was already com- 
puted in Eq. (4.3). Doing the algebra, we arrive at 



where wq is the "renormalized" energy of the 6-particle, 

Wo = Co - ffgrnrA/TT^. (B6) 



Wc sec that the scattering length and effective range ex- 
tracted out of Eq. (B5) 
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(B7) 



coincides with Eqs. (5.13), (5.14). 



2. Hybrid model 

In the literature it is popular to consider Fcshbach- 
resonant interactions together with the interactions via 
a short range potential. The single-particle Hamiltonian 
which captures a combination of these interactions takes 
the form 



(B5) ^ ^Y. ^aluk+eob^b+X a^{0)a{0)+g, {b'<a{0) + a\0)bj 



2mr 

(B8) 

The many-body version would then be a combination 
of a two-channel model with a direct four-fermion point 
interaction scattering term 
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It is instructive to calculate the scattering amplitude of 
the a-particlcs which follow from Eq. (B8). The calcula- 
tion largely parallels that given in Eq. (B2), except ggDgs 
gets replaced hy ggDgg + X. After some algebra we obtain 



fo{k) = 



1 



2A 



(BIO) 



+ ik 



We see that the a-particles scatter in a rather compli- 
cated fashion. If we are only interested in low-energy 
scattering, we can expand the denominator of the scat- 
tering amplitude and read off the scattering length and 
the effective range as 



27ren 



2A 



a = 



27Tg^s 



mr{€o\-gl) TT ' " ml{eo\-glf 

(Bll) 

In principle, wc can now redefine the parameters e^, g, 
and A in such a way that the scattering length and the 



effective range computed here would coincide with the 
ones produced by the pure Fano- Anderson model. We 
see that ultimately including both the 5-like potential 
and the Fano- Anderson term in the Hamiltonian does not 
produce new low energy physics compared to the pure 
Fano- Anderson case, and amounts to just redefining the 
parameters of the Fano- Anderson model - and thus, of 
the two-channel model Eq. (5.6). Its only physical effect 
is to accommodate for Obg ~ d absent in the pure two- 
channel model. 



We also remark that this may sometimes not be true 
in lower dimensions. In ID the inclusion of the contact 

interaction term may change the physics described by the 
two-channel model qualitatively [80]. 
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3. P-wave Fano-Anderson Model 

The p-wave version of the Fano-Anderson model 
Eq. (Bl) is given by 



(bx - tiy^ = 6m=-l- 



(B13) 



k ^ a=l ^ k,a 

(B12) 

Here the a particle scatters in the p-wave channel and 
can convert into a b particle which carries internal an- 
gular momentum i = 1, since the angular momentum 
is conserved. The angular momentum is represented by 
the vector index a. It is related to the states with defi- 
nite projections of angular momentum &m=i, &m=Oj a-nd 
bm=-i via the standard formulae (Ref. [49]) discussed 
in Eqs. (5.35), (5.36), (5.37), which we repeat here for 
convenience 



By construction, a-particles scatter only in the p-wave 
channel. The scattering amplitude^ can be easily calcu- 
lated in the same T-matrix formalism, as in the s-wave 
case, Eqs. (B2), (B5). 

The propagator of the 6-particles is now 



(B14) 



bz = b. 



m=0, 



with D{E) given by Eq. (B3). The T-matrix is given by 
the p-wave version of Eq. (B2), 



Tk,k' = ^ gpkaDk'^gp + ^ gpkaDgpk'^G{k" , E)gpk'^Dgpk'(j + ...= 
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where G{E, k) is the Green's function of the a-particles 
and n(£) is now 



n(^;) = 1 1 
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Here E = k"^ /{2mr). Just like everywhere else through- 
out the paper, we cut off the divergent integral at ~ A. 
Unlike Eq. (3.26), the integral is divergent as and pro- 
duces two cutoff dependent terms. With the help of the 
notations Eq. (5.45), Eq. (5.46), introduced earlier, we 
find the p-wave scattering amplitude 
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This coincides with the result of the many-body calcu- 
lations reported in Eq. (5.48) with the exception of a 
numerical coefficient. This difference is the result of the 
indistinguishability of identical particles which was im- 
portant in Eq. (5.48) but played no role here, in the one- 
body scattering calculation. 



APPENDIX C: DETAILS OF THE p-WAVE 
SADDLE-POINT EQUATION AND FREE 
ENERGY 

The thermodynamics of a p-wave supcrfluid is com- 
pletely determined by the free energy, Eq. (7.5) and the 
corresponding saddle-point equation, Eq. (7.7), derived 
from it. These were expressed in terms of one key tensor 



/^P[B], defined in Eq. (7.8). Here we compute /^y[B 



(T)r 



a/3 



at zero and finite temperatures and thereby obtain the 
corresponding ground-state energy and the free energy. 



1. Zero temperature 



At zero temperature l'^^ [B] reduces to 
(fik kakjj 
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where we used the spectrum i?k, Eq. (7.6). The integral 
is naturally computed in spherical coordinates, with the 
radial part over k conveniently expressed as an integral 
over the free spectrum e/j = k"^ /2m 



Ic,^ = gl 2m\ii\N{\ii\) 
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where we defined a function 
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that arises from an integral over e scaled by the chem- 
ical potential e ^ e/\fi\, N{fi) = m^/^\n\^/y{2^/'^Tr'^) = 
c|/x|^/^ is the density of states, fi = = ±1, Ea = 

(AV2m)/|/i|, and 



Smgf 



= ^((u.kr + (v.k)^). 



(C4) 



(C5) 



Because at large e, I{Q) scales as e^/^, its one set of 
dominant contributions comes from the region of integra- 
tion near the uv cutoff E\. We isolate these I\ contri- 
butions by writing 



HQ) = Ia + SI, 



(C6) 



with 
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Because the remaining contribution 61 is uv-convergent, 
we have extended its uv cutoff Ea to infinity, thereby 
only neglecting insignificant terms that are down by a 
factor of order 0{fil Ea,Q / Ea) < 1. 

Combining the uv contribution Ia{Q) inside Iai3, 
Eq. (C2), and doing the angular integrals we obtain the 
uv contribution 
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where ci and C2 constants were defined in Eqs.(5.45), 
(5.46), and we used three-dimensional spherical averages 
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This confirms the uv contribution to /^^[B] used in the 
main text, Eqs.(7.18) and (7.19). 

The value of the second contribution, 5I{Q) in 
Eq. (C6) critically depends on the sign of fx. In the EEC 
regime of /it < 0, /t = — 1 and the integral in 61 (Q) is 
convergent everywhere, making only a strongly subdom- 
inant 0{fi/ EatQ / Ea) <C 1 contribution to Ia{Q), that 
can be safely neglected. 

In contrast, in the BCS regime of /x > and Q <C 1, 
the integral in 61 (Q), while uv convergent, makes a large 
contribution to Ia{Q) that is, in fact, logarithmically di- 
vergent with a vanishing Q. This large contribution arises 



from a region around e = 1, physically corresponding to 

low-energy excitations near the Fermi surface. 

We focus on the integration above and below the 
Fermi-surface, 61 = 6{^){6I~ + 61^), with 

6I~ = I dig{e) = / dxg{l - x), 
Jo Jo 

= 6r + J- (C12) 

rOO pi pOO 

61+ = / deg{e) = / dxg{l + x) + / dxg{l + x), 
Ji Jo Ji 

= 6i+ + J+, (C13) 

and the integrand g{i) given by Eq. (C9) with jj, = 1. 
Above we separated the dominant logarithmic contribu- 
tion 61"^ out of 61"^, 



61+ (Q) 



f 

Jo 



dx 



[a;2 + g(l±a;)]i/2^ 



(C14) 



= ln^+/-(Q) 



68 



with 
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and the subdominant in small Q contributions 
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that arc finite for Q ^ 0. Taylor-expanding these sub- 
dominant contributions to lowest order in Q and combin- 
ing everything together, we obtain: 



SI{Q) = 0{fi)\n 



64e-ifi/3 

Q 



■e{n)0{QlnQ). (C20) 



Combining this with Ia{Q) inside lajsiQ), Eq. (C2), we 
obtain = 7^^ - 57„^ 6{fi) fif^y2mV^/7r^, where 



SI. 



al3 



I 



kakf3 In (u • k)2 (v . k)2 , (C21) 



and 



1 



B = u + iv = 



B. 



(C22) 



The spherical average in Eq. (C21) is easiest to evaluate 
in the transverse gauge u • v = 0, taking u and v to be 
the kx = ku and ky = ky axes. This reduces it to 



SI. 



al3 



kak/3 In u^kl^ + v'^kl 



= dct> 
Jo Jo 



dO sin Okakfj x 



(C23) 
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In [sin^ 6{u'^ cos^ (j) + v'^ sin^ 0)] , 

= A{u,v)Sal3 + C{u,v)UaUl3 + C{v,u)VaV0, (C25) 

where in the last line we took advantage of the general 
tensor form and u ^ v symmetry of SIap{u, v) to reduce 
its computation to two functions A{u,v) and C{u,v), 
that can be obtained by calculating any two combina- 
tions of components of SIai3{u,v). With this, we obtain 
in the transverse gauge 



SI. 



af3 



I- 



,-4/3 



{u + v) 



(C26) 



-UaUp 



•iu(u + v) " 3v{u + v 
which when combined with Eq. (CIO) gives the saddle- 
point equation Eq. (7.20) used in Sec.VIIB 1. 

By integrating this saddle-point equation, 
^^21 = (e„ _ 2^)B„ - ^7„^[B]B;3 = 0, (C27) 



over Boi we can also obtain the ground-state energy den- 
sity £gs[B]. Utilizing Eq. (C21) to integrate the last 
term we find 



eGs[B] = (wo-2M)(l + C2)|B|2 + 7pC2 



Sep 
5n 



(B-BY 



IB-BI' 



+e(M)37pMyj / ^(B • k)(B . k) (in 



ao(B • k)(B • k) 



(u;o-2A<)(l + C2)|B|2+7pC2^ 



(B.B)Vi|B.B|2 



- l) , (C28) 



a, 13 



(C29) 



which gives the result quoted in Eq. (7.27) of Sec.VIIB 1 
and used to study phase behavior of a p-wave resonant 
Fermi gas. 



2. Finite temperature 

Above analysis can be easily extended to finite tem- 
perature, by analyzing 



7'^) - o2 

— yp 



d^k fcafc/j tanh (^) 
(2^ ' 



(C30) 



69 



9l 



2m\^,\Nm{i<fl + i<:^l), 



where 
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(C31) 



(C32) 



r 



and we have safely Taylor-expanded in Q since Fermi- these together we find 
surface divergences are cutoff by finite T. Combining 



J 



AT) 

-'a/3 



ci + 2mc2 + l@{n)gl2miiN{ii) In ^ 



= Ia0 + f{f^)gpmf,N{f,) In ^5^fs, 



- -=mglc2 ((5a/3B • B -h -Bq,B;3 -I- BfjBa) , 



(C33) 



where /^^ is given in Eq. (CIO), above. As an- 5apQ{fi)2jpfi^ fi/ep ln(/i/T). The resulting free-energy 
ticipated, the nonanalytic (Fermi surface, ai) terms density is given by 
have been replaced by 6a0^Q{fj,)gp2miJ,N{ij,)ln{ij,/T) = 



/[B] 

1 + C2 



^(a>„(T)-2/z)|B, 



02 



(B.B)Vi|B.B|^ 



(C34) 



with 



u)a{T) =0Ja-ai ln(/i/T), 



(C35) 



determining by a;a(Tr) = 2At. 
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